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ABSTRACT 


—— em 





A pulsed mass spectrometer is described which was constructed 
for the measurement of rate constants of ion-molecule reactions at 
thermal energies. The mass spectrometer could be operated with 


1 


ion source pressures up to 10 rately The reactions studied were 
condensation and attachment reactions which can only be observed 
at high pressures. 

Reactions of millitorr traces of methane in 3.5 torr of 
krypton were first studied to check the pulsing technique. The 


values of the rate constants for reactions (1) and (2) (both 


1.0x10 ~ cc/molecule sec.) were found to be close to previously 
+ + 
tL) CH ace Set Fob ciheta + CH 
+ 
(2) CH + CH cae, es | seus 


determined values and thus the validity of the pulse technique 
could be established. 

The ionic condensation reactions (3) and (4) were studied 
in the xenon-sensitized ionization of ethylene to gain information 


on the radiolysis of ethylene. The values for the rate constants kn 


Nd 
+ n - 
4+- ——)- 
(3) Ch H CS Hy Dae Ho +d 


2n 2 Scythe bee 


were found to agree with values which had been estimated from results 


of a previous technique in which the reaction times were not well 


known ky = Re ihe ie ka = inaearorae and Ke = rel eed? (cc/molecule 
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sec.). The values kK for reaction (4) were measured as 


k 


+ k 
C pects oC 3 
kee 24 k+2 “2k+3 


+ 


(4) e355 


k, = eGxto. Kk, = Giteston and kon) 0. 27x10) 


2 
- (cc/molecule sec.). 


These had not been previously measured. The rapid decrease of rate 
constant with increasing ionic size was explained in terms of the 

; 
ion structure. The effect of Xe on the stabilization of the 


- 
condensation product C was observed. Ethylene pressures were 


H 
4°8 
varied from 1x107° Lo oUso tore and xenon pressures from*0.5 to 8 torr. 


The rate constants and equilibrium constant for the attachment 


reaction (5), which appears to be of aeronomic interest, were 
(5) O 8, aga 


Measured over the pressure range of.1 to.8+torr at temperatures 
from 298 to 341°K. The reaction was found to have third order 
kinetics. The rate constant for the forward reaction was 
measured as > g0x107 > ict Greens sec. at 298°K with a 
negative activation energy of 2.0 kcal/mole. The reverse rate 
constant was 2. 44x10 > 7 cc/molecule sec. at 298°K. The values 
of the equilibrium constant 0.385 (standard state 1 torr), 
entropy -22.4 cal/mole deg. and enthalpy -9.98 kcal/mol. 

agree very well with previous results obtained under equilibrium 
conditions. 


Attachment reactions in nitrogen were studied at pressures 


from 0.4 to 3.6 torr at temperatures from 297 to 420°K. A forward 
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rate constant for the attachment (6) was measured as 8. 3x10 7” 


Peep leciie- sec, at 298°K with a negative activation energy 


(6) N + aN 


of 2.3 kcals/mole. These values support previous observations. 


The equilibrium constant could not be measured. 
Under identical conditions as those used above, the 


attachment reaction (7) displayed second order kinetics. The 


(7) Noistn =e oN 


-12 
rate constant was measured as 1.5x10 cc/molecule second. No 


temperature dependence was observed. The rate constants and 
equilibrium constants for the oxygen and nitrogen reactions 


were compared by means of theoretical equations. 
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INTRODUCTION 








isd) wihe oPresentrstudy 

The present study was directed towards the investigation 
of the kinetics of ion-molecule reactions in gases at high 
pressures (1-10 torr) and at thermal energies. The temperature 
of the ion source could be varied from 300° to 450° K. >The 
reaction times could be measured directly by means of a pulsing 
technique. These unusual reaction conditions were designed to 
fill a gap in the field of ion molecule reaction investigations. 
In the following few pages the development of the field of ion 
molecule reactions, and the ae for measurements of ion molecule 


reactions under these conditions will be discussed. 


1:2 Development of Ion-Molecule Reaction Studies 

In an ordinary mass spectrometer used for analysis (isotope 
analysis, qualitative and quantitative analysis of vapours) the 
pressure in the ion source and mass analysis region is kept low 
so that the ions can leave the ion source and be mass analyzed 
without the occurrence of collisions with gas molecules. If 
the pressure of the ion source is increased, collisions of the 
ions with gas molecules begins to occur and these can lead to 
ion molecule reactions which change the nature of the primary 
ions. For example J.J. Thomson (1) with the first mass 


spectrograph observed ions of mass to charge ratio m/e = 3 in 
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+ 
hydrogen gas and reasoned that these were due to H This 


3 e 
ion was also found by Dempster who observed that its intensity 
increased with hydrogen pressure (2). The origin of this ion 


was later correctly explained (3,4) as being due to the 


reaction (1:1): 
+ 
Cie) H + H 2*3.4H +H 


aie é : ' : : 
The presence of the H ion obviously interferes with isotopic 


3 
determinations when the abundance of the ion HD* has to be 
measured. Since at that time interest in mass spectrometric 
measurements was mostly centered on isotopic abundance 
measurements efforts towards improvement of the vacuum systems 
soon led to the elimination of the undesired ion molecule 
reactions. In subsequent mass spectrometric research which 
dealt with the measurements of ionization potentials, ionic 
heats of formation and analysis of vapours the occurrence 

of ion molecule reactions was not desired and thus low ion 
source pressures were maintained. Only during the early 1950's 
was it recognized that ion-molecule reactions might play a 
large role in radiation chemistry. With this in mind Talroze 
and Lyubimova (5), Stevenson and Schissler (6,7) started 
systematic studies. This early work caused a renewed 

interest in these reactions, and the field expanded rapidly. 


Many reactions were observed and their rates measured in 
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conventional mass spectrometers at ion source pressures of the 
order tae nae torr. The basis of these measurements will now 
be briefly described. 

A typical mass spectrometer used for gas analysis is 
shown schematically in figure (1:1). Ions are produced in the 
ion source by bombardment of the gas by an electron beam. These 
primary ions which have the masses of the original molecules 
and fragments thereof are pushed out of the ion source by the 
small electric field ES of the repeller. They are then 
accelerated by the potential V and analysed according to their 
mass to charge ratio m/e and detected. When the instrument is 
used for gas analysis the pressure in the ion source and analyser 
is kept very low (less than one torr) so that no ion molecule 
collisions occur. When the ion source pressure is increased 
to 10K torr, ion-molecule reactions become observable since 
some of the primary ions eo suffer reacting collisions in their 
path 2 on their way to the ion exit slit (see Figure 1:1). 
Considering the general reaction (1:2), one can write equation 
Piss). oor 


(1:2) ewe as: Shira) 


the attenuation of the primary ion, where i is the primary ion 


ae a7 anh 


(1:i) I = 
Pp p 


intensity at distance &. r is the original primary ion intensity, 


Q is the experimental reaction cross section and n is the number 
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Figure 1:1 Schematic diagram of a mass spectrometer. 


Some secondary ion formation shown. 





of molecules per cc in the ion source. Since I -I. =I 


fe) 

p p S 
, : - : + 

where I. is the intensity of the secondary ion S_ one can 


rearrange equation (1l:i) to the form: 


treat) I 


At low pressures, where only a few of the primary ions react and 
the ratio becomes 1/1, < 0.05 equation (1:ii) can be simplified 
to tne form (12111) 
1:iii r 
(1:iii) ae 


a 
p 


= Ond 


All measurements by the early workers like Stevenson (6,7,8), 
Hamill (9), Field and Franklin (10,11) were done under low pressure 
conditions where equation (l:iii) holds. 

These experiments were directed towards two different aspects 
of ion molecule reactions. One was to determine the types of 
reactions occurring when different reactants were present, the 
other was the dependence of the reaction cross section on the 
velocity of the primary ion. 

Gioumousis and Stevenson (12) developed the theory of 
velocity dependence for the reaction cross section. This was 
based upon a previous theory of Langevin who proposed that the 
cross section for collision of an ion and a molecule was 


dependent upon the polarizability a of the molecule. The ion 
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and molecule converge with a relative velocity g. It is assumed 
that the ion induces a dipole in the molecule and the resultant 
force causes the trajectories to converge. If they pass within 
a certain distance the ion and molecule orbit around each other. 
The critical distance within which capture still occurs can be 


used to evaluate a capture cross section se which is given in 


equation (l:iv). 


(l:iv) Oo a Te a 


Where uw is the reduced mass the other symbols having been 
defined previously. 

The primary ion in the ion source of a mass spectrometer 
experiences continuous acceleration by the repeller field of 
strength Ee In order to relate the experimentally observed 
cross section Q (equation 1l:iii) to the Langevin equation (1l:iv) 
Gioumonsis integrated the variable cross section ae over the 
total path & of the ion. Since the velocity of the ion due 
to the repeller field is very much larger than the velocity of 
the molecules, it is set equal to the relative velocity g. This 


averaging process leads to the relationship given in equation (l:v). 


(1:v) QO = 2Te a 2m, 
u Fs 


whe re ve is the mass of the primary ion. Gioumousis in the above 


treatment assumed that reaction occurs at each capturing ion 
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molecule collision, i.e. that the reaction efficiency is unity. 
Using the Langevin cross section dependence (l:iv) Gioumousis 
also derived an equation for the rate constant of an ion molecule 
reaction where both the ion and molecule moved with thermal 


(Maxwell) velocities. This relationship is given in equation (l:vi). 


(1: vi) k = 27e / O 


u 


Comparing (l:vi) and (1l:v) one sees that the thermal rate constant 
k can be evaluated from the value of 0, experimentally determined 


in the presence of a repeller field, by the relationship (1l:vii). 


Cievis) k =/eEF 2 
/ x 

2m 2 
Pp 


It was shown for many reactions (e.g. 1:1) that the 





experimentally determined Q followed the 1/ ES dependence 
predicted by equation (l:v) and that the absolute magnitude 
corresponded very closely to that calculated by the same equation 
(8,13). The dependence on the reduced mass was also observed. 
Equation (l:v) does HouEveaten a temperature dependence for Q. 
Indeed, it was eliminated by the assumption that the relative 
velocity g is only the velocity of the ion. The values of Q 
were found to be independent of temperature. The large values 
of 0 {> 107166547 and up to 10 Fis ines larger than those of 
equivalent reactions of neutrals) and the lack of a temperature 
dependence, lead to the conclusion that the observed ion- 


molecule reactions were non-endothermic and had activation 
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energies close to zero (14). Thermal rate constants had large 
values Mow cc/molecule sec.) when calculated from (1l:vii) 
and these corresponded with the predicted values (l:vi). 

The validity of equation (l:vii) for thermal rate 
constants depends upon the assumed correctness of equations (1: v) 
and (l:vi). However, many reactions were phanamtoe which 0 did 
not obey the 1/ EX dependence of (1:v). It was pointed out by 
Stevenson (14) that the probability that a collision would lead 
to reaction may not be unity or constant. For complex reactions, 
the probability was not independent of the internal coordinates 
of the reactant ion or of the relative velocity g of the colliding 
pair. It was found to decrease with increasing velocity. 

Thus equation (l:vii) and experiments done in the presence 
of a repeller field could not be used to determine rate constants 
for reactions occurring at thermal energies. Since ion molecule 
reactions occur at thermal energies in a large number of chemical 
systems, such as radiation chemistry, the ionosphere,- flames and 
electrical discharges, it would be desirable to determine k under 
conditions where ions and molecules are in thermal equilibrium. 


The pulse technique was developed for this purpose. 


183 Thermal Rate Constants of Ion Molecule Reactions at Low 
Pressures with Pulsed Instruments 


The pulse technique for observing reactions at thermal energies 
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was first developed by Talroze (16,17). The ion source is pulsed 
as follows. The electron beam is switched on as a pulse for a 


very short time and ions are formed. After this there is a delay ta 


during which equipotential conditions are maintained in the ion 
source and the ions may react under thermal conditions. The 
repeller potential is then applied in a pulse te repelling the 


ions. By varying the time t, and measuring cae as a function 


d 


of ta it is possible to find the rate constant k if the concentration 


of gas n is known. The equation given by Talroze (17) is (1:ix) 


(11355) I. 
ea knt, + E(t. 
p 
where E(t, ee, is a function of the pulse widths and is constant 
, 


ity! 

for constant t, and tL. This method is similar to most kinetic 
investigations in which the degree of conversion is determined as 
a function of elapsed time. The pulse technique has since been 
used by other workers notably Hand and Weysenhoff (18), Ryan and 
Futrell (19), and Harrison et al (20,21) to measure rate constants 


at thermal energies. 


1:4 Limitations of Low Pressure Studies 

The study of ion-molecule reactions by the two low pressure 
techniques, the earlier method of continuous repulsion, and the 
later pulsed method, had certain limitations. Since the pressure 


-4 
was low (around 10 torr) and the reaction time was short (about 


=G } : —2 
10 seconds) only reactions with large rate constants (10 to 


-10 
10 cc/molecule sec.) could be observed. The concentration of 
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secondary ions was usually more than Tone times smaller than 
that of the primary ions and thus further reactions could not be 
observed. However much useful information was obtained from 
observation of these fast reactions. It became possible to 


classify the reaction types as follows. 


s 


cies. Types of Reactions Observed from Low Pressure 


Experiments and the Need for High Pressures 

It is widely. accepted now that low energy ion molecule 
Pe nijens occur through the formation of a short-lived collision 
complex. It has been suggested that the transition complex takes 
one of two peraits (22). One is a loose complex in which the ion- 
molecule system is held together by polarization forces. In 
the second the complex is considered to be covalently bonded. 

For example reaction (1:3) involved a loose complex. It 


was shown (23) that the hydrogen atoms are not mixed, that the 


+ ; + 
ee mgr en 
(bs53) CH, + CD, se? CH, CD, 3 CD, Bid CH, 


product is cD, "H rather than CH,D, ch would result from a 


complex of the form (Cs, Dalen 


The more strongly bonded complex is formed when it has 
the stoichiometry which corresponds to that of a stable ionic 


structure. In this case rearrangement does occur as in (1:4). 


+ 4 om Of 
(1:4) ect HS ses Bi ce 
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: a 


fats 


» 
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+ + a s 
33 and te are Similar 


The intensities of the products C 
to those of the fragment ions produced by electron impact 
ionization of butene molecules (11,24). 

The types of ion molecule reactions are classified below 


more with respect to the product formed than to the type of 


complex involved. 


i Simple Charge Transfer 
Superficially in charge transfer the ion abstracts an 
electron from the molecule during the collision and a new ion 


and molecule are formed as in reaction (1:5). 


(125) xxt (2p 


+ 
)} AeecH sa eI3ChH + Kr 
3/2 


4 4 

The requirement for charge exchange to occur at thermal energies 

is that the recombination energy of the ion be greater than the 

ionization potential of the neutral molecule. (In reaction (1:5) 
+ ak2 He 


RE” Kr P = 14,0 eV and L.(CH 
( 3/y) é aA ; 


in energy -AH is 0.88 eV and reaction occurs readily). 


= 13.12 eV. The difference 


When the ion and molecule are the same species, AH = O and 
the reaction is called resonant charge transfer (25). Occasionally 
resonant conditions are imposed on reactions of different species (1:6), 


+ + 
. ul -)> 
(1:6) N, + N, N. + N, 


especially when the numbers of degrees of freedom of the ion and 


molecule are small. This was observed with xenon and ethylene (1:7) 


(26). 
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, (825) asta sya sagrcarit 


De. 


+ + + 
: r, H > ke C aoe OS EY 
Ci) “Ke + CoH, (XeC,, >) Xe CoH, AH 04 e 
ua Charge Transfer Induced Dissociation 
When the change in enthalpy AH is large because of large 


differences in ionization potentials, the new ion has a large 


excess energy and may dissociate (27) as in reaction (1:8). 


. 
‘ 


ser + 
: K > 
(1:5) papas Bae? - CH, CH, crete ees 


+ 
In this case the recombination energy of Kr (14.67 eV) is 1.55 eV 


more than the ionization potential of CH”. This large excitation 
energy causes a C-H bond to break. The rule for this to take 
place is that the recombination energy has to be larger than the 
appearance potential of the fragmentation. In this case the 
appearance potential equals the bond dissociation energy D(CH ,-H) 
plus the ionization potential of CH. 
By a suitable choice of the reactant ion the extent of 
fragmentation of the production may be specified. This is very 
useful if one wants to form only one or two types of ions for 
further reaction. .The choice of ions is limited by the available 


recombination energies and appearance potentials of the ions and 


molecules. 


re AY Abstlacuszon-o1-.Transter of Atoms or Atomic Ions 








ee ee re 


A large majority of ion molecule reactions involve the 


abstraction by. the ion of an ionic or neutral fragment from the 
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molecule, or the transfer by the ion of such a fragment to the 

molecule. Reactions in which the fragment is hydrogen have 

been categorized as shown below. Reactions in which the 

fragment is large are treated as condensation reactions. 
Ion-molecule reactions of organic ions can be classified 

according to the electronic structure of the ion (28). Generally 

organic ions have zero, one or two electron vacancies in the 


valence shell. This can be illustrated by the methane ions. 


O vacancy 1 vacancy 2 vacancies 


a Proton Donation 
Ions which have no vacancies would be expected to 


react by donating a proton (28). 


+ + 
(1:9) CH, +ND, > CH, + ND.H 


These reactions have long been postulated in organic chemistry. 


b Hydrogen Abstraction 
Ions in which there is one vacancy and thus have an 


odd electron number are radical ions. Two reactions of radical 
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ions are hydrogen abstractions and addition to double bonds. 
Addition will be considered later as a condensation reaction. 


Hydrogen abstraction is shown in reactions (1:10) and (1:11) (29). 


= + 
: > H + 
oa a CPA ar CoH, CD, CoH 


+ + 
° > 
Cif.) CD, — NH, CDH + NH, 


It was observed that when proton donation could be distinguished 
from hydrogen abstraction the rate constant of the former was 
about one order magnitude larger than that of the latter. This 


-++- 
is probably due to the extremely electrophilic nature of H . 


c Hydride Ton Abstraction 
Ions in which there are two vacancies in the valence 
orbital tend to abstract H from saturated molecules (28). 


+ + 
Examples of such ions are CH. and CoH. (30). 


+ + 
. — = 
Cle 2) C3H, + neo CoH,» CoHyy + C.H. 


These ions can also react in condensation reactions. 


Using the above criterion of an ion having zero or one 
vacancy one cannot always predict whether proton donation will 
predominate over hydrogen abstraction or vice versa. Many of 
these reactions appear to involve a loosely bound complex and 
the most stable product ion results. Both reactions (1:13) 


and (1:14) give the same product (22), the former by hydrogen 
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+ 
C213) Dron i> CH, ance BD 


+ 
Oo. ($208 
2 2 3 


2 


+ + 
: > 
(1:14) CH) “ D0 H D0 + CH, 


abstraction and the latter by proton donation. One might expect 
: + m tea” P 
hydrogen abstraction by D0 and by CH, Since both ions have 


one electron vacancy. 


d Atom-Ion Interchange Reactions 
Certain reactions between inorganic diatomic gases 
have been labelled atom-ion interchange reactions and the exact 


mechanism left unspecified (17). In these reactions larger elements 


+ + 
Ca) OS + No = NO. +.8 


than hydrogen are transferred. 


iv Condensation Reactions 
The term condensation covers all those reactions in which 
covalent bonds (C-C and C-S) bind the complex and considerable 


rearrangement takes place so that units larger than atoms or 


atomic ions are transferred. The general reaction is given by (1:16) 
+ +% 
(1:16) AB + CD > (ABCD) > products 
ae , : : . : 
where AB is anion with one or two vacancies in the valence orbital. 


This group of reactions is probably the largest involving organic 


ions and molecules. Examples of these reactions are (1:4) in 
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, 1; . ‘ 
which CoH, has one electron vacancy and (1:17) in which CH,” 


has two electron vacancies. 


+ +%* + 
(a3 17) CH, + CH, - (C,H) - CoH. + H, 


All the above types of reactions were observed at low 
pressures. One expects, however, that as the pressure is 
increased, these reactions will be modified and new reactions 


will appear. 


At high pressures the intermediate complex of a condensation 


reaction may be stabilized by collision with the gas molecules. 
The complex has a large excitation energy due to bond formation. 
This energy is distributed in its vibrational and rotational 
modes, and is transformed into translational kinetic energy by 
collision with the gas molecules. Thus the complex becomes 
thermalized. The product of reaction (1:4) is observed to 


5 aie ie : 
Change at high pressures as (CHa) is. stabilized as in.(1:18). 


+% + 
(1:18) (plied wun Teale chen male aha 


Since at high pressures, the concentration of secondary 
ions becomes large, higher order reactions (consecutive reactions) 
are observed. The products of (1:4) or (1:18) can react further: 


+ +. 4% 
(1:19) CoH, UE CoH, =: (CoH, ) 


+% 
(1:20) rTM aya Ct (CH 5) 


4°8 24 doy". 


Ion-molecule reactions which have low reaction rates become 
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observable at higher pressures. A new type of reaction, which 
involves three body collisions and thus has a low rate, is the 


attachment reaction. 


v Attachment Reactions 

Reactions in which the bonding is mainly ionic, in which 
no rearrangement occurs and which involve the reactions of 
inorganic ions to form complexes bonded by covalent or ionic 
bonds, will be called attachment or clustering reactions. 
Examples of attachment reactions involving only ion-dipole 


: ; SF 
Torcesjare, the, clustering weactions, of H,O. around K., (31). 


2 
= M + 
2 —_—— 
ie ad) K + HAO a enuiley 
+ M + 
(1:22) K (HO) + HO weal OS (HAO), 


The reactions exhibit third order kinetics and are only observed 
at high pressures as the excess energy caused by bond formation 
must be removed by a third body to stabilize the complex. An 
example of a reaction in which the complex bonding may be 
partially covalent and partially ionic in character is the 

three body attachment of oxygen. 


+ + 
€)s2 3) 0, hin2O;, =) 0, he Os 


1:6 Previous Experimental Work at High Pressures 
Instrumental Developments 


Many workers recognized the desirability of obtaining results 
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at high pressures and mass spectrometers were modified accordingly. 
Field and Franklin and their co-workers were among the first to 
do so. The ion source of their mass spectrometer was modified 
to decrease the rate of escape of gas by making the electron 
entrance aperture and the ion exit slit smaller. The ion 
source pressure could be increased without increasing the 
analyzer pressure so that no collisions occurred outside the 
ion source. Thus Field and co-workers were able to observe 
high order reactions in methane (32) and in ethylene (33) at 
ee of 0.3 torr. This apparatus was further improved 
and the ions produced in methane were observed at pressures up 
to 2 torr (34,35). These ions could then be compared with 
reactions proposed for conventional radiolysis of methane. 

Disadvantages arise from the use of an electron beam 
in such apparatus. The electron beam may be intense, but the 
electrons are severely scattered and attenuated by the high 
pressure gas. The scatter can be decreased by collimating the 
beam with a magnetic field. However then, the beam, though 
collimated, travels in a distorted path and the beam to ion 
exit slit 1 becomes uncertain. 

The problem of attenuation can be overcome by use of 
more penetrating radiation than electrons. Rudolph and Melton 
used a radiation from po708 ad source of ionization. The 


polonium was painted onto the repeller plate so that the 


fl Ses a7. 
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whole ion-source volume was irradiated. With this apparatus 

they were able to study high order reactions in ethylene (36) 

and acetylene (37) and observe charge exchange reactions 

between rare gas ions and organic molecules such as methane (38) 
and cyanogen. The problem of beam collimation is still present 
with this type of ion source as ions are formed at varying 
distances from the ion exit slit. Since the path lencths of 
primary and secondary ions are not known quantitative measurements 
are precluded. 

Wexler overcame the problem of penetration and collimation 
by using a beam of 2 m.e.v. protons from a Van de Graaf generator 
as the toutzanaesource (39). At pressures of 1.25 torr, high 
order consecutive ion-molecule reactions in ethylene were 
studied (e.g. reactions (1:19) and (1:20)). The ion Gas was 
observed, which confirmed that at such pressures, collisional 
stabilization (1:18) took place. Similar polymerization 
reactions occurred in acetylene. 

One common failure of these three types of instruments 
was that all used a repeller field to extract the ions. 
Consequently, even though high order reactions could be 
observed, the ions had energies larger than thermal. 

With the idea of simulating conventional gas radiolysis 
Kebarle and Godbole (40) built a mass spectrometer with a field 


free ion source which could be operated at pressures up to 200 torr. 
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The experiments of previous workers on the ionic polymerization 
of ethylene were extended to a pressure of 40 torr (41,42). 
However in this case the reactions occurred at thermal energies 
and one would expect the results to be applicable to the 
radiolysis of ethylene. 

Attachment reactions involving ion-dipole interactions 
were also studied with this instrument (43). Such reactions 


were of the type shown in (1:24). 


+ + 
(1:24) eee oer On (HO). 


The instrument was also used to study attachment reactions of 


rare gas ions (44) e.g. (1:25): 
ee) ar’ + i2Zkr > Ban + Ar 


Tres there is no repeller field in this type of ion 
source, one relies on mass flow and diffusion to remove ions. 

As with Melton's apparatus, the reactions times are poorly 
defined. 

The time dependence of a reaction may be observed by means 
of a pulsed stationary afterglow. This technique was used by 
Fite et al (45) and Sayers and Smith (46). In this method 
a gas mixture is subjected to a short pulse of excitation and 
ions are formed in concentrations proportional to the 
constituent concentrations, and react with the neutrals. Ions 


emerging from an aperture in the side of the reaction chamber 
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are analysed by a mass spectrometer. 

The variation of each ion's intensity with time is 
recorded by an oscilloscope attached to the mass spectrometer. 
The ion-molecule reaction rate constants can be calculated 
directly from the time change of intensity of an ion and the 
neutral concentration. Unfortunately since the excitation 
pulse can excite and dissociate the neutral gas, the ions 
May not be reacting with ground state molecules. This technique 
was used to obtain the first values of aeronomic reactions. 
Sayers and Smith (46) observed the reaction of o* ions in 


oxygen and obtained a rate constant which was 


+ + 
: * 4 
(a820) et 0. O 0. 


temperature dependent. Fite et al (45) observed the attachment 


ze 
of N. to N. and calculated a value for the rate constant. 


(inet) N +N 7 iN 


The technique of a flowing afterglow was developed by 
Ferguson, Fehsenfeld and Schmeltekopf (47) to study ion-molecule 
reactions at thermal energies. In this technique, helium which 
has been ionized or excited by a microwave discharge flows 
through a tube at a pressure of about .5 torr. Reactant gases 
are injected downstream from the discharge. The He ions and 
excited He atoms ionize the reactant gas molecules. A second 


gas whose reactions with the primary ions are to be studied 
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is injected further downstream. The ions are sampled mass 
spectrometrically at the end of the flow tube. - The reaction 

time can be calculated from the gas flow rate and path length. 
Since the reactant gas is injected below the discharge, most of the 
neutral molecules will be in the ground state. Ferguson et al 
developed the apparatus so that the reactant gases could be excited 
by a second discharge, just before they entered the stream, and 

so such unstable species as atoms could be added to the reaction 

in a very selective manner. They were able to study many ion-atom 
interchange reactions of aeronomic interest. These reactions 

had not been successfully investigated before. For example, 


the variation of the rate of reaction of (1:28) with the 


(1:28) , eles N, > NO’ +N 


vibrational temperature of N. was measured (48). 


2 
Because the pressure of the carrier gas He is low 

~ .5 torr, and the partial pressures of the reactant gases 
even lower, only reactions with large rate constants 

=9 sagt . 
(10 co. 10 cc/molecule sec.) may be studied. Fortunately 
most ion-atom interchange reactions have large reaction rates 
(49). The problems of mixing of the gas with the flow stream 
and the effect of ambipolar diffusion coefficients and ion 
detection efficiencies limit the accuracy of the rate constants 


to within a factor of two or three. However, this is the 


case for most reactions studied at high pressures. 
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1:7 The Development of the Present Instrument 

The major limitation of the previous high pressure mass 
spectrometer built in this laboratory for the study of thermal 
ion-molecule reactions, was the uncertainty of the reaction time. 
Since the ions were removed by mass Flow and diffusion, the time 
an. ion remained in the ion source could only be estimated. The 
present instrument was gestghed so that reaction times could be 
measured accurately. 

High energy electrons were chosen instead of a particles 

Vitcawrr Catt 
for the source of ionization as an electron beam if conveniently 
pulsed. The energy used (4000v) provided a reasonable penetrating 
power. The pulse technique was similar to those used previously. 
Ions were formed by a short pulse of electrons. At a certain 
time later the ions leaving the ion source were admitted to the 
Mass analyser by means of a second pulse and thus the time 
history of each ion intensity could be obtained by varying 
the time. 

The ion-molecule reactions studied by this technique 
are described in chapters 4, 5 and 6. Since the actual reactions 
are not all -related,the relevant background and survey of 
previous results will be given at the beginning of each chapter. 

General reviews of ion-molecule reactions are given in the 
following references. Low pressure results are discussed in 


(13,15,16 and 22). Some of the recent reviews on high pressure mass 


spectrometry and thermal ion-molecule reactions are references 
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(39,42,49,50 and 51). The role of ion-molecule reactions 


in radiation chemistry has been reviewed by Ausloos (52). 
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2 EXPERIMENTAL 
2:1 Requirements and Design of the Present Instrument 

This instrument was designed with several requirements in 
mind. We wanted to operate the ion source at pressures of several 
torr and this required a large capacity pumping system to keep the 
mass analyser at low pressure. 

A convenient source of ionizing radiation was required. A 
high energy beam of electrons was chosen as it was easy to pulse 
and reasonably safe to operate. Its use does not have the problems 
associated with the use of radioactive particle sources such as 
Polonium 210 which gradually spreads in vacuum and leads to 
contamination of the apparatus. An electron beam is conveniently 
pulsed by a pair of electrostatic deflection plates. 

Finally a mass analyser was required. A quadrupole mass 
spectrometer was chosen as it has a high transmission and thus 
greater sensitivity than a magnetic instrument. The absence of 
a magnetic field makes ion beam collimation easier. The quadrupole 
unit also has a relatively low cost and is compact so that it is 
easily mounted in a vacuum chamber. The main components of the 


instrument are shown in figure (2:1). 


Aza the Ion Source 
The high pressure ion source is shown in figure (2:2). It 
was made from a 1s inch long, 3 inch diameter cylinder of non- 


magnetic stainless steel (type 304), through which a 5/8 inch 
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Figure 2:1 


Electron filament. 

Filament support plate. 

Electron extraction electrodes. 

Grounded cone (electron acceleration electrode). 
Collimating electrode. 

Solenoid coil providing ariet magnetic field for 
electron focussing. 

T.V. yoke for horizontal and vertical motion of the 
electron beam. 

Deflection plates for electron pulsing. 
Fluorescent beam focussing screen. 

Electron entrance cone carrying slit leak. 

Ion exit cone carrying circular leak. 

Ion source with heaters. 

Electron trap. 

Cylindrical shielding screen. 

Ion acceleration and focussing electrodes. 
Quadrupole mass analyser. 

Secondary electron multiplier. 


Gas inlet tubes. 
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Section 
Scale 1:1 


Figure 2:2 Top view and section of the high pressure ion source, 
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channel was bored perpendicular to the axis to provide a gas 
chamber. On the left-hand side was mounted a cylindrical electron 
entrance cone [11] which was sealed to the body of the ion source 
by a Viton "A" Q-ring and held in place by a threaded steel ring. 
The right-hand end of the chamber was plugged by a steel insert 
which contained a Kovar glass high resistance electrical feed 
through [14] on which was mounted a’ Faraday cup electron trap. 
The trap was made from two concentric cylinders of .001" 
Tantalum spot welded. The trap was made from two concentric 
cylinders so that most of the secondary electrons emitted by 
impact of the primary electron beam, would not be able to 

escape from the trap. 

Two gas inlet tubes [23,24] were sealed into the ion 
source (by an O-ring seal) so that gas could be flowed through. 
In most of these experiments, however, the exit tube [24] was 
sealed and gas flow was not used. 

The electron beam entered the ion source through a small 
slit mounted on the cone [11]. The slit was made by welding 
stainless steel razor blades to a circular disc of .007" thick 
Tantalum which contained a 1.7 mm (1/16") hole. The final size 
of the slit was 1.70 mm x 0.040 mm. The Tantalum disc was 
welded to a narrow ring of platinum which in turn was soft 
soldered to the stainless steel electron entrance cone [ll]. 

The ions left the ion source through a small hole 


mounted in a cone [12] which was sealed by an O-ring into 
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the bottom of the source. The ion exit hole was a 0.120 mm 
diameter hole drilled with a fine needle point in a 0.001" 
Tantalum aie! This in turn was spot welded onto a platinum 
washer which was soft soldered onto the top of the ion exit 
cone. The ion exit slit cone was held in place by a threaded 
ring. The inside of the ring was bevelled so that the high vacuum 
Side of the ion exit cone consisted of a smooth conical surface 
with a solid angle of 40°. 

Below the ion source was mounted an electrostatic shield 
[15] which provided a boundary for the electric field between 
the ion source and the ion acceleration electrodes (see figure 
2:1). The shield was a thin wall cylinder slotted longitudinally 
so that the'gas pumping speed was as high as possible, without 
deforming the electric field. 

In front of the ion source entrance slit there was placed 
an electron beam focus plate [10]. This plate contained a 1/4" 
hole axial with the electron entrance cone. The plate was 
coated with a mixture of phosphor, sodium silicate, and sodium 
chloride. The electron beam caused the phosphor to glow 
so the beam could be focussed. The sodium chloride was added 
to make the phosphor coating conductive to prevent surface 
charging. The electron focus plate was extended down in front 
of the electrostatic shield [15] to prevent spurious 


ionization, caused by stray electrons from the electron beam. 
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The ion source was heated by six heaters made from 28 
gauge nichrome wire, spiral wound on ceramic centres and insulated 
from the metal block by quartz tube insulators. Two iron- 
constantan thermocouples [T.C.] were inserted into the ion source. 
One was close to the heaters so that it could be used to control 
the temperature. The second was inserted in a well close to the 
gas chamber so that the gas temperature could be read. It was 
assumed that the gas reached the temperature of the walls of 
the ion source chamber. The thermocouples were mounted in short 
lengths of 1/8" copper tube which fit the wells tightly. The 
thermocouple joint was electrically insulated from the copper 
by thin strips of mica. The electrical insulation was only 
sufficient for a few hundred volts. However it provided 
electrical insulation between the ion source and the heater 
temperature controller. 

The ion source was mounted on the vacuum chamber port 
flange by a cylindrical steel support [13a] which was electrically 
insulated by a disc of Araldite epoxy resin. The Araldite 
plastic could withstand a temperature of approximately 100°C 
before plastic deformation took place. The steel support was 
slotted to decrease heat conduction. The Araldite insulator 
was bolted directly to the flange which remained at room 
temperature and thus its temperature remained approximately 


the same. Nine electrical feed throughs were provided on the 
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flange. The gas lead-in tubes [23,24] were heated by nichrome 
wire heaters wrapped around the glass. An iron constantan 


thermocouple was attached to the inlet tube [23] by glass tape. 


2:3 Temperature Control and Measurement 

About 40 watts were required to heat honor source to its 
maximum allowable temperature of 150°C. The temperature was 
controlled by manual setting of the heater voltage. It was found 
that this steady state method gave a much more constant temperature 
than when a temperature controller (A.P.I. Instrument Co. Model 905B) 
was used. The heater voltage was supplied by an auto transformer 
erie) connected to a constant voltage transformer. 

The thermo-electric e.m.f. of the thermocouple close to the 
gas chamber was measured by a potentiometer. Ice water was used 
for the reference thermocouple. The temperature could be read to 


within one degree centigrade. 


2:4 Vacuum Chamber 

The main vacuum chamber was an 8 inch diameter stainless 
steel tube (figure 2:1) about 12 inches long. This was pumped by 
a National Research Corporation (HS6-1500) 6 inch, 1500 liters per 
second, diffusion pump through a water-cooled optical baffle. The 
pumping speed at the top of the baffle was quoted to be 7502%/second. 


One can calculate the pumping speed at the ion source (at the end 
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32 
of the chamber) if the conductance of the lead is known by 


equation (2:i) from Dushman (53, p. 130) 


(2:1) 


where S is the pumping speed at the ion source and a3 is the 

speed at the pump (750%/sec.). F, the consuctance of the lead, 

may be obtained from a table as in reference (53, P.97) Of g(>4, Pibols 
In this case F = 1000%/sec. and thus S = 4002/sec. 

At the top of the pumping lead four 4 inch pipes were welded 
into the § inch lead, mutually perpendicular, to provide mounting 
flanges (see figure 2:1). The one designed to hold the quadrupole 
probe was made "T" shaped. A second high capacity diffusion pump 
(NRC VHS-4) was mounted on this flange so that the quadrupole was 
differentially pumped. An ion gauge tube was mounted on the flange 
which is shown unused in figure (2:1). 

The leakage of gas from the ion source was very high. Gas 
leaked through the ion exit hole and the electron entrance slit. 
Since the latter was larger most of the gas escaped through it. 

The conductance Fo of an aperture under molecular flow conditions 


is given by (2:ii) (53), where U is the average velocity of the 
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U inmem/sec. Fo is in units of cc/sec. The two apertures had 
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—A Fades. ; - 
areas of 1.1x10 cm (ion exit leak) and 6.8x10 rs ake (electron 


entrance slit). At 300°K the velocity of N, gas is 4.7710" cm/sec. 


2 
The calculated conductance of both leaks would be 9.3 cc/sec. The 
measured conductance was 8.5 cc/sec. This indicates that the 
flow out of the ion source was essentially molecular. 

The pumping speed at the ion source (aheuslay at the ion 


gauge which was close to the ion source) can be calculated by use 
of equation (2:iii), where Fo is the total conductance of the leak, 
Or pam ia op | F 
P the source pressure and he the ion gauge pressure. At 300°K 
Fo was measured as 8.5 for Noe We was fe x18? torr when P. was 
1.1 torr. If these values are substituted into equation (2:iii), 
then S = 560/sec. which is close to the value of 400%/sec. 
calculated before. Since the ion gauge is closer to the pump 
one would expect that the pumping speed obtained on the basis 
of the ion gauge reading will be somewhat larger than the actual 


value at the ion source. 


2:5 Gas Inlet System 

The gas inlet system is shown in figure (2:3). The construction 
materials were stainless steel and glass. The vacuum line was 
surrounded by an oven which could be heated up to 150°C to bake 


out impurities. The pressure in the ion source and the gas inlet 
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system manifold was measured by means of an Atlas MCT manometer 
which read pressures up to 20 torr. Higher pressures up to 200 
torr were measured by the Wallace Tiernan mechanical diaphragm 
gauge. The gas handling plant could be pumped by either of two 
pumping systems. The rough pump was used to remove large amounts 
of gas. The "clean" pump was used to remove the last traces of 
gas. Since the clean pump did not pass large quantities of gas, 
it had very low back streaming and thus the impurity levels in 
the gas inlet system and ion source could be kept very low. 
Large conductance valves were placed in the line between the 
storage bulb and the manifold and the manifold and the ion 
source so that there would not be any appreciable pressure drop 
between the bulb and the ion source. The diameter of the 
manifold and the tube to the manifold was 1.5 cm. The large 
conductance valves had 3/4 inch diameter cross sections. The 
small valves had about 3/16 inch diameter openings. 

During construction of the inlet system the 4 liter bulb 
was calibrated by weighing the amount of water it could hold. 
The volumes of the other parts of the inlet systems were 
measured by expansion of gas from the 42 bulb. The calibration 


was required so that gas mixtures could be made. 


2:6 The Electron Gun 
The electron gun is shown in figure (2:1). Electrons 


emitted from the filament [1] were accelerated slightly by 
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the field between the filament support plate [2] and the 
extraction electrode [3]. They were then accelerated and 
focussed slightly by the focus electrode [4] and finally 
accelerated to the grounded cone [5]. The potentials on the 
electrodes were as follows. The fittament support plate [2] 
potential was -4000 volts. The centre of the filament also 
had this potential. The extractor [3] potential was -3990 volts. 
The focus [4] potential could be varied from -3500 to -3900 
volts. The gun was designed with a large distance between the 
group af electrodes at high voltage [1,2,3 and 4] and the 
grounded cone [5] so that it could be operated at voltages 
up to 25 kv. The electron beam was focussed by the solenoid 
coil [7] (1000 turns of 22 AWG copper wire with 0.2 to 1.5 
amps generally required for focussing). Horizontal and 
vertical deflection of the beam was obtained from a television 
tube yoke which was also outside the vacuum envelope. The 
vacuum envelope at this point was made from brass so that the 
magnetic fields would not be seriously disturbed. The solenoid 
was controlled by a low voltage (7 volts) gas chromatograph dc 
supply modified so that two outputs were available. | 

The deflection plates are [9] and [9a]. A cross section 
of the tube at this point is shown figure (2:4). Electrostatic 
plots were made by means of conductive paper to find the 
required number of guard plates so that the field would 


be uniform. Three were required. The electrostatic potentials 


— PA: 
soljoae see, 4 ine bid te ats ea ? a wp bos 


a 
*/ 


brie 


Ljeveotamando oes (er lov Ss 


nies CS): wth dues ure qromsdih ‘oat nid 


Lhd 


. —— a 
beta toleeee rept ~ ew 


S doce Drhewel i, Ont ik kad 2s aE 












a) . 4 y ve 


ee 


\ 


: 


“ow vod rs 


Pa Ants. 
TAT ‘emettesto evel amt aa on 
_ _ 


(Heitanoy bag 2.18) en matuao on tt ” 


’ Aa yidmre © if? : oy} Lov vy 
ne net at 
j $q) fr 
5 J 4 - Le & rut , * 
¥ ‘ : 






| blugm, 2 iced te 


+ 2 


aa 


rane oe 
Th Bae bot bi woe-hool See oneal 






a oe 
' pe 7 


t ’ > ae sacl Ti 
Ware ‘ : 5.274 SEF 
+ LE f Sir552 in 
4 ‘ Wy ret of Bod aff 


ty 
aco ls oo MmonOopv ayt+ eee sate 


6 Sd MOURA olugths peut raion 1 eh a Ld 


a? boski2kp vied suai nif? ate LO aes 
5 ; 7 : WS 
vi T) at Gag S08 vor E 


_, 


i P vie © 


shdphisvs stew: ‘adic sub. ae re 


’ 


is ~) 


>" oe ie se 
: ieerens sand et tye oe: ‘evils ip an | 
; By 


re! om 
‘ 
Seth © 
é 
1.4 
aoa f 
d rus a 
~~ 
{ 
é r, (Vs 
aad ry 
HOGS 
‘ 7 rt J 
: . 
r4 
he? 
4 
al » 
j * ape 
a Ps 
i Tet 
Ps ) . 
fs. Ons) 
a 
¢ 
je, 2 er 
” \ 
‘ 
a! 
: 
+ rt 
bed 
P 
oy ‘ 
oe 4m 
oy 
Kiso DER A 
x 
of" 
z 
aw 7 
2 
ferry! 
Preah 
g. 
_—s 
; { 
: 


Bi wie? it ran Paaeied 








0 volts 
guard 


plate 





deflection a é 
plate (9a) vacuum! case 


ERT ESE 


Teflon support 


Figure 2:4 Cross section of deflection plates showing guard plates 


and electric field potential lines. Section (A-A) of 


figure (2:1), Scale 4:1. 
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are shown in figure (2:4). One plate [9a] was kept at 
O potential and the pulse applied to the other. A network of 
resistors (total resistance 10° ohm) was placed across the 
guard plates so that there was a uniform voltage gradient 
across them. The 10° ohm resistance was also used to terminate 
the coaxial cable from the pulse generator so that the voltage 
pulses experienced by the plates were as square .as possible. 
Electrode [6] was a collimation apperture to prevent stray 
electrons from hitting the deflection plates. Electrode [10] 
is the focus plate. The beam could be focussed to aspot 
from 2 to 5 mm in diameter (depending upon the position of 
eViagent used) and deflected about 8 mm from the centre. 

The voltage required to deflect the electron beam is 


found as follows. 





The angle is found from equation (2:iv) (55). 


(2:iv) tan 9 = Z2 Vp 
2y we 


where Z is the length of the plates and ye the electron acceleration 
voltage across the plates and Ve the electron acceleration voltage. 


Since tan 8 is also related to the amount of deflection a and 
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39 
the distance b from the centre of the plates to the focus plate, 


G2: Vv) tan oo = a 
b 


The voltage required to deflect the beam can be calculated from 


fore ee et aos ecm, oOo em, 2 = 15 cp, yi =~2~em; ips 4000 volts, 


(2:vi) V =» 2V + 
Z 


to 
0) 
om 


then ve = 60 volts. The pulse generator output could be varied from 
O to 40’ volts and thus if the beam was centred on the focus plate 
hole 5 mm diameter it could be moved completely clear (about 8 mm). 
Two materials were used for the electron gun filament. The 
first was tungsten ribbon Vv .005" x .030". This was used for the 
experiments involving CH, and Kee CoH 


and Xo The life of the 


4 4 


filament was reduced drastically by oxygen. In this case a thoriated 
iridium filament was used. This filament material was very 
stable in oxygen. It also performed well with nitrogen although 
initially the presence of Ny decreased the emission. 
2:7 Ion Acceleration and Gating 
The ion acceleration electrodes are shown as [16] to [20] 
in figures (2:1) and (2:5). Two sets of electrodes were used. 
The first is shown in figure (2:5a) and the second in (2:5b). 
Ions leaving the high pressure ion source exit [12] were accelerated 


by a negative potential to the quadrupole entrance cone [16]. This 
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cone was mounted upon the existing quadrupole mass spectrometer 
ion source plates [17] to [20]. The ions passed through this 
auxiliary ion source into the quadrupole. 

This auxiliary source was used for low pressure gas analysis. 
[17] is what is normally called the case of an ion source. [18] is 
the electron filament support and reflector (Pierce gun type) 
and [18a] the normal trap. The ion acceleration potential (about 
7 volts) was applied between [17] and [20] the quadrupole entrance 
aperture. Electrode[19] was an ion focussing electrode. 

The electrode potentials are shown in Table 2:1. The ion 
source potential was kept at +7.5 volts as this was the value 
which gave maximum signal on the high mass range of the quadrupole 
analyser. 

For the earlier experiments, the studies on methane krypton, 
ethylene-xenon and oxygen, the electrodes in figure (2:5a) were 
used. In this case the electrodes [19] and the bottom half of 
[17] were pulsed to "gate" the entrance of ions into the mass 
analyzer. When the gate was closed electrodes [17] and [19] were 
at +40 V. The pulse from the generator (-40 V) returned this 
to zero so that the ions could pass through for the duration of 
the pulse. 

Because of the presence of a battery (40 volts) in the 
pulse line it was found that the pulses were not square due to 


the large capacitance. The second group of electrodes shown 
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Potentials on Ion Acceleration Plates 


Electrode 


When auxiliary ion source used. 


.) 


16 
iY, 
18 
19 
20 


When high pressure source used as (2: 
12, eS 
14 
16 
Le 
18 
Vie 
20 


When high pressure source used as (2: 
b2 15 
14 
16 
l6a 


léc 
5 
18 
19 
20 


Electrode Potentials 


Volts 


a. 
Shed) 
-50 
QO grounded 


5a) 
at 
= ot ek, 
A907 So R225 
0 gate open (pulsed) 
+40 gate closed 


0 grounded 


5b) 

nal fs 

+7.5 

-67.5 

-67.5 gate open (pulsed) 
0 gate closed 


Os OC: 7 Ca © 


grounded 
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in figure (2:5b) were designed so that the capacitance of the 
system could be reduced to obtain square pulses and also so 
that the flight time of the ion from the ion source to the 
ion gating position was reduced. 

The ions could be deflected by the pair of plates [16a] 


and [l6b]. These were two half cylinders mounted inside the 


~ 4 


cone. The electrode [16b] was kept at the cone potential (-67 
volts) all other electrodes [l6c] to [20] were kept at ground 
potential. During the interval between pulses, electrode [16a] 
was a. at zero potential. All ions would be deflected and 
discharged on [16b] or [l6c]. A negative pulse (-67 volts) 

was applied to [16a] so that the space inside the cone became 
field free and the ions passed through into the quadrupole 


analyser. 


2:8 Travel Time of the Ion from the Ion Source 
to the "Gate" Electrodes 

It was necessary to know the travel time of an ion from 
the ion source to the point at which the ion "gate" potential was 
applied. This travel time was calculated from the potentials 
of the electrodes. 

A simple model was used at first for the electrodes in 
figure (2:5a). We assumed that the ion accelerated across 


a uniform field between the ion exit cone [12] and the quadrupole 
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entrance cone [16] and then travelled at constant velocity 
through the cone to the ion gate [17,19] where it was stopped. 
It was observed that the two distances were about the same. 

For an ion accelerating in a uniform field the terminal velocity 


UL is given by (2:vii), e/m is the charge to mass ratio and Vis 


(22711) Vv = 2eV 


the potential of the field. Since the average velocity is half 


the final velocity,the time is given by (2:viii) where s is the 


(237511) vel, Peay hs 


distance between the cones. Since the ion travelled at constant 


velocity U.. through the quadrupole cone of length s', this time 


is given by (2:ix). Since the two distances were equal, s' = s, 
(2:ix) th 25 ost. 
Ve 


the total flight time t, is given by (2:x). This formula was 


Fe 


(2:x) ty == 3s / m 
2eV 


used to find the flight times used for the methane-krypton 


results. The flight times are shown below in Table 2:2. 
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TABLE A °2 


Flight Times of Tons Calculated by Equation (2:x) 


for Electrodes Shown in Figure (2:5a). Acceleration voltage = -225 volts 
Ion (m/e) Flight time (micro sec.) Value used 
“3 
Kr (84) 5.0 6 
he , 
CH, (15) ae 53 3 
4 
CH (16) . ak 3 
ae 
CH, (17) 2.45 3 
C,H,” (29) aed 3 
205 ; 
+ 
oO. (32) 320 4 
+ 
0, (64) 5.40 5 





— 


Sees com, .V = 232eVolLts nse ee eee me m= g/6.023x10°" 


A second method was used to calculate the flight time. The 
assumption that the field was uniform was not correct as the field 
did not penetrate the ion exit cone evenly. The cross section of 
the ion source and electrodes was drawn with silver conductive 
paint upon conductance paper. Electrostatic plots were obtained 
for both electrode systems. These are shown on figure (2:5). 

The total potential was divided into ten equal units. The time 
the ion took to travel from one equipotential line n to the 


next (n+l) was calculated. The total time is given by the 
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summation of these times by equation (2:xi), where the voltages 


(2: xi) Cee) ee / ih (Va os J Yard Joe tebeeo 10 
Av 


at an equipotential level n is given by Vn! and distance between 
n and n-l is qd: AV is one tenth, in this case, of the total 


potential V and equals Sena a 


~] 
It was found that since the equipotential lines were more 
widely separated at the beginning of the ion's path, the times 


were slightly larger. The difference was less than 1 micro 


second. The travel times for nitrogen were calculated this way. 


TABI Ee 2S 
Flight times of Nitrogen Ions Calculated by Equation (2:xi) 


Electrode Arrangement Shown in Figure (2:5b) 





Acceleration Voltage = -67 volts 
Ion Flight time (micro sec.) 
Nn" Bue 
+ 
4.5 
me 
+ 
5.4 
Na 
Net 6.4 
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2:9 PulsefCircuits 
tnegorder of the pulse generators is shown in figure (2:6). 

The order of the pulses is shown on figure (2:7). The first 
pulse generator (General Radio, Type No. 1217-C), was used to 
generate the pulse of width At. applied to the electron beam 
deflection plates. The synchronization signal was fed into the 
oscilloscope (Tektronix BBG), and generator number 2 (Datapulse 
Model 100). This was used as a variable delay, to change the 
times t" between the pulses. This delay time was the reaction 


time t plus the flight time t The output of this generator 


&° 
(actually the delayed synchronisation pulse) was used to trigger 
the third generator (Datapulse Model 100) which produced the ion 
gate pulse of width At. . The delay time t" was measured 
visually by means of the oscilloscope as the vernier scales 

on the pulse generators were not accurate. The technique of 
using an ion gate to establish the ion reaction time was chosen 
as it allowed the mass spectrometer to operate continuously. 
Since the signal amplifier ee slow response time, the signal 


indicated on the chart recorder was the integrated signal from 


IMany successive pulses. 


2:10 Quadrupole Mass Analyser and Detection System 
The Measurement of the Transmission of the Analyser 


The quadrupole mass spectrometer was an Electron Associates 


Inc. Model 200. Three mass ranges were provided, Low (1-50), 
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Figure 2:6 Pulse generation circuit. 
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Medium (10 to 150) and High (10 to 500). The quadrupole probe 
was powered by an R.F./D.C. Generator which was part of the unit. 
Ion detection was by means of an electron multiplier from which 
the signal was fed into a Keithley Model 416 Picoammeter, and 
thence into a high speed oscillographic recorder (Sanborn 

4500 series). Since the theory of operation of a quadrupole 
analyser may be required to explain the operation it is given 
below. 

The quadrupole mass spectrometer was developed by Paul 
and co-workers (56,57,58). It consists of four long parallel 
rods (figure 2:8a). present Seri vars cross section of the rods 
should be hyperbolic, but a circular cross section may be used. 
Opposrte pairs are connected together electrically. A dc 
voltage U and an rf voltage Ve cos Wt are applied to the poles 


so that the positive pair of electrodes has a potential 
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+(U + Wee cos Wt) and the negative pair a potential -(U + pe COG OL). 


The electrodes are separated by a distance 2x, and the potential / 
at any point is then given by equation (2:xii), where x is the 


2 
(2:xii) = (5 -y ) WU + “A cos wt) 
2 


na 

fe) 
direction of the positive rods, y is that of the negative rods and 
zis the longitudinal axis. Ions are accelerated into the rf-de 


field along the z axis. The motions of the ions are described 


by equations (2:xiii) to (2:xv). 
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Figure 2:8 (a) Quadrupole rods. The hyperbolic cross section is 


approximated by a circular cross section. 
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(b) Stability diagram of the Mathieu paramaters a and q. 
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(2:xv) is most easily integrated to give dz = constant. This 
dt 

means that the ions travel at constant forward velocity through 

the quadrupole rod structure. Equations (2:xiii) and (2:xiv) 

are solved by use of a Mathieu function. The Mathieu equations 


Gf motion are (2:xvi) and (2:xvii). 


2 


(a-avi) ax + (a+ 2q cos 26) x = 0 
ap? 
oe, 2 
(2:xvii) dy - (a+ 2q cos 20) y = 0 
St ? 


Three dimensionless parameters are defined. 


(2: xVLil) od = wt 


Lax xis) a = S8elJ 


(2:xx) q = 4eV5 





When an ion enters the quadrupole structure its trajectory 


is controlled by these parameters. There is a region of stable 
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oscillations where the amplitude remains less than 25° This 
region is shown on the stability diagram shown in figure (2:8b). 
Ions with values for the parameters outside this region follow 
a trajectory of increasing amplitude and strike either the x 
rods when gq is greater than the x stability limit or the y 
rods when a is greater than the y stability limit. Thus only 
ions with values of a, gq, U, ea W and O will pass through 
and be detected. The number of ions with different masses 
which are detected depends upon the ratio a/q. If a/q is 
small (the dashed line a/q), many ions of similar mass (m), 
Ms m.) will be detected and the resolution is low. If a/q 
is large, as is shown by the scan line (the solid a/q line) 
on figure (2:8b), perhaps only one mass m) will be stable at 
a time and the resolution will be high. Resolution increases 
with the a/q ratio and theoretically becomes infinite at 
a/q = 0.336 (a = 0.23699 and q = 0.70600). The quadrupole is 
operated as a mass spectrometer with a/q just smaller than 
0.336. Mass scanning may be accomplished in two ways, w is 
kept constant while U and tc are varied so that the ratio 
AN is kept constant, or U and ee are kept constant while w 
is varied. The mass spectrometer used in this work used the 
first method. The transmission of the quadrupole is large 
because no ion defining apertures are required. 

The transmission of an ion is dependent upon the Ue 


ratio (UV, = a/2q). Von Zahn et al (59) assumed that q was 
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kept constant, at the apex of the stability diagram (figure 2:8b) 
and that any deviations from a/q = constant, were then dependent 
upon a. He developed an equation (2:xxi) for the observed 


intensities I, and I, of two ions after mass analysis. 


ie 2 
re) 
2: i ~ : a 
(2 2ex1) aT 3 I, aie a,) where aaa 23699 
pee —5 Bs a 
To to (a - a) 
max 2 


ae and "i are the original intensities of the two ions before 
analysis. Thus Sea and ele are the transmission factors 

Ore tire two - Lonss. If the ratio Vis is constant as U and Ne are 
Changed a/q will remain constant (a/q = 2U/V .) and thus one would 
expect constant transmission for all ions. The ratio ve was 
measured for the low mass range of the present instrument. The 
Farie did nee remain constant as U and ve increased. Figure (2:9a) 
shows the variation of U/Ve as Ne was increased. The ordinate is 
expressed as a fraction of ve aan where Ve re. is the highest rf 
voltage for the particular mass range. Since it was difficult 

to measure Na accurately, the measured ratio u/s was normalized 


so that at .05 Ve ; Lae =a /2q = 0.16784. The variation 


max max max 


of intensity with ue was calculated by means of the von Zahn 
equation (2:xxi). The parameter a was calculated from the Give 
ratio with gq assumed to be constant at 0.7060. The intensity 
is shown as a relative intensity (with I at 0.5 Me oe taken as 
unity) in figure (2:9b). One can see that the transmission of 


ions of low and high masses would be less than ions of 
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intermediate mass. The equation (2:xxi) gives a linear 
relationship between transmission and a/q (q constant). 
Brubaker and Tuul (60) reported plots of transmission factors 
versus a/q. These plots show that for their instrument the 
transmission decreased slightly as a/q increased from 0 to 
perhaps 0.330 and then suddenly decreased as the ratio 
approached the maximum volume 0.336. Thus small changes in 

a/q would change the transmission factor by a very large amount. 
This was found to be the case for the present instrument. The 
transmission was very dependent upon small changes in the ratio. 
Brunée et al reported (61) that even at constant ee ratio 

the transmission decreased with increasing m/e and thus the 
equation (2:xvi), used to generate the theoretical transmission 
in figure (2:9b) is probably too simple. 

Brunée et al (61) also described a method by which the 
quadrupole could be used to collect the total ionization. If 
the dc voltage U is switched off, a/q = 0, (a/q = 0 is the ordinate of 
the stability diagram (figure 2:8b)), and there is no discrimination 
of ions. Acertain small rf voltage Ne is required to focus the 
ions. All ions of mass greater than that corresponding to 
the small value of ue are collected. The present instrument 
was modified to permit this method of operation. 

The transmission factors were measured in the following 


manner. The ion source exit leak of .120 mm diameter was 
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replaced by a large leak V 1 mm diameter. A repeller was 
welded to the electron trap so that the ion source could be 
operated as a conventional low pressure source. The only 
difference was that the electron energy was 4000 eV instead 
of 70 eV. The quadrupole entrance cone ([16] of figure 2:5a) 
was connected through a battery (-90 volts) to a sensitive 
amplifier (Keithley 610B electrometer). The ion current 
hitting the cone was measured. It was found that when the 
potential on the gate electrodes [17] and [19] was positive 
(+40 V) so that the gate was closed, the current measured on 
the cone increased. This change in current AI was taken to 
be the current of the ions passing through the cone into ane 
quadrupole analyser when the ion gate was open. The ion 
current at the exit of the analyser was divided by AI to 
give a transmission factor. Only the results in which the 
transmission was constant as AI was changed, by changing the 
ion source pressure, were used. 

The gases used to calibrate the mass spectrometer were 
chosen so that electron impact produced ions of one mass 
or ions which were spread over a narrow mass range. The 
rare gases were used first. The intensity of the major ion 
was corrected for isotope ratios and for contributions due to 
doubly charged species. The cross sections for production 


by 4000 V electrons of multiply charged ions, relative to 
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singly charged ions were taken from tables in Schram (62). 


They were usually less than 10%. The diatomic gases 0O,, N 


Ro 2 


etc. were used next. The intensity of the molecular ion was 
corrected for the contribution of atomic ions to the total 
ion current AI. Finally hydrocarbon gases were used. These 
were chosen such that a group of ions of similar mass were 
produced in large abundances. For example CH, is ionized to 


4 


ctipl and CH,” whose intensities total over 90% of the total 
ionization. Other hydrocarbon ions were used. The relative 
intensities of the major ions were taken to be essentially 
the same as those reported in the API tables for 70 eV electrons 
(Kebarle and Godbole (63) showed that the fragmentation pattern 
remained essentially constant as the electron voltage was 
increased). By this method it was possible to establish 
transmission curves for the low and medium mass ranges. 

The high mass range was calibrated by a slightly different 
technique. The transmission factors for three masses were found 


+ + 
m/e 84 (Kr’), m/e 92 (C.H,-CH,”) and 132 (xe"). A spectrum of 


S 


c,cl. was run using the auxiliary low pressure ion source with 


the electron energy set at 70 eV. This was compared to a 
spectrum obtained by an A.E.I. MS-2 magnetic mass spectrometer. 
The relative factors of the ions of C,Cl, obtained by comparing 


the quadrupole spectrum to the MS-2 spectrum were normalized 


to the factors obtained for m/e 84, 92 and 132. 
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The transmission curves for the three mass ranges are 
shown in figure (2:10). They were plotted in logarithmic 
form for convenience. These curves are remarkably similar to 
the theoretical curve, figure (2:9b). 

The transmission factors were checked later when the 
ion acceleration electrodes shown in figure (2:5b) were installed. 
The ion current was measured on the deflection plate [16a]. The 
ion source potential was +7.4 V and the cone [16] potential 
-67.5 V. The deflection plate [16a] had a potential of -135 volts. 
The current was measured by the Keithley 610 B electrometer. 
Essentially the same results were obtained by this method and 
thus the use of AI, as the current passing through the cone 
was verified. 

The total ionization was collected by means of the 
quadrupole, when the dc voltage U was switched off. The 
transmission of the total ionization of Noy with respect to 
the current measured on the deflection plate [16a] was about 
30%, compared to 15% for N. when the ions were mass analyzed. 

This could be interpreted two ways. The statement that when 

U = 0 all ions are focussed by a small rf voltage Vor may 

not be strictly true. Alternatively there could be a considerable 
loss of ions in the path from the cone [16] to the quadrupole 
entrance aperture [JO]. This is more likely. It would 


mean that the transmission of the quadrupole analyzer itself 


rf 


ru 


. 
. 


| a oes a 
on apsije x note. W204 ‘jes baat rmaos tt 3 acd kee 


sidsiebianes & 7 ney geaitt » esbvkane 
» % . * ‘ 


uy 
Loe 
an 
af 
m 
2 
> 
re | 
r 
a 
oie 
e* 
a 
Fa 
% 
ib 
a 
= 
© 
- 














_ P 
+ 4c *prte videdtamet exe aavtlo sandr 8 =««S 
| ad = ole & es « 


re 
i4e:S) sxopit «evige 
a 
34 «ndy waded Betoed® etsew evasos2 lett 0 
w Wag 
|, | , - l . as 
ete Cee e 4 cele he Stupt-> rt “aworne 
*- 
a 44 
Sry eo. | iy 1 2D0k : 
Jz7103 0) 7S ' iS &) ; . * 
s sedi’ | 
5 
r i re. ; 
ov @&6t- ao, Ls isoprog (le } 
‘es } 
~ ob a - 
@ 
al 4 d r wu sy 
one arsitow atl i oe 
4 ” * 
} e 
> en ' whit + Lies is > st 
z \ a 
ia 26 ere vq Barcel foo eawutic 
= =) a> & P 


- he 


TT 136 béitosiw: zaw Uv spsbiiow ob: ait ‘pie .< 
aay YS) 


? 
Pare 
a ‘ 
ee 
a ts 


jicts wav leek) ats! noizoge ab wes rule ; 
4 7 é 


no dw tan? Jusmetete oct eye ithe tlle cf 


wl a ee ~ @e 


4 


- 





* 


, yaw. V apazton 7 rhans a en 
2 bor 


. ate parte up ai a (oi) samt sd 
ae 
pleow ra tied ante Jn 


tae Oras 


ie Sa. es 7 
Mawsd as oe si a ” 


' ae, 
<. 


a hs ae Fil 





60 


"00S Sut4RZes uoT{ANnTOSSy “SZTOA 7°Z ebeRTOA UuoTReAZETEDNe UOT 


*zesATeue sseu otTodniapenb syj z0Z ssew snszea szOROeJ uoTSSTusuerzy OLT?Z eanbTty 


002 2/001 wea 





| Papeete = j eT RE Nec cerTe Y 


RSE PY ra TT ROE 


PRES <2 a Re mer 


CR LUN tee 


SS sere 





RE ST 


uBiy WINIpOWw MO} 


ebubs sSsDW 


aL 





NVaL 


- 


~ 


2? 


a 


> 


S| 


ep 


<< 


«. 


(°%o) NO 





when run in the total collection mode was close to 50%. 
Probably the answer lies between, some ions were lost in 
transit from the cone [16] to the entrance aperture [20] 
and the quadrupole did not pass all ions in the total 


collection mode. 


2:11 Discrimination in Detection by: the Electron Multiplier 
The ion signals in the present apparatus were detected 
by a 14 stage, copper beryllium secondary electron multipler 
(SEM). Unfortunately an electron multiplier does not detect 
‘all ions with equal efficiency. The number of electrons 
ejected from the first dynode by the collision of the ion 
depends upon its mass and its chemical composition. The total 
amplification is proportional to the number of electrons from 
the first dynode. It has been reported (64) that the secondary 
electron current is proportional to the velocity of the ion. 
Now if ions of different masses are accelerated towards the 
first dynode by its negative potential, they all gain the same 
energies and thus their velocities will be proportional tc 1//Mass. 
The relative electron currents for two ions of masses my and ™., 


will therefore be related by equation (2:xxii). 
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Stanton,Chupka and Inghram (65) showed that the chemical 
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composition of the ion may affect the amplification. It 

appears that multi-atomic ions break apart on collision with the 
first dynode and act as a group of particles of the same energy. 
The amplification thus increases with the complexity of the ion. 
For example, Stanton et al found ieee had an amplification 
factor 2.2 times that of Eee 

The amplification factors for the various ions were 
measured in the present apparatus. First the amplified ion 
Signal was recorded when the SEM was operated in the normal 
fashion. The signal was again recorded when the first dynode 
of the SEM had been connected to the recorder so that it acted 
as eco cup. The ratio of the first signal to the second 
was the amplification factor. The amplification factors are 
shown plotted versus m/e in figure (2:11). One can see that 
the chemistry of the ion has a larger influence than does the 
mass upon the amplification factor. Equation (2:xxii) has 
been plotted so that it passes through Kr" m/e 84. The 
atomic ions fit this line more closely than do the multi atomic 
ions. 

An excellent discussion of the operation of secondary 
electron multipliers has been given by Beynon (64). The problem 
of amplification is discussed in detail. 

The measured intensities of the ions were corrected for 


discrimination by the quadrupole analyser and the electron 
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multiplier. Each ion intensity I was multiplied bya factor F. 
This factor was calculated for each ionic species by equation 
(229xXeEI12) . 


a 





(2:xxi ils) apie ee cere 
Quadrupole Sensitivity Factor x SEM Amplification 
For most of the results relative factors rather than absolute factors 
were used. Consequently the absolute intensities are not given in 


figures in which the ion intensity is the ordinate. 


2:12 Collision of Ions Outside the Ion Source 

The collision of ions along the trajectories from the ion 
source to the point of detection may cause discrimination. Two 
experiments were carried out to determine the extent of collision. 
The total ion current was measured as a function of the ion source 
pressure. It was measured on the cone deflection plate [l6a] and 
by means of the guadrupole analyser operated in the total ionization 
mode. The results are shown in figure (2:12a). One would expect 
the signal would increase as the pressure increased provided that 
at all pressures the number of ions leaving the source was 
proportional to the number formed by electron impact. After the 
initial rise the cone signal became reasonably constant, and 
then rose slightly as the pressure outside the ion source i.e. 
in the vacuum chamber reached ane = torr. The signal measured 


by the quadrupole, however, decreased slightly. The decrease 
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Total ion current as measured on the cone (electrode 16a) 
and by the quadrupole operated at a/q = O versus pressure 
of nitrogen in the ion source. 


a. Normal operating conditions. 


b. Operation in the presence of unwanted gas leakage. 
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is probably due to collisions in the quadrupole tube as the 
flight time of an ion through the quadrupole analyser is 
relatively long. For Ny with an ion source potential of 
+7.4 volts and a cone potential of -67 volts, the flight time 
of the ion from the ion source to the cone can be calculated 
as 4.5 micro seconds, while the flight time through the 
quadrupoles is 25 micro seconds. In figure 6:12b) are shown 
similar plots in which the leakage of gas from the ion source 
was much larger than normal (F = 25 cc/sec. compared to 9 cc/sec.). 
This was due to an unwanted leak. In this case the exterior 
pressure increased much more rapidly as the ion source pressure 
increased. One can see that collisions inside the quadrupole 
tube reduced the signal by a large factor than before. Also 
the pressure of esti torr in the vacuum chamber was reached 
at a lower ion source pressure. Since the cone current shown 
in both diagrams (a and b) of figure (2:13) increase rapidly 
at pressures above fo. * torr in the vacuum chamber, the 
increase is probably due to formation of ions in the region 
between the cone and the ion source BES ray electrons from 
the electron beam. Normally the ion source operated under 
conditions where the pressure of ios torr was not reached, 
and thus the only ions observed would be those produced in 
the ion source. 


: ; + . + 
In the second experiment the ions N, and 0. were produced 
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Ion signals produced by the auxiliary ion 
source as the vacuum chamber pressure P 
was increased. 
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Ratio n,/ 0. versus pressure. 
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in the low pressure ion source. The intensities of the two 
Signals were observed as the vacuum chamber pressure was 
increased. Air was admitted to the chamber in two ways. It 
was admitted via the high pressure ion source and also by 
allowing it to leak through one of the mounting flanges. 

The results are shown in figure (2:13a). One can see that the 
ion signals increased until the pressure reached aie and 
then decreased. This decrease must be due to collisions 
inside the quadrupole analyser. The ratio N, /0," is plotted 
versus pressure in figure (2:13b). The ratio remained 
reasonably constant as the pressure was increased to 1x10 
torr above which it decreased rapidly. Since No" has a 


larger ionization potential than those of the other major 


+ 
and H_O, N 


Sonstituterits “of vair pu -eJ 20 4 
4 2 zZ 


would be lost by 


charge exchange collisions at a greater rate than would O 


2 


This is the cause of the decrease in the N, /0,” ratio at the 
higher pressures. From this experiment it was concluded that 
discrimination of ions due to collision, were not important 
if the pressure in the vacuum chamber was kept below 1x10 7 
torr. Most of the experiments reported later, were done 
under conditions in which the vacuum chamber pressure was 


kept below inane Toor. 
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2:13 Purification of Gases 

Krypton, xenon, neon, oxygen and nitrogen were obtained from 
1 liter bulbs manufactured by the Air Reducation Company. These 
contained very small amounts of impurities (about 10 ppm). Nitrogen 
was also obtained from other sources and purified. It was found 
that the nitrogen experiments were more consistent when N, was 
continually flowed into the ion source. The Airco bulbs did not 
hold a sufficient amount for flow experiments. Nitrogen from a 


tank. (containing 2:ppt-0O, andi 7 ppt HO) was passed through a 


Ps 


molecular sieve (Type 5A), to remove HO, through sodium-potassium 


2 


alloy to remove O, and again through molecular sieve (Type 3A), 


2 


which was kept at -126°C by a methyl cyclohexane-liquid nitrogen 


slush bath. It was possible to reduce the O, and HO impurities 


2 
to less than 50 ppm by this method. 


Ethylene was obtained from a Phillips research grade 
cylinder of purity 99.94 mole percent. Methane, from a tank of 


Matheson "Ultra Pure" which was labelled to contain 50 ppm CO., 


and 1.7 ppm C.H. was 


240 ueppm 10 He 


13 ppm N SxS. ‘ppm -C 


o rhe 38 
passed through a trap of 3A molecular sieve which was cooled to 


-196°C by liquid nitrogen and stored in a 1% bulb. 


2:14 Experimental Procedure 
Appropriate gas mixtures were first made in the gas handling 


plant Ghiguree2: 6)c im Inothe early experiments (CH ,-Kr, chapter 4) 
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and (CjH,-Xe, chapter 5), the mixture flowed from the 4 liter 
bulb into the ion source. It was found that the pressure 
decreased very rapidly and since it took several minutes to 
measure the intensity curve of each ion more reproducible 
results were obtained by remaking the mixture before measuring 
the next ion. The technique was improved by making large 
quantities of the gas mixture and storing it in an external 
bulb. This bulb was attached to the gas handling plant by a 
Granville Phillips variable leak so that the pressure in the 
G.H.P. was kept constant. The 4 liter bulb was used as a 
ballast volume. The gas mixture was analysed before and after 
the high pressure measurements by means of the auxiliary low 
pressure filament. It was found that gas mixtures of gases 


with widely different molecular weights (such as CoH and Xe) 


4 
did not fractionate while passing through the variable leak 
as they might be expected to do. Since for these mixtures 
one gas was only a fraction (%v 1/107) of the major constituent, 
the latter presumably controlled the flow characteristics. 
When the pure gases 0. and Ny were used, they were flowed at 
constant pressure directly from the source through the variable 
leak into the gas handling plant. 

The ion intensity curves were obtained in the following 


manner. The gas was irradiated continuously by the electron 


beam and a spectrum was recorded. This spectrum was used to 
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establish which ions were present and whether there were 
impurities present. The electron beam was switched to the 
pulse mode and a second spectrum at constant ion collection 
was recorded. The ion gate was then switched to the pulse 
mode, the mass spectrometer was adjusted so that only one 
Mass was recorded. The ion intensity was measured by 
increasing the delay time in steps. -The time was measured 
by the oscilloscope. At each time, the chart recorder was 
moved forward and stopped, so that the curves appeared as 
step functions on the chart paper. This method of obtaining 
the intensity curves was used as the amplifier was slow and 
integrated the signal from many successive pulses. The ion 
intensity curves were replotted on graph paper (at times 
corrected for the flight time tes with the appropriate 
factors applied to correct for the intensities for efficiency 


Of collection. 
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if. 
3 PHYSICAL CONDITIONS IN THE ION SOURCE; THE LOSS OF 


IONS FROM ION-ELECTRON RECOMBINATION, FROM DIFFUSION TO 
THE WALLS AND FROM MASS FLOW OUT OF THE ION SOURCE. 


The ions formed by the electron beam are lost by three 
mechanisms, by recombination with electrons, by diffusion to the 
walls and by removal from the ion source by flow of the gas. 

Each of these mechanisms has a different effect upon the ion 
intensity. The Pesan nation of ions with electrons is a second 
order mechanism. Since the concentrations are equal, the loss 

of ionization will be second order and not pseudo-first order. 

The losses of ions by diffusion and mass flow are both first order. 

Sampling of ions from any mass spectrometer ion source 
May involve discrimination against ions of different masses. 

This problem is especially complicated in sampling from high 
pressure ion sources since the ions are removed from a region at 
pressures of several torr into a region in which the pressure 
should be less than tone torr so that mass analysis is not affected 
by ion-molecule collisions. The ideal requirements are that the 
sample be representative of the ionic constituents and that the 
reaction system be undisturbed by the process of sampling. 
Unfortunately these cannot be fulfilled completely as the above 
processes, recombination, diffusion and mass flow affect the ions 
in different ways. The effect of these processes on ion sampling 
and on the suitability of the system for kinetic studies will be 


considered below. 


“1 


- ' 
a " ¥) 
e "eS 
of 
A A Tin 
_ 
>| 


eer Sat 


~ 
; iy a¥) » We. 4 {) yh Woy) . 
= oe sana eel i 
i ie e = 


: ’ Bry fT 
~ Seaeeht aged ‘sell ete mn ad wreseelo ort “yd a ply 


? 


, NOSIS 


-* 
[3ib so wed 


iis 
q 

‘A ‘ 
- ; UA 


-Spaune MOL BHP. wx amoecayaa 
wiv tony REA OS {MOU ET RNCOSN MT INASP 4 
f oJ al gent Mc sat ie 


gie ddiw iti moos 






















ay i 
ber mn ps : 
ont “ig 

sas) | 


Oe ae | 


fy 


sot 
* 


a 7 iy 


* 
F + 
_— % oo 
nis 2 £ el SaOTIDSLO + Brier 26 Oi Je mat 
ig | | ay 7? 
+ goule aie ao 
; it ,is ts enpissxtoscan® St gee Mima ff 
SO4 aut 4 is | cay nee 
- : 4 *? _—" “ ha 1o . bieness ad fikw * oy - 
Pip Gx) EE a tet ci2 Om sa‘. ; { 2 : nie 
. it _ i 
S$ dyod esp woll vane Ons ropen Tae “a ent 20» 
, eae > ey ‘3 i FOG 2=te y , 22 i ; ) 
ae a ia ; 
| etimostoeqe Kean yas Text Maw Toor 
} ty _, ited (a. 7TSMOs i2e 12 ez 4 
i ‘= we J 7 
az ; 
2 b to ehoi 2an » so its 
28 fois, ne ie \e?* , 
gh mont ealiqned.ct Seasusignes spre 3 
- = tf. < 
haveme: ete snot pga 


+, gpiber, S me 


aks toidw at, ie 


er 
oo 


= aX ie ? +4 


heads atte seu atl Ssi#ye= 


pag 


eis Smelt ast © tags. rapes fae 


off 4 
o9t. wah 


pebbeynee. 20 eagootg 8 


‘6 tid vee aiesetyats sopantet 6 


aero t sii toe ti & wok. aaem bits apie * 


ms vain 10 


eahnane i RO > sieianosg 
a ee s 


ark ' Sean ap oz x03 


abs. ear 4.3m 


Pe 


sotid. bas ajanns 192009 o£nas oe 














i 
y 7 4 #) 
! Ae 


73 


3:1 Recombination of Ions and Electrons in the Beam 
i Steady State Concentration of Ions Under 
Continuous Illumination 
An estimate of the extent of ion-electron recombination may 
be obtained by an approximate calculation of the concentration of 
ions formed in the beam. Figure (3:1) shows an enlarged diagram of 
the ionizing region of the ion source. The ion current I, may be 


calculated from the electron current - by means of equation (1:11) 


Cre id) I 


where Q is the cross section, n the gas density and & the path 
length. The electron current I, was measured on the trap at 
about 1UA when there was no gas in the ion source. This 
. ‘4 | 

corresponds to 6x10 f electrons/sec. The cross section Q for 

-17 vl 
N. gas for 4000 volt electrons may be taken (62) as 3x10 rou tih F 

2 
16 : 
At a pressure of 1 torr at 300°K,n=3x10 molecules/sec. Since 
we want the concentration of ions, 2 may be taken as 1 cm. If 
these values are substituted into equation (1:1ii) I. is calculated 
as Wex10*> ion/sec. If the average thickness of the beam is 
about 1 mm, and its width is governed by the slit width 1.7 mn, 
2 2 

the cross section of the beam is about 2x10 cm. Thus the 


number of ions/cc ies produced each second would be 6x1loO . 


This is a very high density of ions and many would be lost by 
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recombination. Since the secondary electrons also cause 
ionization the actual density will be higher. The change in the 
ion concentration I due to recombination is given by (3:i) where a 


tor kL) ar 


_—— oo + 1 
dt oT fe] I. 


is the recombination coefficient which may be taken as agnon! cc/ion 
Ce. : : ; 

sec. “for N. ions (66) )""Since*the *formation*of an -ionpby electron 

impact produces a new electron,the electron concentration [e] may 


be taken to be equal to the ion concentration I. Under steady 


state conditions dI/dt = 0 and equation (3:i) reduced to (3:ii). 


ne Zee 
C3521) O Tepes WL : 


The concentration of ions under steady state conditions I, is 


gaven by o(3:ii31). 


Caeini) I a.) TES 


1 ON i tne ions/cc sec. and @= cesT cc/ion sec. then the 
s 


; a 
steady state concentration oe 4.5x10 : Tops /cc. 


44. Tifetime,of an.Ion.in the- Beam 
The average life of an ion in the beam is given by the ratio 
of the steady state concentration to the concentration of new ions 


produced each second. This is given by equation (3:iv). When the 


(3:iv) % See 
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above values for I and rT are substituted into this equation 
the average lifetime of an ion in the beam is calculated as 


ez 7x01? seconds. 


iii Time for Removal of an Ion from the Beam by Diffusion 
At high ion densities (greater than 10° ions/cc) ions 
are removed from the beam by ambipolar di ffusson. (677 <p y206).. 
The time required for the removal of an ion from the beam is 
related to the ambipolar diffusion coefficient dD. and the 
q@istance ad by equation’ (3:v) (67, p.200;° 68, p."493). 


2 
%— 

tha elena) oe 

| D 


(3:v) 


+8) 


d may be taken as 0.5 mm, half the width of the beam. The value 
of the ambipolar diffusion coefficient is somewhat unknown. We 
might take a value of 100 ent yeec! from a table given by McDaniel 


.o8,~p 15 16); is then calculated as 25 micro seconds. 


"aire 

This value has the same order of magnitude as the lifetime of 

the ion in the electron beam. Consequently one would expect 

that recombination would have an effect on the sampled ion intensity. 
An experiment was carried out so that the effect o£f*2or-— 

electron recombination might be observed. The electron beam was 

pulsed and the ions leaving the ion source were collected 


continuously. The pulse length was changed while the pulse 


repeat frequency remained constant. The increase in intensity 
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eae with the increase in pulse length was measured. The pulse 
repeat time was set at 2000 micro sec. for two runs and at 4000 
for a third, as the ion intensity decayed after the pulse to 
very small values within 1000 micro seconds. Consequently the 
major effect on intensity would be brought about by the pulse 
length (from 10U sec. to 1000u sec.) The ree are shown on 


figure (3:2). The observed intensity I , as a fraction of 


obs 
the intensity under constant illumination I,, is plotted versus 
the pulse width which is Wore essed as a fraction f of the pulse 
repeat time. This is called the illumination fraction where 
the value £ = 1.0 corresponds to constant illumination. The 
plots ee shown in log form (figure 3:2) as the solid lines 

a, b and c. If the ions in the beam were lost only by first 
order processes, the intensity would increase directly 
proportional to the illumination time. The dashed line d 
indicates this situation. It can be seen however that the 

ion signal increased at a rate, such that the log-log slope 

is less than one. This indicates that processes greater than 
first order are responsible for loss of ions. The part of the 
time fraction scale between about 0.3 and 1.0, where the slope 
increases,corresponds to the situation where the ion signal 
does not decay between pulses. This would cause an overlap 

and the change in intensity would not be due to the increase 


in pulse length. The data of 0.92 torr with 4000 micro second 
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Figure’ 3:2 Plot of cbserved total ion signal I obs normalized to ion 
signal under constant irradiation I, versus pulse width 
expressed as a fraction (f = At /t,) of the pulse repeat 


time. 
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repeat time line c of figure (3:2) was plotted in a slightly 
different manner in figure (3:3a). The normalized intensity 
was divided by the time fraction f. When the pulse widths were 
short this time-normalized intensity was much larger than when 
the pulse widths were long. This same data is shown in a log 
form figure (3:3b). 

A theoretical calculation was carried out. It was assumed 
that equation (3:i) held that ions were only produced by a linear 
relationship (Tt and lost by a second order process (-at7t) . 
Equation (3:1) may be integrated to give the following equation 


(3:vi) (67, ps 79) where I, is the intensity at time t and I> 


1 
: = 7 spat #8 
(S2Vi) I, = ze oo 1) 

3, Ae his 
(c co + 1) 


the steady state intensity given by equation (3:iii). 


(32211) Ee Nea 


Equation (3:vi) was used to generate the curve showing the change 
in intensity with time, after the beam was switched on. This is 
shown in figure (3:4) as the solid line. The beam was considered 
to decay by a second order process after the pulse ended. The 
equation used is (3:vii). The decay of the intensity I after 


a 


t 


T _ et ee 


(1 + ar, t) 


(aevi2) 


each pulse is shown as dashed lines. The areas under the curves, 
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Figure 3:3 a. Normalized ion signal I. 


/I,, divided by the 
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illumination fraction f versus f. 


b. Above data on log-log plot. 
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which show the increase in intensity during the pulse, and the 
decrease after the pulse, were measured geometrically. These 
areas would normally correspond to the ion signal received by 
the mass spectrometer. These areas were normalized to the area 
of constant illumination. The results were plotted versus pulse 
length (as a fraction of the pulse repeat time) and are shown 
by line e in figure (3:2). One can see that the general shape 
of the curve is similar to the observed values. One would not 
expect it to agree exactly as many approximations were made in 
the Per eniationd: The decay was assumed to be entirely second 
order and the geometric areas were only measured approximately. 
Figure (3:4) also shows that the steady state condition 
was not established until the electron beam had been on about 


200 micro seconds. This time is much longer than the pulse 


width (10 micro seconds) used in the later experiments. 


iv Charge Density in Volume Sampled by the Leak 
The density of charges in the sample leaving the leak 


may be calculated from the current measured at the ion acceleration 


cone. The density is given by equation (3:viii). n is the charge 
(3: viii) n = qo 
Fe 


density, I_ the cone current and Fo the conductance of the leak. 
Cc 


-~9 : , 
The cone current was measured as 1x10 amps which corresponds 
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to 6x07 ions/second. The conductance of the leak was about 
1.1 torr cc/second. The charge density is then calculated 
as Re exi0- ions/cc at 1 torr. This is only one order of 
magnitude less than the calculated steady state concentration 
of ions in the beam. It indicates that the ions are probably 
not sampled directly from the beam but from just outside. 
3:2 Processes Occurring in the Region From 
Which Ions Were Observed 

4 The Effect of First and Second Order Loss Mechanisms 

Each short pulse of electrons produces a group of ions 
which decay during the time between electron pulses. It would 
be interesting to determine the mode of this decay. A schematic 
diagram of the ion decay after each ionizing pulse isshown in 
Sucre (3:54), 

First consider that the decay occurs by a first order 
mechanism. On each pulse an initial ionization I, is produced 


and it decays according to (3:ix) where I, is the intensity at 


= 
(a:3x) Ai = I_e 

time t after the pulse and k the rate constant for decay. The 

pulse width At. is assumed to be infinitely narrow, i.e. At 0. 

Let the residual intensity remaining from the previous pulse = C. 


fie total intensity at each pulse is A. 


(3:x) Bae dee © 


we 



























Nite ae abe 3° ooma s9ubuis ast) ha 


*Heab BAW 


fhuidn lio DRS rash BL eaten aia es 


et a oe 


7 5 p 7 i ' i Fz 
* oe ean: i a af ait fad) 
: eke 
(it > age SOTA » arti 2 ys “ule bers luc ~oo 
= gpeepeiaan et nc nity 2a/% sansa ihrer Fo aed 
Wi , )4 bibars G % } ; a } 
> wrth’ as 
carasoe pa eed Sted Teed oss a i08 ab ; 
‘ x * © 7 9 a6, x 
_ * : ra * ‘ es * ; ve Se ti . ' | 
: ; a! ‘1 ¢ ue : 
| res as aps 
i? ae ; ; i “4 g ee si ad 
. « : : lay’. i, sl a 
; a ‘Se 
72RD 2 i tana t pores PAs 3° pes 
8 en ane te a . ee ply de Sin 
Smoke tim HED mh BODO LO Lg Brow ei ie |. 10 aig 
Cc 2h aft ee - ae 
I u pening ae rDole tteewted smh wth? ® - 
stoma ££ ) Bend | pf / | 
‘ ‘, Gt i dees 
jeunes. & Jieah abrty to ebom eid ociopateb ote Wad ‘3 
! Silage ih * ; . 
; ' a = be 
ay aeoree? oeidk giiminol dose 10248, ac ad hor Y 
Wi) ae 
ie “eh 
a _ 
| a bhaaon S82. 
“A “s nf one fy ertt be) 200d 2 
~hie aetit = ye e1oeo YB OS! aa : 
oop eet wales tent iaisiat cf eas A aes 
22M eel” . 
: f 
is ywpeiotosinl att Al) oe aenety (thse oa 
48 | ; : 
; md 


Pak, it ' cary i ¥ 2: ore . 

rs * | ‘ : Oo ay ar ; 

7 eS. eta ei = re ani nh cal 
eit wee ie sits ga noe oath no 
—_ far. wie Ped ase pio | 
08, $f, oy out (gheatnatat ok oF Deanpanth : 

| 

= citi lai ie ; ave T 


a 


Figure 3:5 
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Schematic diagram of electron pulse and ion intensity curves 


at constant pulse width and variable repeat time a 


Normalized observed ion signal I ps’ too! divided by illumination 


fraction f, plotted versus the pulse repeat time tL. 
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85 
Thus the observed intensity I! at any time t is given by (3:xi). 
(3:x1) am Te: Are 


mie value of C is given by (3:xii) where t. is the time between 


pulses. 
Se x22) C = Ae i 


If equations (3:ix) and (3:xii) are. substituted into (3:xi) the 
observed intensity is given by equation (3:xiii). 


(325444) es ee ee 


t emer ere 
kt 
CRETE, 
Under conditions of continuous collection, the ion signal is 
proportional to the area under the decay curve. This may be 


found by integrating equation (3:xiii) between the limits o 


and 1B equation (3:xiv). 





c. - 
(3: xiv) ; uy rea bi Pas 
I. ze fe) bb 
(een) k 
re) ° 
(3:xv) -- Area = bo 
k 


Thus when the decay is first order the area is independent of 
the pulse repetition rate. However the observed signal is 
dependent upon the initial intensity wey and the intensity of 
the electron beam which is proportional to 1/t,- If one 


corrects the observed intensity to constant illumination then 


the signal is again independent of a 
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This calculation may be repeated for the Situation in 
which the decay of ionization is second order. The general 
equation for second order decay is given Oye (os xVaL) ome The] 
parameters were defined earlier. The integrated equation which 


(3: xvi) ti = we 


PO ef) 
fe) 


is proportional to the area is given by (3:xvii). 





(Sev t4) )  Atea t= a 


a at 
= 


In this case the area is dependent upon Ee and increases as t 
3 
increases. 
he vores onar rrom 0.92. torr No in which the major ions 
+ + 
were N, and Ny was measured. The pulse width At. was kept 
constant at 10 micro seconds. The results are shown on figure 
(3:5b). The observed ion intensity was normalized to constant 
illumination by dividing by the fraction of time f = (At_/t,). 
The initial rise indicates that the decay of the pulses at 
times less than 200 micro seconds is second order, whereas at 
long reaction times the decay has decreased to first order. A 
: ; ‘ eRe, + 
plot of the normalized ion signal from Xe gas (Xe and Xe, ) 
at 3.90 torr, with 20 micro seconds pulse width,also shown in 


figure (3:5b), demonstrates that the behaviour was constant for 


different gases. 
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ii Analysis of the Ion Intensity Curves by 
Means of the Einstein Equation 
The total ion intensity curves obtained for N, gas were 
analysed by means of the Einstein equation (3:xviii) (68, jaa 492; 
OFF PewlIo)eerMis equations for tea heaudauite ofjascloud of 
" ane /4Dt 


(3:xviii) Nee eet 


27Dt 
particles by diffusion in one dimension through a gas. Ny is 
the number of particles located at the origin of a one dimensional 
coordinate system at time t=0. N is the number of ions at a 
point x some time t later, D is the diffusion coefficient. At 
any time t a plot of N as a function of x has the shape of a 
Gaussian error curve. 

For the purpose of this calculation the electron beam was 
considered to be an infinitely thin plane. The ions would 
diffuse from the plane effectively only in one dimension and 
thus the above htt may be used. The distance x remains 
constant in the present situation. The variation of N with 
time was determined, by means of equation (3:xviii). Various 
combinations of D and x were used in the calculations. The 
results were plotted in semi log form as the initial intensity No 
was not known and the plots could be moved up or down until the 
maxima had the same height. These were compared with the 
experimental ion intensity curves which were also plotted 


in semi log form. 
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Figure (3:6) shows some of the measured total ionization 
curves and figure (3:8) some theoretical curves. It can be 
seen that the results do not agree very well. .The ion curve 
Speatned som h.55 torr N, at 385°K is the only one which has 
the same general shape. Very large values of the diffusion 
coefficient D were required to obtain ion intensity plots 
even close to the experimental ones. The position of the 
maximum of the ion intensity curve deserves some comment. 

The maximum time ue is found by taking the derivative with 
respect to t, of Sey (3:xvi). The maximum time is given 


‘by equation (3:xix). Since the diffusion coefficient D is 


(S221) t = 


inversely proportional to pressure (67, p. 201). One would 
expect that eee pressure increased that the time would 
increase proportionally. The time of the maximum of the 
experimental ion intensity curves were found to increase 
slightly, but certainly not proportionally to the pressure. 
If the measured distance x, between the plane of the slit and 
the ion exit aperture, is substituted into equation (3:xvii), 
large values of D (500 Eo Pere to 1500 er eee) are predicted 
by the values of the time 't_- 

These calculations confirm that at short times, the loss 


of ions was by means other than diffusion, i.e. recombination, 
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CES EE SS BS EE 


q 2:05 torr 
b 3:40 torr 
C 2°65° torr 





d 0-96 torr 
100 
10 
eee be ee mee — 
20 30 40 90 
time (ysec.) 
Figure 3:6 Observed total ion intensity versus time curves for 


pure No at 380°K, x = 0.67 mn. 
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c 0-067 50 
d 0-167 1500 
e 0-167 500 
f 0-167 100 


a Rene 


oot) 100 
time (ysec.) 


Theoretical curves calculated from the Einstein 


Equation (3:xviii) 


on 









(Havtners Coopadx 
| Ue e0-0 9 
‘pot yeoo 4 
Won... Poo 

per Baro boo my 

‘108 Oto “9 rie 
ee 


92 


as the initial loss was very rapid. At long reaction times 
2 30 micro seconds (from figure 3:6) the decrease appeared to 
be due to diffusion,as the log plots (figure 3:6) show almost 


a linear decrease in intensity. 


iii Effect of Mass Flow 

If the ion is removed from the ion source entirely by 
mass flow, the time required for an ion to travel from its 
point of origin to the ion exit aperture may be calculated by 


means of an equation (3:xx) developed by Kebarle et al (42). 


(3: xx) t = an (rr >) 
aEy ; 


This is the equation for flow from a hemisphere above the 
leak. Te ee conductance Fo = 1.13 cc/sec. for nitrogen. 
yr is the radius of the hemisphere and ro the radius of the 
leak. If r is taken as 0.67 mm and ry as .060 mm, the time 
t is calculated as 560 micro seconds. When r = 1.67 mm 

t = 8.5 milli seconds. This would indicate that the effect 
of mass flow towards the removal of ions is not as effective 
as recombination or diffusion of the ions. Since a large 
amount of gas flowed out of the electron entrance slit, one 
might expect some ions to be removed also. The same equation 
May be used. r = 5 mm and ro is negligible in comparison. 


At a distance of 5 mm the time for outflow will not be 
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affected by the length of the slit, and the flow pattern will 
be essentially spherical and equation (3:xx) may be used. 
When Fo = 8.1 cc/sec. and r = 5 mm, t = 30 milli seconds 


which is also a long time with respect to the ion decay time. 


iv Average Lifetime of an Ion in the Ion Source 


The average life of an ion was calculated from the ion 


intensity curve. This curve was obtained for N, gas. atwW.55° torr 


at a temperature of 395°K. The average time t was calculated 
from the values of the intensity I. Measured at many times t,. 
The time axis was divided into time intervals At each 5 micro 
seconds long. At each time interval At, the height I, was 
measured. The average value of the function where I = f(t) 

is given < (3:xxi). 


(3:xxi) E> SinBeathere? 
ak ab 


POAtKOL, 
at 
Since I = £(t) the average time t times the average height I 


equals the area under the curve i.e. equation (3"xxii). 


(3:xxii) t:Tt = - At, a 


Thus the average time is given by (3:xxiii). 


= 2 
(3:xxiii) Cos cle At. I.) 
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Peo aetorr N, (figure 3:7) were tabulated and the value for t 


was found to be 160 micro seconds. 


sas Conclusion 

The above order of magnitude tay Poa and experiments 
show that the ion intensity was affected by three loss mechanisms 
recombination, denito lanadnceiedat and mass flow. 

Under continuous illumination it appeared that the 
majority of the ions were lost by recombination in the beam. 
Since this is a second order process, ions undergoing first 
order chemical reactions in the beam would exhibit a kinetic 
order which is not first order. The average lifetime of an 
ion in the beam was calculated to be 70 micro seconds, at 
l torr. The time for diffusion out of the beam was calculated 
at about 25 micro seconds. Since these are the same order,a 
considerable concentration of ions would diffuse from the 
beam. It was shown however, that the ion density in the 
region sampled by the leak was an order of magnitude less 
than the steady state concentration of ions in the beam. The 
observed ions must therefore have been sampled from a region 
just outside the beam and consequently one would expect that 
recombination of ions and electrons in the beam would not 
adversely affect the kinetic studies. At higher pressures 


the average lifetime of an ion in the beam would decrease 
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whereas the time for the ions to diffuse from the beam would 
increase. Consequently even though the ion concentration in 
the beam increased the number diffusing out may not increase. 
This might explain the observed ion intensity versus pressure 
plot (figure 2:12a). The ion intensity increased up to about 
1 torr and then remained constant. The loss NG oe from 
a predicted continual increase with pressure could not be 
explained completely by loss through ion-molecule collisions 
outside the ion source. Thus the ion intensity measured 
outside the ion source was probably representative of the ion 
concentration in the sampled volume. 

When the electron beam was pulsed, the processes occurring 
appeared ee be somewhat different. If the rate of production 
or ions is nies ions/cc sec. then eeehG,: ions/cc would be 
produced in 10 micro seconds (the usual time used in most of 
these experiments). This is less than the steady state 
concentration. The ion intensity plots (3:6) indicate that the 
loss of ions occurs initially at a rate greater than first order 
and then decays to a first order rate. This was also indicated 
by figure (3:5b). One might estimate from figure (3:6) that 
the change of rate took place at perhaps 20U seconds, or from 
figure (3:6), that it took place at 200 micro seconds. Since 
figure (3:6) shows the ions intensity curves directly the value 


of 20 micro seconds is probably the better one. The ions were 
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removed from the ion source by the first order processes, 
diffusion and mass flow. The diffusion coefficients obtained 
from the time of maxima of the ion intensity curves were 
extraordinarily high consequently the Einstein equation 
probably does not describe the situation properly. Since the 
maxima were independent of pressure one might conclude that 
ions were formed close to the leak and that these were 
removed in a very short time. Since the travel time of the 
gas was large the ions were probably removed mostly by 
diffusion. Unfortunately the diffusion coefficient is mass 
dependent. Equation (3:xxiv) PS 201) shows the relationship 
between the diffusion coefficient of an ion and the masses m 


of the ion and M of the gas. 


(3S Exxtv) Da/M+t+m 





Thus D will increase with the mass of the ion and decrease with 
increasing mass of the gas. Consequently there would be some 
mass discrimination and this will affect the kinetic measurements. 
The decay of the total ionization will affect the kinetic 
measurements only if the decay is second order. If the decay 
is first order the kinetic measurements will be correct. 
Consider that the total number of ions per pulse = ee 


Since this decays by a first order mechanism the observed 


R ‘ : , ‘ . t 
intensity at time t, with no reactions occurring is I, . The 
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presence of a reactant gas, which reacts with the ions by a 
first order mechanism, will cause the decrease in ions to be 
faster. If one takes the measured concentration of TI * at 
fe) 
F , \= 
some time,t i.e. Io , and adds a small concentration of reactant 
: ; : t ; : : 
gas the intensity will fall to I. - Since the mechanism is 


first order the intensities may be expressed by equation (3:xxv) 


(3:xxv) Qn sa = kot 


t 
ai 
c 


where k is the reaction rate constant and the concentration of 


reactant gas. To measure k take times t, and t.. 














1 ie 
9 i ie 
F ; a = Ph be a 
(3:xxvi) ke(t, t,) Ln I, Qn I, 
7 t2 aR 
c c 
c VE 
3 ee “Y = 2 
(3: xxvii) ke(t, t)) Qn I, I. 1 
cto ie Taek 
c fe) 


If the ion intensities are normalized by dividing out by the 

total ionization Io at all times the fractional ion intensities i 
; ; ; : Lyall 

will give the correct rate constant since i, = I. /1, « Thus the 

decay of the total ionization by a first order mechanism will 

not disturb the kinetic studies. This technique was used in 


most of the kinetic studies described in the next sections. 
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4. REACTIONS IN METHANE - KRYPTON 


4.1 Introduction - Purpose of Studying CH ,-Kr 
System, and Previous Results 

The technique of using a pulsed electron beam and gates 
ion beam to observe the rate of [es of intensity, of ions 
produced at very high pressures was new, and consequently some 
chemical verification of the method was required. Since the 
pulse technique was designed to measure absolute reaction rates, 
an ionic reaction with a well known rate constant was needed. 
All well known ion-molecule reaction rate constants are in the 
order of vale to toa cc/molecule second and were measured at 
low pressure. Probably the most frequently measured rate 
constants are those for reactions (4:1) and (4:2) in methane. 


These reactions are frequently used for calibration purposes (69). 


k 
a 1 + 
(4:1) Gi, LOCH gee aeecront ee) ROMS 
+ Ky + 
(4:2) CHa. °C, Seem beg + 3H. 


They were first observed by Talroze and Lyubimova (5) and by 
Stevenson and Schissler (6,7) who measured their cross sections. 
Similar measurements were made by many other workers in conventional 
mass spectrometers. These early results are contained in a review 
by Field et al (33). Unfortunately since the ions were reacting 
under the influence of a repeller field they did not have thermal 


energies. Rate constants for thermal ions, calculated from these 
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cross sections by the equation (l:vii) (see section 1:2), were 


found to be slightly dependent upon the repeller field strength, 


(1: vii) wae jee 
= & 2mp 


and therefore the use of these values for thermal conditions is 





somewhat ambiguous. 

Values of the rate constants ae thermal energies have since 
been determined at low pressures by means of a pulsing technique 
(16,18,70,71,72,73,). This technique was described earlier in 
section (1:3). 

The results for ky (reaction 4:1) and ky (reaction 4:2) 
are summarized in Table 4:1. One can see that the different 
workers have obtained very similar results, since those values 
for ky obtained by pulse methods agree within 10%. Fewer values 
of ky are available. When CH, is ionized by electron impact 
the intensity of ae is less than that of gee Thus Kk, is 
more difficult to measure and the uncertainty of its measurement 


is larger than that of k However the value of 0.79 for Kk, 


1° 


compared to 0.99 for k, from the repeller method (reference 33) 


ui 
shows that the two reactions have about the same rates. 


From the data in Table 4:1, one can calculate a suitable 


pressure at which the reactions may be conveniently observed. 
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TABLE 4:1 
Rate Constants for Reactions (4:1) and (4:2) 


Measured at Low Pressures ~ (Literature Values) 


Technique Rate constants (cc/molecule sec.)x10° Reference 
a . Ko 
Repeller 0.99 0.79 (33) other values 


are quoted in 
this reference 


Pulsed . ad te (16) 
TAP oe (70) 
phe 0.70 C71) 
eae (18) 
0.96 (72) 

0.86 (73) 


The half life of a bimolecular reaction such as (4:1) is given 


by equation (4:1) ( 74). 


(4:1) TO = ub 
k, [CH] 


where [CH] is the concentration (in molecules/cc,if ky is in 
cc/molecules sec.) In this instrument it was most convenient 
to set the electron pulse duration at 10 micro seconds. Since 
the pulse width must be short with respect to the reaction 


time a convenient half life of the reaction might be set at 
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50 micro seconds. Using equation (4:1) and a value of Ix1o77 


ror kis one obtained a CH, concentration of aeieee molecules/cc. 
At 300°K this is about 10° torr. 

Our intention was to calibrate the ion source at high 
pressures, of the order of 1 to 5 torr. To get such pressures, 
and yet still have the pressure of CH, about 10> torr a charge 
exchange dilution technique (which is also termed inert gas 
sensitization (42)) was used. The reactant gas is mixed with a 
large excess of the inert gas. A rare gas or one of the 
dintoms me see such as nitrogen is used as the inert gas 
sensitizer. If the ionization cross sections of the two gases 
are not too different, the concentration of ions produced from 
the reactant and the inert gas will be approximately proportional 
to their partial pressures when the mixture is irradiated. In 
our example these pressures are ‘is ae and 1 torr and would lead 
to the ion ratio of 1:1000. However if the recombination energy 
(R.E.) of inert sensitizer gas is larger than the ionization 
potential Ty of the reactant gas most of the ionization 
will be transferred to the latter through charge transfer 
TeertiOma (seCcionm bso, 2 sihd. 

Krypton was chosen as the inert gas. The Kr ion has 
two states. The Sucys with a recombination energy of 14.00 eV 
and the “Py, with a R.E. of 14.67 eV. These react with CH 


4 


+ + 
to produce CH, and CH, respectively (27). The overall 
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reaction scheme may then be written as follows: 


k 
4 + 
(4:3) Kr oy, + CH, aes CH, + Kr 
k 
+ 2 4 + 
(4:4) Kr (oP) lt om ———- CH +H + Kr 
% 4 | 
+ 1 + 
: — 
(4:1) CH, + CH, CH,. + CH, 
+ Ky + 
: —— 
(4:2) CH, + CH, CoH, + H, 
+ Ke 
(4:5) CH, + CH, —— No reaction or very slow 
k 
+ 6 : 
(4:6) CoH, + CHy aan. No reaction or very slow 


The unreactivity of Tian and Cue ions in methane was 
previously established by Field et al (34), by Kebarle (75) and 
by Wexler (39). This then sets a limit to the reaction sequence. 

The presence of both CHE and ere is unfortunate since 
it complicates the kinetics. Two other gases with recombination 
energies larger than the ionization potential of methane, and 
which might have been conveniently used, are xenon and argon. 

The choice of sensitizer gas was made by referring to several 
studied of charge transfer reactions of the rate gases with 

methane (27,38,/76,/77). The results of von Koch (27) were the 

most useful. Figure (4:1A) shows a breakdown curve for CH, 

taken from this paper. This curve indicates relative concentrations 
of the ions that are produced by charge exchange, as the 
recombination energy of the primary ion is increased. 


2 


The use of xenon (R.E. Py = 13.44 eV) would appear 
2 


; : + 
ideal, since by exchange only CH, is produced. Unfortunately 
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Figure 4:1 From von Koch reference (27). 
(a) Breakdown curve of methane; the mass spectrum 


as a function of the energy absorbed. 
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(b) Distribution of energies given CH 


4 by electron 


impact. 
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however Field and Franklin (78) have reported that in mixtures 
: + + + + 
of xenon and methane the ions xeCH, F XeCH, , XeCH , XeC and 
XeH” are observed. These ions would complicate the kinetics 
exceedingly. Argon might have been chosen since it produces 
+ : 
a large intensity of CH, , but the rate constant for reaction 
a. 
(4:2) is less well known and CH, is also produced. 


. 


As was mentioned earlier, the break down curve indicates 


+ 2 ; 
that of the two states of Kr , the Sy. reacts with CH, to 
2 
+ + 
give CH, and the te) to give CH, . Cermak and Herman (77) 
2 


reported the relative yields to be 0.65 and 0.35 respectively. 
Krypton ions react with krypton atoms to form a molecular 


ion, by means of a three body reaction (4:7). 


+ + 
(4:7) RY. Tau2ikr 7 Kr, + Kr 


This reaction is not expected to interfere, since at a relative 
+ 

pressure of methane to krypton of 1:1000, Kr ions will be 

removed much more quickly by the charge exchange reactions 


(4:3) and (4:4). 


4.2 Methane-Krypton Experimental 


The variation with time, of the intensities of the ions 


+ 


+ + 
a ee and Con was measured for mixtures of 


Kr’, CH.", CH 
: 4 3 . oa 


3 4 


methane and krypton. The methane concentration was varied over 


3 


~ =3 
a range of 0.68x10 to 3.4x10 torr in a constant pressure of 


3.4 torr krypton. All experiments were done at room temperature 


~~ 
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(i.e. 297 to 300°K). The electron pulse and ion gate widths 

were set at 10 micro seconds. The pulse repeat time was set 

at 2000 micro seconds so that the ionization produced by each 
electron pulse had decayed completely between pulses. All 


other conditions were the same as described previously (section 2). 


4.3 Results 

i Qualitative Observations 

The variation of intensity of the major ions onsered in 
a mixture of eases Lee ara CH, ieee COT nes Shown in 
figure (4:2). This represents a set of typical results. The 
plot Eino shows the change of total ion intensity with time. 
As can be seen, the decrease of total ion intensity was very 
rapid for the first 50 micro seconds. The decrease then slowed 
down so that the ion signal became negligibly small after 2000 
micro seconds. This decrease in total ionization in the ion 
source was due to ion-electron recombination, ion diffusion to 
the wall, and mass outflow through the leak. These factors and 
their significance to the kinetic measurements were described 
eerliem in section 3:3). 

One can see qualitatively in this figure how the reaction 
sequence is carried. The Kr" ions initially formed travelled 
towards the ion source exit and there was an initial increase 


s ata . * . . . . 
in the Kr intensity, and hence the cOtal ion intensity, since 
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Figure 4:2 Raw data showing decay of kr , CH, » CH and total ionization. 
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i + + 
at short times (0-10 micro sec.) Kr was the major ion. Kr was also 
being produced in this time by the 10 micro sec. electron pulse. As 
+ : ; 
the Kr ions reacted they disappeared more rapidly than the 
oh bel ; + ao : 
total ionization and CH and CH appeared. These in turn 


a a 


+. + 
decreased as CH, and CoH, were formed. 


ii Preliminary Estimate of Rate Constants 
In order to gain a preliminary set of values for the rate 
constant one may try to select conditions where the decrease of 
an ion intensity due to reaction is faster than the decrease due 
to diffusion and mass outflow. It was assumed that this was the 
. ae oF oe 5 ; 
case for the reactions of Kr , CH, and CH, - One may derive an 


4+ 
equation for the loss of Kr as follows. For the general 


bimolecular reaction (4:8) 


(4:8) e + CH, —— products 
+ + 

ay LUGHAD ee S'S) (Gea bm cece 

(4: ii) at = 4 


Since the concentration of ions is much less than the neutral 
concentrations i.e. [CH] >> ee Ghe we may consider that [CH] 
is constant, so that the reaction becomes pseudo-first order. 
We can rearrange (4:ii) and integrate (4:iii), Let [m7] be 


+ 
written as Ff . 


(4:iii) aileeeek cH pdt 


I mca" (CH 1 € 
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This can be written alternatively in an exponential form. 


+ = t 
(aby) oe if Spe eane 


From equation (4:iv), if log, {T'] is plotted versus t the 
slope becomes ~k[CH,]/2.303. The equation ees ie practice then 
becomes (4:vi), if equation (4:iii) is integrated between times 
Aa and t,- 

(4:vi) k \= -2.303 log, (I¢,/It,) 


[CH] (t, a t,) 


: + 
One may also use this equation for the loss of CH and 


; 4 
a, Cars . ‘ : 
CH, ‘if an additional assumption is made; that reactions (4:3) 
and (4:4) had gone essentially to completion before reactions 
(4:1) and (4:2) occurred so that these ions were not being 
produced at a time when their decrease due to reaction gave a 
linear log plot. 
The log plots of the intensities (corrected for mass 
Hh ae eee : oP + oF 
spectrometer sensitivities) of the ions Kr , CH, and CH, 
for a mixture of tddxion: torr CH, and 3.4 torr Kr are shown 
in figure (4:3). The respective rate constants are calculated 


: and pheasion- cc/molecule seconds. 


as 1.97x10 >, 1.50x10- 
+ 

The rate constant measured from the decrease of Kr is a 

composite of those for reactions (4:3) and (4:4) as the two 


+ F ; : . 
states of Kr were indistinguishable experimentally. 


The preliminary values for the rate constants for three 


am few. 
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methane concentrations are given in Table 4.2. 


TABLE 432 


Rate constants from the decay of 


Kr 
Pere = 5. torr 


Pp CH, x10" jelowaie 


The values have the 
by other techniques 
factor of two, The 


the reproducibility 


+ + + 
CH CH, 


4 - preliminary va 


lues 


“Rate constants (cc/molecule eaeehelor 


k (Kr’) Tec oc 
Rx (4:3), (4:4) Rx (4:1) Rx (4:2) 
2.21 1.60 1.94 
2.15 ess 1.85 
1.97 1256 iced 
1999 1937 1.48 


same order of magnitudes as those measured 


(Table 4:1) but are different by almost a 


EWOm DUM Smelt 0.68x10°> ‘opene 


was good, as the difference 


runs -is less than 103%. 


The total ionization is also plotted on 


taken from figure (4:1). 


One can see that the 


time is almost as large as the decrease of the 


intensities, 


CH, show that 


4 


in duplicate 
figure (4:3), 
decrease with 


individual ion 


and since this will make each ion intensity 


decrease appear to be faster than that from reaction alone, 
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these rate constants will be too large. 

One can see also that these rate constant rates are not 
independent of methane pressure. There were two possible 
explanations; firstly that the reactions were not really 
pseudo first order, or secondly that the decrease of ion 
intensity due to mass flow was not the same at each methane 
concentration. Careful comparison of total ionization curves 
at different CH, pressures (not shown) showed this was the 


case. Consequently: a normalization procedure was applied as 


follows. 


iii Corrected Values of Rate Constants 

The data shown in figure (4:3) were normalized by plotting 
each ion intensity at a time t as a fraction of the total ion 
intensity at that time. The normalized curves are shown in 
figure (4:4). The reaction sequence can be seen quite clearly 
in this figure. 

Plots similar to (4:2) and (4:4) were prepared for other 
data resulting from the irradiation of mixtures 0.68x10 > and 
3.4x10 > torr methane in 3.4 torr Kr (these are not shown) 
subsequently the fractioned intensities will be labelled with 
asmall i (e.g. i By, 

Kr 


The data in figure (4:4) represents a system of consecutive 


(pseudo) first order reactions, since the concentration of CH) 
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Figure 4:4 Ion intensity curves from (4:3) normalized total ionization. 
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is constant. Reactions (4:3), (4:1) and (4:5) and reactions 
(4:4), (4:2) and (4:6) can be treated as two parallel reaction 
chains. Mathematical formulae are readily available (79). 
First consider one of the chains. This is a single two stage 


‘ F ‘ : + 
reaction with a single starting substance kr . 





+ 2 cA + 
(4:3) Kr ny, 2 es aoa petal 4 
+ bak + 
(4:1) CH, + CH RS CH, + CH, 
. ’ ‘ 
(4:5) CH, + CH, > no reaction 


The concentrations of the ions are given at any time t by 


“equations (4:vii) to (4:1ix). 


re) ~k,[CH,] a 


(4: vii) tk 8/, = dey 3/, 5 


(4: viii) L aoteR Levba / US, ook ; 
4 2 Lk, -k,) 
k, shee 88 | 
(k,-k,) 
(4:ix) 1 gat i ateg eee et ent) 
5 y) 2 4 


A similar set of reactions may be written for the chain starting 


: +2 
with Kr Pos 
78 


+ 2 2 
The initial abundances of Kr in the Ps/ and the ae may 
2 


uaa ; 
be expressed as fractions fy and f of the total Kr intensity 


: 3h: 
i° +. At t = 0 the intensity 1 1 and thus 


+= j = 
Kr Kr ttotal 


di Lo 
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fe) f fe) 2 
i cane dete. ae ? i i 
it P35) Xt and i... ( Py) Xt At any time t, i,t is 


given by (4:x). 


(4:x) ae OL CP 


he Pals. +(°P, ) 
Kr K Kr ls 


2 


Since one cannot differentiate between the two states of krypton 


_ + 
with the mass spectrometer, one can obtain ler by 


substituting the two form of (4<Va3) intow(4:x) 


, ’ Hy =e CH bt -k [CH Jt 
(4: xi) gee = fy Cnes 4 f, e 4 4 


The values for f, and f, will be given by the relative abundances 


1 


+ + 
at long times of the ions CH. and CoH. (see figure (4:4). The 


observed relative abundances, when averaged over several runs, 
were found to be in the ratio 0.75 to 0.25. These relative 


abundances for ag and cH,” produced initially by charge exchange 


+: : 
from Kr are similar to expected values. The theoretical ratio 
depends upon the transition probabilities. These are given by 


(2S+1) where S is the spin CZ ox jee ihe sat.o then is 


si ioe +P, =94:2 = 0.67 : 0.33 similar to that obtained 
Kr 3/, Es % 


by Cermak and Herman (77) of 0.65 and 0.35. 


Le k, = ky equation (4:xi) can be reduced to (4:v). Plots 


Of tog i aes versus t for the three methane concentrations are 


K 

shown in figure (4:5). From the straight line portions of the 
+ 

plots values for the rate constant for the disappearance of Kr 


are obtained. Now from equation (4:xi) log it would only 


‘ : : ae 
give a straight line over the whole range of Kr if k, = ky: 


7 ae 3 ay 


ss ~ tl at -* , + 
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Figure 4:5 Decay of i, versus Le 





However the linear portion in figure (4:5) does not extend over 
the whole Kr’ range, but from perhaps 50% to 1%. This is 
consistent with either kK, = kK, or k, < kK, since if reaction 
(4:4) were faster it would be complete in a shorter time than 
meaction (2: 3)\i.¢.”inathe range of Kr’ from 100% to 50%. In 
figure (4:5) the initial slope on the 0.7 and regis ta te tore CHy 


runs is approximately twice that of the linear portions. The 


+ 
maximum for the time dependent curve of CH 


3 always appeared at 


-- 
a delay time shorter than that of CH (see for example figure 


(4:3) or figure (4:4). Both of these observations are consistent 


Wiin-k, <‘k Thus the value of the rate constant calculated 


3 4° 
from the linear portions of the plots of figure (4:5) should 
be wthat of reaction (4:3) (i.e. Kk, = Pons cc/molecule sec.) 

4 
but a reliable value could not be obtained. Runs at lower 


concentrations (than 0.7 millitorr) of CH, where the portions 


100% to 60% of it should have lasted longer gave unreliable 


; ip Ge , 
results since under these conditions the decrease of Kr intensity 


due to reaction became much too slow compared with the decrease 
due to other processes (diffusion, etc.). Thus there were 
large errors in the fractional intensities. 

Equation (4:viii) for the concentration of the second 
ion in the consecutive reaction sequence is used for the CH 


4 


; ' , sae ; : me 
ion. The intensity of CH, is given by equation (4:xii). 
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k, appears to be approximately twice k(%2.6x10 cc/molecule sec.) 
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eb 28 6b geet ag a ee [evs tcH lt : ork cH It] 
CH, af ee 
Ths 3 
which is derived from (4:viii) where add = f xi. The factor 
f, = 0.75 is that used in equation (4:xi) as the initial concentration 


ui 


+2 oe ‘ ead ape 
Cres oh. 4.) nh Set Per vequdcion (4:x1311) for CH is used where 





3/, 3 
£, | 0.25. 
Ciao, Ee me, CALS Jerk glcult qk, cH) ¢| 
CH 2 
3 K-k, 


The rate constant for reaction (4:1) can be obtained from 


+ 
through equation (4:xii). The simplest 


the time dependence of i 
CH, 


evaluation of kj arises when kK) is less than k.. Then the term 


eotegl® is very small; that is reaction (4:3) may be considered 


to be essentially complete. In this case a plot of log (Log +) 


4 
versus t should give a linear portion with slope = sill 
PE OS, 
Such logarithmic plots of doy + for the three pressures of 


4 
methane are shown in figure (4:6). 


om , : 
Using a similar argument for CH, , that reaction (4:4) is 
“k, (CH, I]t oe 


essentially complete before (4:2) is observed, the terme 


: eras : v 
equation (4:xiii) becomes very small anda plot of logi versus t 


RCH 4] 
will have a slope approximately - 2 4° . These plots are 
22000 


shown in figure (4:7). 
The rate constants obtained from plots such as those in 
figures (4:5), (4:6) and (4:7) are summarized in Table 4.3. 


Duplicate runs of each methane pressure are shown. 
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Decay of the Ions kr cH and cate 
Pressures Rate constants (cc/molecule Eats: 
CH, (mi lili torr) Kp torr k Cet) ky Gog +) K Gog +) 
4 3 

0.68 o74 a 1.04 pee oe o4 Go 102-89) 

Oa7 3.4 OF Ovc 0.8 

1.4 3.4 TE AR ois Dae O78 7 Oy Ose 

i! 3.4 1.04 0.96 0.95 

3.4 ee Seek be 50.99 c (Oc 29 

3.4 3.4 123 0.95 Les 

Average values ral on POM O59 O59 

a from figure (4:5) 

b from figure (4:6) 

c from figure (4:7) 

d Average values transferred to Table 4.4. 

iv “Preparation of Theoretical Plots 

As a final check the equations (4:xi) to (4:xv) were used 
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to generate theoretical ion concentration curves FOr 7 ee 
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ae + ; : 
CH and C,H. . Fyrom equation (4:ix) one may obtain the two 


5 Pie’ 
new equations (4:xiv) and (4:xv) required for i_ + andi - 
CH Beg | 
5 oa0 
Lat ; i Bs ok (Ch it 
(4: xiv) don + = f) le, e 3 4 ae +| 
5 4 
(4: xv) i Oe Clg ie Be 
CH, = f, Ee, eek 6 aCe 


Curves with various combinations of values of Ky rkorkyrkye 


fi and f. were prepared. These were then compared with 
figure (4:4). Figure (4:8) shows the curve which fitted the 


best. The values used to obtain figure (4:8) are listed in 


column b-of Table 4:4. 


TABLE 4:4 
-Thermal Rate Constants for Reactions 


+ = + 
or KE CH, and CH, with CH... 
A 38 


Rate’ constant 


Reacti 
ion (cc/mol.sec.)x10” 
ee Ue ee ig 
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xx (2p ) + CH y 4. CH y Per i Zao 
8 4 3 
CH + + CH ee CH t CH. 1.0 1.1 
; ; 7 ‘ 0.9 ‘ landed: 
k 
CH + cH, —> CH +H ee 7 5S eae 
3 4 2°5 2 ; ; ‘i aoe 





a Values from averages in Table 4:3. 
b Values used for theoretical plot figure (4:8). 
c Averages of literature values Table 4:1. 
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Bec Theoretical ion intensity curves. Compared to figure (4:4). 





4.4 Observations and Discussion 

The values of the rate constants obtained from the 
theoretical plot are probably the most accurate since they 
correspond to the most complete use of the experimental data. 
The values for Kk and kK, oF om) © cc/molecule second of 
column a of Table 4:4 will be low since they were calculated 


without accounting for the rates of formation of the reactant 


1 + 
ions CH and CH. . 


4 ce) 
The present results for ky and ky (column b) agree very 
well with the previous results (column c). This agreement is 


not necessarily to be expected even if both sets of data are 
correct since the ions reacting in the two systems have different 


; P : é ea : 
excitation energies. Consider CH, in reaction (4:1) for example. 


When CH, is ionized by 50 to 100 volt electrons, as in the low 


suas : ‘ F ‘ ‘ 
pressure studies, the CH ion is produced with a distribution 


+ 


of excitation energies. Von Koch (27) has estimated this 
distribution function which is shown in figure (4:1 b). The 
distribution function was determined by comparing the relative 
intensities of the ions produced by electron impact (50 to 100 
volt electrons) to the "breakdown curve" (figure 4:1 a) which 


shows the mass spectrum of CH, as a function of the energy 


4 


absorbed by the molecule. Since the energy required to produce 


Ss is 1.2 eV (14.3 - 13.12) above that required to produce 


+ Tg s 
CH, , the distribution function shows that the CH, ions 
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involved in the reaction studied (by electron impact) at low 
pressure, would have, on the average, an excitation energy of 
perhaps 0.6 eV. 


aah 
eh exnCH ion, produced by charge exchange from Ke P 


4 3/, 


would have an excitation energy of 0.88 eV (14.00 - 13.12). 


However under the high pressure conditions used in the present 


work the cH,” ion will undergo about 1000 collisions with Kr 


atoms before meeting a CH, molecule. It is not certain that 
1000 collisions will remove all the excitation energy, as the 
Peeency of energy transfer in a Che =EKEeColaSLOn. 1Sanot 
known. One would expect that since the collision complex would 
be held together by polarization, it would be long-lived and 
allow energy exchange to occur. The rate of exchange would be 
much larger than that for collisions of neutral species for 
which the collisions are short lived. Very small concentrations 
of an ion of mass corresponding to (KrCH,)" were observed. Thus 
it is likely that in the present experiments the cH ions were 
thermalized. It is remarkable therefore that the values for Kk) 
agree exactly. If the assumption, that the CH,” ion is thermalized 


at high pressure, is correct, it indicates that the rate of 


reaction (4:1) to produce CH is not affected by the presence 


5 
of some internal excess energy. This situation was implied by 


the Gioumonsis-Stevenson equation (equation 1:vi) which defines 


the rate constant in terms of the polarization forces and not 
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on the internal energies of the ion and molecule. The complex 
between cH and CH, was shown earlier (section 1:5) to be held 
together by polarization forces. 

The values in Table 4:4 indicate that the measured rate 
constant for the reaction (4:2) of CH, was larger from the present 
high pressure determination than from the previous low pressure 
determinations. This appears when the rate constants k, and k 


1 2 


are compared. At high pressure the data indicates kK) = Kos 


whereas at low pressure ky > ko: If one assumed that this 
difference is real, and is not due to the limits of experimental 
accuracy, one might explain it as follows. One might assume that 
the same conditions existed for cH,” as for CHa that is, at 

low pressure the reacting CH, had some excitation energy, but 

at high pressure was thermalized. It was previously remarked 
(section 1:5) that the complex (CH,-CH,)” was different from 
(CH,-CH,)”, that it was bound by covalent bonds as well as 
polarization forces. The presence of excess energy might inhibit 
the formation of a covalently bound complex and thus decrease the 
rate constant. 

It is interesting that the rate constant for the charge 
exchange reaction (4:4) in which Ke Fy, reacts to give CH,” is 
almost twice as large as the reaction (4:3) in which Kr “P3,, 
re It may be possible to explain this in terms 
of the Massey near-adiabatic hypothesis (80). This hypothesis 
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states that for charge exchange from an ion to a molecule, 
having a low relative velocity, the rate will be small unless the 
change in excess energy AE is small. AE is the energy that 
must be converted into translational energy during the collision. 
For atoms and atomic ions, AE would be due to the difference in 
electronic levels. However when molecules are involved, the 
vibrational and rotational levels would be used, so that AE 
would depend upon them as well as upon the electronic levels. 
Conseguently, some knowledge of the cH,” electronic and 
corresponding vibrational levels is needed,so that the Massey 
adiabatic hypothesis may be applied. 

| Baker et al (81) measured the photoelectron spectrum 
of methane and observed it contained two maxima at approximately 
13.7 and 14.7 eV (compared to 12.70 eV, the adiabatic ionization 
a) 


potential of CH . They attributed the results to "transitions 


4 


from the methane molecular ground state to ionic states of 


symmetries lower than T 1.2.60 and Coy! formed by the 


a 3Vv 
removal of an electron from the highly bonding owe orbital." 


The C3y state is that in which there are three equivalent C-H 


bonds, and one non-equivalent C-H bond lying on the rotation 


wh 

axis. Baker et al also state that the CH, appears when the 
; 2 

critical amount of vibration energy is supplied to the C3, A 


ground state so that the non-equivalent C-H bond breaks. 


; a5 
Since charge exchange from Kr to methane produced two 
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different ions we might consider that two distinct processes 


ade y 
ion in the Ta 


+ eng 
occur, #that Kr aby produces initially a CH, 
2 


F +2 nit > ara “hg : 
state, while Kr F produces initially a CHy ion in the ae 
2 


state, which since it is also vibrationally excited breaks 
down to cH,” and H. The rates of these two reactions would 
depend upon the probability of transition from the methane 
molecule ground state (T, symmetry) to the two ionic states 
t and Coy! both vibrationally excited. One must therefore expand 
the application of the Massey hypothesis by considering also 
Franck-Condon factors. 

In order to explain the results, one might propose a simple 
potential energy schematic diagram for methane, figure (4:9). 
. The methane is considered, for simplicity, to be a diatomic 
molecule HC-H. One would expect that since the CHT and the 


cH, "T have the same symmetry that the bond length H,C-H would 


be essentially the same. The potential energy functions are 
shown vertically above one another. The transition of highest 
probability would then be between the lowest levels. Since, 
the ey leak ion has a non-equivalent C-H bond, and since it 
more easily absorbs vibrational energy one might expect it 

to be larger than the C-H bonds of the T3 state. This one is 
shown in figure (4:9). The probability would be large for 
transition from the low vibrational level of CHT; to quite 


+- 
high vibrational levels of CHy Cay" 
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Figure 4:9 


Approximate schematic of potential energy functions of 
methane if the molecule is considered to be a diatomic 


HCH 
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When Kr 25 charge exchanges to give CHy Ta: 


2 


difference between the recombination energies (14.00 eV) 

and the energy of the lowest state of SE (is, 12°ev). is 

0.88 eV. However since this may not be taken up as vibrational 
energy (according to the present hypothesis) the excess energy 


AE which must be transformed into translational energy is large. 


the difference 


+2 + 
For the case of Kr P, exchanging to give CH, C. , 
% : 4 3V 


between the R.E. and the lowest level of Gh,” cee is even larger 
(1.0 eV). However since the transition would more likely be to a 
high vibrational level, most of this energy would become 
vibrational energy, and AE would be small. Consequently, by 

the ester hypothesis, this latter reaction would have the 

larger rate constant. This agrees with the experimental 
observation. Figure (4:9) could also be used to explain why 


+2 ; : , 
Kr , with a recombination energy of 14.0 eV, would not 


1s 
2 


produce CH in the lowest vibrational level (about 3.7 eV). 


47a 
Since the transition would not be vertical, it would not be 
allowed by the Franck-Condon restriction. 

One might make one final comment concerning the reliability 
of the present pulse technique. We have shown that the results 
were consistent, and agreed an expectation from the results 
of other workers, and thus the pulse technique may be considered 


sufficiently reliable to be applied to kinetic studies of 


unknown systems. 
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5 REACTIONS IN THE XENON-SENSITIZED IONIZATION OF ETHYLENE 








5:1 Discussion of Previous Investigations and the 
Purpose of the Present Study 
Ion-molecule reactions in ethylene have been studied 
by many workers (39,41,42,82,83) who wanted to obtain information 
to explain the processes occurring in ethylene radiolysis. Of 
great interest was the ionic polymerization of ethylene to 


form large ions of the type CH (83,84). Electron impact 


2n 
of ethylene produces many primary fragment ions (32), a large 
number of reactions occur to form secondary and higher order 
ions. The reactions are simplified when xenon-sensitized 
ionization is used, as the number of different primary ions is 
reduced (42,85). 
At pressures of several torr of xenon the majority of the 
. hr + ar ; 
ions are Xe and Xe, A Xe, is formed by the third order attachment 


reaction (5:1). Both of these ions can ionize ethylene. Since 


+ + 
(S23) Xe + 2Xe > Xe, + Xe 
+ 2 : f F ; 
Xe 3 has a recombination energy (12.16 eV) which is larger 
2 
OAs, : : . - 
than the ionization potential for the formation of C_H (D0 SL eV, 


268 


the latter are produced by the charge exchange reaction (5:2). 


ae + 
, + ens: 


; + 
Charge tranfer studied (17,77,86) have shown that Xe aie 
2 


with a recombination energy of 13.44 eV, reacts with ethylene 
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ae 
to-wtorm,.C Ho ‘J (Ip. 13..93.,eV)4and CoH (Ip 13.5 eV) by reactions 


Fake) 2 


(5:3) and (5:4). Both reactions are endothermic. These studies 


+ 2 + 
Cy. 3) Xe Beer Cou a Ty + H + Xe 


2 4 Re 


wv 


+ 2 + 
(5:4) Xe Pie + Cl 2 at Cs + H + Xe 


2 4 ey 


wv 


were done with ions with energies slightly above thermal. 

+ 
Karachevtsev et al (87) observed that charge transfer from Xe 
ions at close to thermal energies produced the ions in the 


GT), ‘crs ("0-05 and cu. Cue 


following concentrations CoH, oH of, 


At slightly higher energies Cermak and Herman (77) quoted the 
relative abundances as 0.68; 0.09 and 0.23. Szabo (86) found 


Similar values 0.60; 0.09 and 0.31. The theoretical prediction 


+ 2 


: wai + 2 
for the intensities of Xe P34 and Xe Py produced by electron 
2 


: , : +; 
impact are 0.67 and 0.33. The relative intensity of CoE, to 


: ae at: 
the sum of the intensities of CoH, and eine observed by 


these workers thus supports the reaction mechanisms given 


above. 


+- 
The molecular ion Xe, (R.E. 11.16 eV) may also produce 


+ 
CoH, by reaction (5:5). 


+ + 
: > 
(5:5) Xe, + CoH, CH, + 2Xe 


Kebarle et al (42,85) observed that, under thermal 


conditions at high xenon pressures (up to 40 torr), charge 


- 
transfer to ethylene produced essentially only CH, ects 


would indicate that the exothermic reactions (5:2) and (5:5) 
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occurred and that the endothermic reactions (5:3) and (5:4) 
did not occur. Thus at high xenon pressures reaction 


7 A ; 
sequences originating from CoH may be investigated, 


4 
The following reaction sequence, for reactions of ethylene 


in xenon, based upon a similar sequence proposed by Field (32) 


for reactions in pure ethylene, was proposed by Kebarle et al (42). 


k * 
+ ~ 
(5:6) Cla wince, —b_. (C,H,.) (eddition) 
k 
se Ss ae 
(5:7) (CH yes Se SS C,H + Xe (stabilization) 
Sek. CH. + CH, 
(5:8) (C,H, ) yg (decomposition) 
CHL + H 
+ x 
(5:9) C,H. SNOUT 2 LS (charge removal) 
+ Koa on 
(5:10) C,H, + CoH, —— (CoH ) (addition) 
* Boa mR ae 
° tS ; + 
(5S 31T) CoH) > + CoH, (CoH ) (addition) 
etc. 


Kebarle et al (42) observed that, as the ion size increased, 
the reaction rate constants for further addition decreased. They 
estimated the ratios of the rate constants as Ko o/*Ka A 0.002 and 


kK 6/*oa ¥ 0.4 by an indirect method of measurement. The values 


for Koa and ke 5 were then calculated by assuming that ke had the 


Gioumousis Stevenson value of Patil cc/molecule sec. 


, + ; ‘ : : ae 
The ion CoH. , which is the major ion formed on decomposition 


+ * 
of (C,H. ) ,. may react further with C.H 


4°8 oH, and undergo addition 


reactions, 
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ae wate ant ah 


E35 


bse 13} ery Wee Geet cpostang ty" 
a5 24 5°9 
etc 
(5:13) ont + CH, > of lla 
etc. 
+ % 4+ % 
(5:14) Ceeeermronas ACS 
etc. 


+ | per 
These ions ORS and CoH 3 could also result from decomposition 


+ * + * 
of the excited intermediates (CoH, and (CHG ) by reactions 


parallel to (5:8). The rate constants for the reactions (5:12) 
to (5:14) could not be measured by Kebarle et al (42,85). The 
present study was undertaken to measure directly the rate 


constants for reactions (5:6) and (5:10) to (5:14). 


ase Experimental 


; A 4 i ; Te 
The time dependencies of the intensities of the ions CoH, 


4 + = + 
op Nell ana C,H, and C_H 


a : 
CAH, , CoH > 385 7H) 3 were measured for various 


mixtures of xenon and ethylene. For the study of CoH, the 


ethylene pressure was varied over the range 0.74 to a0 


torr,. and the xenon pressure over the range 0.42 to 7.4 torr. 
+ 
CoH, was studied at xenon pressures around 3.8 torr and methane 
—3 +4 
pressures from 1.5 to 4.7x10 torr. The reactions of CoH, 


+ 
and CH. were fast and consequently the ethylene pressure 


had to be low. For the ions which had lower rate constants 


the ethylene pressure was increased up to almost 1 torr 
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while the total pressure was kept constant at about 4.0 torr. 

All runs were made at room temperature (297 to 300°K). 
The pulse widths At. and At, were set at 10 micro seconds 
and the pulse repeat time at 2000 micro seconds so that the 
ion intensity decayed essentially to zero between pulses. 

The experimental work was done before the normalization 
technique described previously (section 4) had been developed. 
Consequently only the intensities of the reactant ions were 
measured and two assumptions were made; firstly, that the rate 
of decrease of intensity of the ion through reaction was much 
faster than the decrease in total ionization by diffusion 
and mass flow,and secondly, that the formation of an ion by a 
previous reaction was completed before the time during which 
the decrease of its intensity was observed. The results 
obtained were similar to those shown in figure (4:3) for the 
methane ions. It was shown in the methane study that the rate 
constants calculated when these assumptions were made, were 


within a factor of two of the corrected values. 


5:3 Results and Calculation of Rate Constants 


The decrease of Cue intensity with time for several 


pressures is shown in figure (5:1). The reaction was assumed 


to be pseudo first order since the concentration of CoH, was 


much larger than that of Sane. The rate constants were 
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calculated by means of the first order equation (5:i) in 
(5:2) Kk & 2,..9033 log ir oMtr 
(tg - ty) [C,H] 


which the ion intensities I, and I, were measured from the log 
2 1 


plots sagurer (5:1). The values for Kee the calculated rate 


+ : : 
constant for the decay of CoH, by reaction (5:6), are given 
in Tablet. 1. 

The ion intensity curves for the other ions are given 


in figures (5:2) to (5:6). The rate constant values, tabulated 


in Tables 5.2 to 5.6, were also calculated by means of equation (5:i). 


TABLE 5:1 
Rate Constant k. for the Decay of oa 
a a a DE 
+ k6 +, % 
(5:6) CoH, ay CoH, (CH ) 

Xe pressure CoH, pressure ke 
(torr) beGut ex 10°) (cc/molecule sec.x 10”) 
0.42 8.4 ate 
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TABLE 5.2 
Rate Constant k,)_ for the Decay of C,H, 
mi em een * (cn, 7) - 
4°8 2 4 G6 IZ 
Xe pressure CoH, pressure | Kio 
(torr) (torr) (cc/molecule sec.x ees 
134 0.035 | (Fo 4 
1 Beg 0713 4.0 
abe 0.28 Ppa, 
20 0.42 32d 
ee 0.64 3.39. 
SINE 0.66 Sau 
ae 2 0.49 4.22 
Sis) 0.70 geo2 
3.6 0.288 See 
356 0.288 3:0 
oy oer 4.5 
ERE 0.078 tee 
S70 23k ae 
4.5 OLS 4.5 
Average value of k = reeleraksie a cc/molecule sec. 
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TASLE—5..3 
4 
Rate Constant kj, for the Decay of C Hy5 
(S421) Seta Feely 7+ ctx, *)" 
: Gnl2 2 4 8°14 

Xe pressure CoH, pressure ki 
(torr) CEorr) (cc/molecule séc.x Tort 
ek 0.38 n He oS 
220 0.56 a2 
3.0 0.42 1207 
nade f OT62 ote | 
ar: 02256 1.24 

-12 a 
Average value of Kia = 1.14x10 cc/molecule sec. 
TABLE 5.4 
+ 
Rate Constant Ki for the Decay of C,H. 
REGIE avleier 25) Bilin gigs de ome Ye 
. i 2 4 sms, 

Xe pressure CoH, pressure Ki5 
(torr) (torr x 10°) (cc/molecule sec. x 10”) 
360) = UR 0.60 
3.8 15 0.62 
3.0 a HP) 9 eer bo" 
i a 4.7 0.90 

-9 
Average value of k SOL cc/molecule sec. 
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TABLE 5:5 
a 
cod 
Rate Constant k,, for the Decay of C.H, 
ee AS Ses ne a 
" 59 2 a4 plas 
Xe pressure CoH, pressure Ki3 
(torr) (torr) (cc/molecule sec. 10) 
; 0.44 0.49 
2.4 0.336 0.60 
- 0.50 0.41 
2.5 0.50 Oo. 39 
oa 0.38 tele 
Ze 0.54 Ono 
Ze Oseed 0.63 
-12 
Average value of Ki3 = 0.46x10 cc/molecule sec. 
TABLE 5:6 
+ 
Rate Constant Kia for the decay of CH) 3 
Gulia nH. at ne Se (ef...) 
. 7e3 214 Prony 
Xe pressure CoH, pressure Kia 
12 
(torr) (torr) : (cc/molecule sec. x10” ) 
1.8 0.36 0.34 
RR 0.294 0.30 
25 0.43 0.21 
oe 0.45 Op2L 
ype | 0.46 Dele 
2635 0.33 0.26 
pp D220 0.44 
-12 
Average value for k 5022 1 eLO cc/molecule sec. 
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5:4 Discussion of the Results 

The measured rate constant Ke was found to be dependent 
upon the pressure of xenon. The rate constant increased with 
xenon pressure as can be seen from figure (p20). This may be 


explained in the following manner. Ke was calculated from the 


ag : 
assumption that the change in CoH, concentration only depended 


. 


upon its loss by reaction (5:6) and could be represented by 
equation (5:ii). However perhaps the reverse of reaction (5:6) 


(5:44) dic,H ey a 


dt 


eae cco Ae Coes 


should be included,as in (5764); 


k 
+ a st as 
: — 
(5:6a) CoH, = CoH, : (C,H, ) 
r 
+ 
The change in concentration of CoH, is given by equation (5:iii). 
ae? + % + = ai 
oe be ® dicjH, ] = k (C,H) [C,H J] + k (C,H, ] 
dt 


: ar? Vite : ; 
The concentration of oud a is unknown. The change in concentration 


+ * + * 
SMG) 49s given by equation (Siiv), since (C,H ) 3s involvediin 


4°8 8 
+% 5 + +* 
(5:iv) dIC,H, ] = k (C,H) (CH, ] - k ICoH, ] 
dt 
+% +%* 
“k, (C,H, ] - k_ [xe] [C,H ] 


reactions (5:6a), (5:7) and (5:8). ke is the rate constant for 


le a ' 
decomposition back to CoH, whereas ks is the constant for 
os + + , F 4g 
decomposition to C_H and ©) H-. If the excited ion (C,H ) 
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is assumed to be in a steady state concentration, i.e. 
+% : : 
alC,H, ]/dt=0 equation (5:iv) may be rearranged so that 
+%* 
(C,H. ] is expressed by (5:v). 


ee y) icuoeeP = k (CH) (cH) 
me 4'8 eke. A 2 ik 
ka + kL + ke [Xe] 


This may be substituted into equation (5:iii) which is rearranged. 


i + Ever a + 
(5:vi) dic,H, ].= KS uf x2 [CHA] [C,H ] 
dt k, + ke. + Ke [Xe] 


Thus the experimentally measured k_ is given by equation (5:vii). 


6 


(Sivii) Ke = kK. La k. 





Ks + k + ke [xe] 


When the xenon pressure is large we may assume k [xe] >> (kj+k)). 
Thus kK. becomes defined by (5:viii) where p is the xenon pressure. 
(Se Voi) k~ =k 
pee 
From figure (5:7) Ke at high Xe pressure is 2.1x10” cc/molecule sec. 
which is the value of k. 
When the xenon pressure is zero, equation (5:vii) reduces 


to -(5<2ix) 


(5'si3) ag 8 ol a4 
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This can be rearranged to give (5:x). 


(5:x) oe vet kk 
a = 
ke ka 

po 


When p?0 Ke from figure (5:7) has a value of about 0,7x10 > (elo J 


molecule sec. Thus equation (5:x) can be used to find the ratio 


of k /Kg equation (5:xii). 


: =9 
(5:xi) ai |, are kK 
0. 7x10 ae 
d 
(5: xii) a = 2 
ca 
* 
This indicates that at low pressures, twice as much (Ce) 
+ + + 
decomposes back to CoH, , aS decomposes to form C,H, and CHW : 


An approximate theoretical prediction of the change of 


the experimental k_ with xenon pressure may be obtained by 


6 
substituting the values for ka k /Ka and KS back into equation 


(52711)... First.a value .of Ke is required. Equation (5:vii) may 


be xewritten as (5:xiii). 


(ov et) k. = k di 1 


If the values are substituted into (5:xiii), k. is given as 


6 


follows. 
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(5:xiv) vr aghs S. PeGe eet, i} 
St} 
1.5 x, [xe] 
k 
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Arbitrarily, the experimental value of Ke (of 1.66x10°) at 


4 torr Xe was substituted into this equation, and the ratio 
k/k,. calculated as eee tees cc/molecule sec. The curve on 
figure (5:7) was calculated from (5:xv). 


(5:xv) ke = Ze loClOh aI! ccedin' 1 
il Gpegpertalal-ey pes 


Since the accuracy of the data is probably not better than a 


149 


factor of two, the agreement of the theoretical and the experimental 


curves is acceptable in that it shows that at low xenon pressures 


a large amount of the excited (C,H 


8 28 


at higher pressures it was stabilized to C)H . Values for Kae 


k. and Ke could not be determined from the present data. A value 


for k from which is and Ke could be calculated, might be 


a’ 


obtained from the increase in intensity of C However in 


yt 
i ie ate 
the present study the increase was at short reaction times 


when the increase in total ionization was rapid. Even if this 


ae + 
problem was not present the rapid reaction of CoH, with CoH, 


would interfere with the calculation of Kae as the rate of 


nF - 
loss of C,H. is almost as large as the rate of formation. The 


ratio kK: calculated as D acwian. cc/molecule from the 


ratios ko/*, and kK, is close to the previously obtained 


+. %* + 
) returned to C._H, , whereas 
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value (42,85), which was estimated as 2 to hes Ciel cc/molecule 


for larger ions i.e. the CE lg Ca homologs. 


Values for Ka have been estimated by other workers (32) 
| -1 ¥ -1 ‘ 
k, = 4x10° sec ~ and (80) k, = 5x10 sec ~. These would give 


values of ke We 2e100- cc/molecule sec. Tiernan and Futrell (83) 


calculated Ke as wy teaion” cc/molecule sec. by the means of the 


Gioumovsis-Stevenson equation (12). Values of 10’ eae for 


. 


* 

k, mean that at pressures of the order of 1 torr (Cn) would 
have a lifetime large enough for a portion of the molecules 
to be stabilized. This was observed previously (42,85,88,89) 
and in the present study (figure 5:7). 

For the following discussion it is assumed that the 
Measured rate constants for reactions of the larger ions 
bea 

ion 


were equal to the rate constants for the addition of CH, 


by reaction il.e. Kio = Koa and Kul = K3a° This would occur at 
the high pressure limit and consequently the assumption may not 
be completely correct as the xenon pressures were not very large 
i.e. <4 torr. The present study confirms the observation by 


Kebarle et al (42,85) that the rate constants for the addition 


of C_H, decrease: as the size of the ion increases. The 


2 4 
+ -9 + 
observed values were, for C_H eo) ce, LOX CH 
2 4 a 4°8 
OMe toe and for Cn.) kc. MP1: i4xi0 (all cc/molecule 
oy et pace ae es abel aga he Py es 


sec.). The ratio ky o/*s = 0.002 is the same as that estimated 


by Kebarle et al (42, ). The ratio KS 55 = 0.26 also agrees 
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with the previously determined value of 0.4. 
It was suggested (42,85) that the large decrease in rate 
constants for the addition reactions must be due to structural 


; ie Sg * 
effects, the largest decrease occurs in the reactivity of CAH 


+ 5 

compared to CoH, (k, /*, = 0.002). One would expect the 
pe Ge ; 
transition state of (CH, ) to have a linear structure. 


+ 
ee eae Hae oe CE 


This could be converted into a stable 1-butene ion by two 


1,2 hydride ion shifts or one 1,3 hydride ion shift. 


HC et —— CH. Sat 


* 
zs) 


Since the (C,H is highly excited a 2-butene might result. 


8 


< ° 
ee, ee 
H.C c CH, 


The reaction of a stable 2-butene ion with CoH, would 


+, % ; 
produce an ion (CH ) with the charge on a terminal CH, group. 


Ze 2 


Since a primary carbonium ion is less stable than a tertiary or 
secondary ion,this transition state would have a high energy, and 


+ 
thus the rate for reaction of the stable CAH, would be expected 


to be low. Kebarle et al (42) measured the cross section for 


the reaction of 2-butene ions with CoH, in a conventional low 


pressure mass spectrometer and observed the value was about 500 


P ‘ : eee 
times smaller than the cross section for reaction of CoH, with 


C,H, under similar conditions. Tiernan et al (83) recently 
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confirmed this measurement. They found the 2-butene ions 
reacted at a rate 1200 times smaller than did the ethylene 
ions. They supported the proposal that the decrease in 
reactivity appears to be due to the tendency of the charge on 
the ion to seeka secondary or tertiary position. The ions 


+ 
H would thus become large branched polymers. 


larger than C, 9 


The rate constants for the further reaction of the ions 


“ + 
CoH. . C.Hy and C_H 


7 fee showed a parallel decrease in value as 


the complexity of the ion increased. The values were CoH, 
12 


-9 + - + 
12 OSLO pp. LOX CoH, Ki3 = Q.46x10 and for CHB, 3 


-12 
14 0.27x10 ° 


k 


k The ratios are thus kj 3/K12 = 0.0006 and 


Ki a/*13 = 0.6. The reactivity of C Fs is almost as high as 


BYs 
that of Stee This was predicted by Kebarle et al (42,85). 


+ 
The CoH. ion resulting from the fragmentation of the l-butene 


structure would probably have the charge on a terminal CH. 


group. Consequently, there would not be a large energy increase 


+ 
CH, = CH Stikhe we 


4+ %* 
in the formation of the excited (CH, ) ain which the charge 


would initially also be on a terminal CH, group, and the 
ey 
reactivity of CoH, might be expected to be of the same order 


+ * 
as that of ony However (CH, ) could rearrange and be 


+. * ; : 
deactivated like (CH, ) to form a secondary carbonium ion 


and thus have a low rate for further reaction. 


Because the rate constants for the reaction of the ethylene 
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primary ions with ethylene were observed to be large,Lampe (90) 
suggested that ion-molecule reactions could play a role in the 
radiolysis of ethylene. Ausloos et al (91) considered the 
contribution of ion-molecule reactions to the formation of the 
small molecular products He and CoH, and concluded that ion- 
molecule reactions were only partly responsible. Wagner (84) 
found that when CoH, was irradiated at -196°C, low molecular 


weight branched polymers were formed. He suggested the 


; ts tdi dee 
polymerization was initiated by CoH » propagated by ionic 


4 
ve 


condensation reactions giving large molecular ions (C Bon ; 


and terminated by ion electron recombination. This is similar 
to the reaction sequence suggested by Field (32) for pressures 


less than one torr. 


+ +.% + % +, % 
> > > 
CoH, (C,H, ) (Cui, ) (CoH. ) ete. 
+ + + 
+ ~ + 
wake CoHy avons 


In this case the chain propogation is by the molecular ions 
ve . ‘ + 
I dd 
CHon (n even), and the carbonium ions Copy (k odd) are 
produced by decomposition. Kebarle et al (42,85) supported 
this mechanism in their study of ethylene at pressures of less 
than 1 torr, in the presence of xenon up to 40 torr. 


On the other hand, Wexler and Marshall (82) proposed 


that, initially the carbonium ion oat was formed by the 
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+ * : : 
decomposition of (C,H ) and then chain propogation occurred 


8 


through the carbonium ions C (k odd). 


nt 
ge 


+ +. % + + 4 
> > - - 
CoH, (CH, ) CH C,H, CoH, 3 etc. 


Since in the present study we found that the rate constants 
for the addition of CoH, to the carbonium ions were similar in 


Magnitude to the rate constants for addition of CH, to the 
molecular ions, it appears that the low pressure (.1 to 10 torr) 
gas phase polymerization of ethylene involves both mechanisms. 
That is, that the chain propogation occurs through the molecular 
ions and, some of these break down to the carbonium ions which 
also propogate the polymerization. Kebarle et al (42,85) 
observed that the rate constants for the larger ions Cy and 


Ch0 were smaller (less than 10ae cc/molecule sec.) than 


those for the C, and Ce ions, reported here, and that the 


high pressure spectrum of ethylene appeared to terminate 

Close to m/e 200. If the rate constant is taken as eree ccs 
molecule sec. and the concentration as Sx10.! molecules/cc 

(10 torr) the half life ee = 1/k [C,H ,]) becomes about 3x10" 
seconds. The average life time of a charge would be about 
1072 seconds. Since these are of the same order of magnitude © 
the polymer size would be prevented from increasing as ions 


would be lost by charge neutralization. At a pressure of 


one atmosphere, fss105° molecules/cc) the half life of a 
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reaction with a rate of foace cc/molecule sec. would be 


about Sn1008 seconds and if the charge lifetime were 10> seconds, 


the polymerization could continue. 
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6. ION-MOLECULE REACTIONS IN PURE OXYGEN AND PURE NITROGEN 
+ 


t- + + 
H Hh ME + ep d a = 
THE ATTACHMENT 0. aus On an N. N, N 


6.1 The Importance of These Reactions to Various Fields 
In this section, studies of the attachment reactions (6:1) 


and (6:2) in pure oxygen and pure nitrogen will be discussed. 


4 + 
: —> 
(6:1) Coe eee 4 
+ + 
: —_> 
(6:2) Nan tallo Ny 


The formation of X and Na have previously been studied 
quite extensively for several reasons: to provide data from 
which theoretical predictions of the structures of the complexes 
could be made, to determine the charge carriers in gas discharge 
and mobility experiments and in atmospheric electricity, and to 
establish reaction sequences occurring in the ionosphere. A 
brief description of the role attachment reactions, such as 
these, play in the determination of ion mobilities and in the 
ionosphere will be given to illustrate their importance. 

The mobility of an ion is determined by measuring its 
velocity as it drifts, under the influence of an electric field, 
through a gas at high pressure (vV 1 torr). Mobility is defined 
as the ratio of velocity to field strength. The mobility is 
dependent upon the mass of the ion and upon the ratio of 
field strength to gas pressure, E/p. In the early experiments 


the ions were not mass analysed and the identities were guessed 
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or deduced from other information such as quantum mechanical 
predictions. Consequently there arose apparently inconsistent 
results. The drift velocities of many ions increase in a 
uniform manner as E/p is changed. However it was found that 
thismaidunot occurs foren 3 in N. gas. Varney (92) proposed 


2 2 
that reaction (6:2) took place as it would explain his results. 
This was later confirmed when the ions from a drift tube 
containing nitrogen were analysed mass spectrometrically and 
a 
Ny was observed. 
Since the mobility of an ion which may undergo attachment 
reactions depends upon the time spent as the original ion or as 
a complex the rates of formation of the complex must be known 
before accurate mobilities can be determined. Knowledge of 
these rates is also necessary to formulate and explain ion 
mobility theories. 
The products of attachment reactions have also been observed 

in the Earth's ionosphere by sampling its ionic constituents with 

; + + 
a rocket-bornemass spectrometer (40). The ions #0 ; HO. and 
ions with masses greater than m/e = 45 were present in the 
D region of the ionosphere. The D region is the lowest region 
of the ionosphere in which there is an appreciable concentration 
of ions and electrons (93). The pressure in this region is 


high (from approximately .01 to several torr) and thus 


attachment reactions, which frequently involve three body 
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collisions become as important as two body collisions in 
changing the ionic composition. The identities of the major 
ions in the ionosphere are required for calculation of ion- 
electron recombination rates. These rates are used to correlate 
the measured electron density with the observed ion densities 
and rates of production from the Sun's radiation. The rates 

of ion-electron recombination reactions depend upon the structure 
of the ions. Since large ions e.g. complex ions, can release 
the energy of ion-electron recombination by dissociation they 
have larger recombination coefficients than do smaller ions. 

For example, at 300°K the dissociative recombination rate of 
NA is a factor of 10 larger than that of No (94). 

Attachment reaction studies provide essential information 
about the structures, stabilities, rates of formation and rates 
of loss of ion complexes in the atmosphere and can provide 
useful data to assist in the construction of ionosphere models. 
These models in turn may be used to predict the effect of 
chemical reactions on changes in electron density. The electron 
density is a basic property of the ionosphere and controls many 
physical phenomena such as the propogation of radio waves. 


+ 


a 
— 
6.2 The Attachment of Oxygen 0. + 0. oO, 


i Previous Work on this Reaction 


Ion molecule attachment reactions of oxygen have been 
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less extensively studied than those of other small molecules, 


especially nitrogen. It is not exactly clear why Ny has been 


studied more than oa but it is probably due to short lifetimes 


of the filament in oxygen. Oxygen is much more reactive than 
nitrogen towards electron filament materials such as tungsten. 
Because of our instrumental design the exposure of the filament 
to gas escaping from the ion source is minimized and this makes 
experiments with high ion source pressures possible. 


+ 
Brederlow (95) first detected the ion O in the positive 


4 
column of an electric discharge in oxygen at high pressure and 
current densities (2.5 torr, 100 mA) by means of a low resolution 
mass spectrometer (resolution 1 in 6). Since Brederlow was 
studying ionic diffusion, specifically the transition from 

free diffusion to ambipolar diffusion as oxygen pressure and 
discharge currents were changed, the kinetics of reaction (6:1) 
were not studied. In another oxygen discharge experiment 
Knewstubb et al (96) observed the ions in the negative column 


: “ ; 
and measured very small concentrations of O In this 


ae: 
experiment the pressure and current were low (0.4 torr and 
0.4 mA) and thus reaction (6:1) was not a major reaction. 

The only previous extensive study of the oxygen system 
was done by Conway and Yang who studied the equilibrium but 


not the kinetics of reaction (6:1). These authors studied 


the formation of O,” by reaction (6:1) (97) and also the 
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+ 
6 to O10 by the general 


formation of larger clusters (98) O 
reaction (6:3). In Conway's apparatus purified oxygen was 


+ + 
° —_—_-> = - 
(67,3) Onn 3 0. Ont? (n is) 


flowed through the ion source and ionization was caused by 

beta radiation from a titanium tritide source. This source 
overcame the problems of using an electron filament. The ions 
diffused towards the wall and passed through two 25 micron 
diameter holes into the evacuated mass analyzer region. The 
ions were then accelerated and analysed. Conway estimated 

that the reaction time was about 16a sec. ae the reaction 
(6:1) had reached equilibrium. The equilibrium constant for 
reaction (6:1) was calculated for pressured from 1 to 9 torr of 
oxygen, over a temperature range of 259 to 344°K. The 
equilibrium constant Kp was evaluated by equation (6:1) 


where Poo is the pressure of oxygen in torr, and joy and jo. 
4 TP 
*O4 
(Orn) Kp = a 
io. e Po. 

the observed ion currents, which were assumed to be proportional 
to the partial pressures of the ions. 

From a van't Hoff plot (of log Kp versus 10°/T) Conway 


and Yang obtained values for the enthalpy and entropy, i.e. 


AH = -9.60 kcal/mole AS = -20.6 kcal/mole °K. 


Theoretical calculations were made of the structure of O, . 
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Figure 6:1 Observed ion intensity curves of 0, and 0, in 0. at 298°K, 


pressures from 0.9 to 4 torr. 
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Three possible structures were proposed which gave values 
of the entropy in agreement with the experimental value. 
The present study was made to measure the kinetics of 


reaction (6:1) and to establish its kinetic order. 


ii Results 

The changes of ion intensities with time were measured 
for pressures from 1 to 8 torr over a temperature range 298 to 
36)°R . Figure (6:1) shows such time dependence curves obtained 
in oxygen at 298°K at pressures from .9 to 4 torr. The total 
. e - : ; = + 
ionization curve was determined by summing the 0. ana the 0, 


7 + : , vey 
intensities. (The O and 0, contributions were negligible 


and not measured). The intensities of oo) and on were then 
replotted fas fractions i of the total ionization,by the general 
procedure which was validated earlier in section (4:4). These 
normalized curves of the same data are shown in figure (6:2). 
Data obtained for reaction (6:1) at several higher 
temperatures are shown as normalized ion intensity curves in 
figures (6:3) to (6:7). The results for several pressures 


at each temperature have been superimposed to reduce the 


number of figures. 


Inspection of the graphs figure (6:2) shows qualitatively 


how the reaction progressed with time. Initially the intensity 


~ , 
of O v decreased rapidly while that of 0, increased,since the 
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forward reaction was predominant. As the reaction time increased 
the effect of the reverse reaction began to decrease the rate 

of the overall forward reaction until equilibrium was reached, 
and then the ion intensities remained constant. As the pressure 
increased the equilibrium went to the right i.e. more oy was 
formed. Also the effect of the reverse reaction became greater 
so that equilibrium was reached in a shorter time. 

The higher temperature runs, (figures (6:4) to (6:7), 
show that the equilibrium shifted to the left i.e. Kp decreased 
with increasing temperature. The smaller rate of decrease of 
ioo showed that the forward rate constant decreased. Finally 
equilibrium was reached at shorter times and thus the reverse 
reaction must have increased as the temperature increased. 

An experimental problem occurred in that at long reaction 
times the total ion intensity was quite small and thus a large 
error is present in the normalized ion” and le This is*’not 
immediately apparent when the data is presented in the normalized 
form, but it accounts for the experimental scatter of the 
results at long reaction times (> 100 U sec.). At the higher 
pressures (figure (6:2) 3 torr and 4 torr) this problem was 
less acute. To reduce the error the signal at long reaction 
times was increased by increasing the width of the ion gate. 
Since the ratio Te was constant at long reaction times, 


increasing the width should not have altered the ratio, as it 
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iii Evaluation of the Equilibrium Constant 

From the figures one can immediately evaluate the equilibrium 
constant, and this was done first so that the results could be 
compared to Conway's. The intensities io, and io, taken at 
long reaction times (> 100 U sec.) where the reaction was at 
equilibrium were substituted into equation (6:i). A van't Hoff 
plot was prepared by plotting log Kp versus 10°/, figure: (6:8). 
The solid line is a least squares line while the dashed 
line represents the data of Yang and Conway (97 ). The results 
are in very good agreement. 

Thermodynamic properties can be calculated from the change 
of Kp with temperature from the following equations. The Gibbs: 
free energy is related to the equilibrium constant by equation 


(6:ii) and to the enthalpy and entropy by (6:iii). These two 
(6:ii) AG = - RT 2.3 log Kp 
(62111) AG = AH - fT AS 


equations can be combined to give the relationship between AH, 


AS and Kp (6:iv). 


(6:iv) As = AH + 2.303R log Kp 
A 


AH was obtained from the van't Hoff plot by setting the 


slope equal to -AH/2.303R. AS was calculated from (6:iv) at 
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T = 298°K. Equation (6:i) gave Kp in units of torr? i.e. the 
standard state was 1 torr. To get AS° in units based upon one 
atmosphere as the standard state Kp in eoreae must be multiplied 
by 760 torr/atmosphere. Since AH is taken from the slope, the 
units of Kp do not affect its units. 

From figure (6:8) the following values were obtained, 
AH = 9.98 kcal/mole, AS = -22.4.cal/mole °K. These are 
very close to those obtained by Yang and Conway (97). 
(AH = -7.60 kcal/mole, AS = -20.6 cal/mole °K). 

The Kp values, measured here, were slightly dependent 
-upon pressure. Figure (6:9) shows the pressure dependence of 
Kp at room temperature (Vv 298°K). One can see that the values 
increased with pressure up to about 3 torr and were then 
constant ian increasing pressure. The first four points were 
taken from the data of figure (6:2) at 298°K and the later 
four points from the data of figure (6:3) at 297°K. The 
points may not agree exactly as one would expect the 297°K 
points to be slightly higher. However one can see that at the 
lower pressures, 0.9 and 1.8 torr, equilibrium has not really 
been established when the ion intensities had decayed to 
small Lee One might speculate that since at low pressures 
the ion signal decayed very rapidly the ions were removed by 
diffusion and mass flow in a time shorter than that required 
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for equilibrium to be established. If 05 and aye were 
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Figure 6:9 Change of Kp with oxygen pressure. 
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removed at about the same rates (by mass flow) it would appear 
that equilibrium had been established. These two values were 

not used on the van't Hoff plot of figure (6:8). At higher 
temperatures the values of Kp were constant over the range 2 to 4 
torr as equilibrium was achieved in shorter times than at 

lower temperatures. Curves at 1 torr were not measured as 


equilibrium conditions were not expected. 


iv Evaluation of the Rate Constant for the Forward Reaction 
Values of the rate constant for the forward reaction were 
calculated from the decay of io,” at short reaction times (from 
the normalized curves). We wanted to establish the order of 
reaction (6:1). The reaction may be rewritten in the general form (6:1a) 


teeta ners ao ven_tn, o Perey 


2 Pah Giaseree 2 


where n = 1 or 2 for the second order and third order cases 


: : “ay ; 
respectively. The change in O is given by (6:v). At short 


2 


. é Sar 
reaction times [0, ] is small and one may assume that the change 


(6:v) d[0,7] 


Sa + n + oe 
at “4 ke10, J [0,1 + k, [0,1] [03] 


: oe 
in 0, is only produced by the forward reaction. If[0,"] = 0, 


equation (6:v) may be rewritten. This equation may be integrated 


(6:vi) a[o..*] 
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Piro t.. to t,) and rearranged. Substituting 2.303 log x for nx 


i 


we get the equation which was used to evaluate ke. 
Bie) ify ) 

O O 
f te vs 2 ne 


JUS) LO 


(6: vii) k, = 


n 
(t. - ak [0.] 


We have also made the assumption that the ion intensifies ape 
and at were proportional to the concentrations [0,7] and Cale 

Preliminary values for ke were calculated from io,” on 
figure (6:2), at the steepest part of its decay by use of 
equation (6:vii). The second and third order values (n=1 and 
n=2 respectively) are shown on figure (6:10). Since the second 
order constant varies linearly with pressure and is zero for 
zero pressure the reaction is clearly third order up to at 
least 4 torr. Third order rate constants were then calculated 
for all the temperatures given in the data on figures (6:2) 
Soa(G= 7) 5 

The preliminary values for the third order rate constant 
were plotted as an Arrhenius plot. The Arrhenius equation (6:viii) 
May be rewritten in a logarithmic form (6:ix) for the forward 


(6:viii) SY Bee ee, 


(6:1ix) log ky = logA - E -/2. 303 RT 


reaction. A plot of log ke versus 1/T will give a straight line 
of slope -E-/2 .303 R. An activation energy E, = -5.2 kcal/mole 


was calculated from the least squares line. This value is rather 
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Figure 6:10 


Third and second order 
versus oxygen pressure 


(a) preliminary values 


(b) corrected values. 





rate constants k 
for reaction (6:la). 
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large for an ion-molecule reaction, especially since most 
second order reactions have zero activation energy. 

Several experimental effects may explain the large 
Magnitude. Firstly, since the ionization pulse had a duration 


‘ =P é ; 2 
CO. LOPMLCrO seconds, 0 was still being produced after zero 


2 


: ‘ ‘ "lia 
reaction time. Theoretically the steepest part of the 10, 


decrease curve should have been at zero time. However since 
om was still produced at this time, the steepest part of 

the curve did not occur until after a time of about 15 micro 
seconds. However at this time the reverse reaction began to 
contribute significantly since on was not zero. Thus the 
measured preliminary values of Ke at any temperature would be 
smaller than the actual values. At higher temperatures (341°K) 
the reverse reaction had an even greater effect than at lower 
temperatures 298°K, since the equilibrium was reached more 
quickly, in 30 micro seconds rather than in 100 micro seconds 
and the values of k¢ would be even smaller. The log k¢ 
(preliminary) versus 1/T plot would then have a steeper slope 
than that of the true values and thus give a large value of E,. 
A lesser problem also occurred which caused the reaction time 
of 15 micro seconds to be chosen. The change in total ionization 
intensity was rapid in the first 10 micro seconds and thus the 


, Siw at 
relative values of 105 and 104 had large errors. 


The rate constant kre was re-evaluated by taking into 
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account the effect of the reverse reaction. The equation 


(Ges), Lor Kp is developed in Appendix I. 


+ 
2, 303m log [ij (14+C) .- TC] 





(6:x) kr = - 


(1+C) | i [0.] 


0 MS ee 5 A, ee 
The constant C = 10, id, taken at long reaction times i.e. 


eq eq 
when the reaction has reached equilibrium. A plot of 


e305 loglig,” (1+C) - C] versus time has a slope of k[O,]" 


fe 


where n=l for the second order rate constant and n=2 for the 
third order constant. Such plots were prepared for all the 
data in figures (6.2) to (6:7). As predicted by the theory, 


the plots are straight lines. Such a plot for the data at 


298°K (figure 6:2) is shown in figure (6:11). The line 


‘ + ; 
should pass through O at t=O. However since O. was still 


2 


being produced during the first 10 micro seconds the line 


did not always do so. This should not affect the value of kr 
as it was taken from the slope of a line through many points. 
The corrected values of ks at 298°K, for the second and third 


order case were plotted versus pressure, figure (6:10b), and 


the third order mechanism was confirmed. 


The preliminary and corrected values for the third order 
constant are shown in Table 6:1. One can see that when the 


effect of the reverse reaction was taken into account the values 


of Ke became slightly larger. 
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TABLE 6:1 


Third Order Rate Constants for the Reaction 


k 
+ i + 
ee) ° 
0, + 20, Oo, + 0, at 298°K 
y 2 30 
0. Pressure Rate constants (cc /molecule' sec.)x10 
torr k,. by (6:vii) k¢ by (6:x) 
0.95 2.28 2n49 
9 2-03 200 
Ba “neu 2.84 
4.0 FSS 2e20u 


An Arrhenius plot, figure (6:12), was prepared from the 
rate constants obtained by the second method and equation (6:x). 
The least squares line gave an activation energy of -2.0 kcal/ 
mole. This is consistent with values obtained for radical 
recombination reactions involving stabilization by a third body. 

Not all the graphs in figures (6:2) to (6:7) could be 
recalculated using the plot of 2.303 loglio,” (14+C) - C] Since 

| 1+C 

in the cases where equilibrium was reached quickly the ratio 
C= in, (eq) ig, (eq) was large and the corresponding correction 
to io,” became large. Some of these plots did not give straight 


lines from which values of ke could be calculated. These were 


discarded. The points used on figure (6:12) are quite 
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representative of the data. 


v Evaluation of the Rate Constant for the Reverse Reaction 
The forward and reverse rate constants are related to the 
equilibrium constant K,which is expressed in terms of concentrations 


rather than in terms of pressure, as is Kp. The relationship between 
kK 
(6: xi) Kune f 
ky 
the two equilibrium constants is given by (6:xii) where Av is the 


change in the number of moles (AV = 1 for reaction 6:1) and R in 


(6:xi4) K = kp (RT)4Y 


appropriate units is 1fos4etore torr cc/molecule °K. The reverse 
rate constant can then be calculated from (6:xiii). ky calculated 


we 
POsxii). ks AVE MOSER CEST Ta 
Kp T x 1.034x10 cc/molecule sec. 


at 298° K and 314° K_ then has the values 2 4AxLOLA: PA i fe ea 
-1 -12 -1 -1 
sec ana 25x10 cc/molecule sec . These lead to a value 


of the activation energy of the reverse reaction of 7.50 kcal/mole. 


A summary of the results is given in Table 6.2. 
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Constants for the Attachment Reaction (6:1) 


EXPERIMENTAL RESULTS 


AH = -9.98 + .5 kcal/mole. at 298°K 
AS = -22.4 cal/mole deg. 
1S tae hate 8 kcal/mole 
BS = 7.5 kcal/mole 
Ke -= p o0x107- > ecefnatenates sec. at 298°K 
= 127x100 nee mee eee sec. at 400°K (extrapolated) 
k, os Pte eae cc/molecule sec. at 298°K 
= Cine Assia cc/molecule sec. at 400°K (extrapolated) 
K = ees es (standard state lcc/molecule) at 298°K 
= 1.94x10 7° (standard state lcc/molecule) at 400°K 
Kp = 3.85x10 ~ (standard state 1 torr) at 298°K 


= 4.7x10"> (standard state 1 torr) at 400°K 
PREVIOUS RESULTS (Yang and Conway) 
AH = -9.60 kcal/mole at 298°K 


AS = -20.6 cal/mole deg. 
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vi Other Reactions, Possible Interference 
; ar + : 
The ions O and OF were also observed. It was possible 
that they might interfere with the rate constant measurement for 


reaction (6:1). Oo” was produced from O, by reaction (6:4). 


2 


4+ 
(6:4) 0. fe > OO tre +70 (ALP. = 18.9. eV) 


: + : , : 
At lower pressures i.e. O to 0.4 torr O was a major ion in 


3 
the spectrum. Its formation has been studied by Curran (99). 
Since he measured its appearance potential as 17.04.05 eV as 
comparec to, 13.61 (R.E.) for nw angi 20 @V (1. Pe.) tor cine 


he suggested that it would be formed from an excited state of 


yg 
07 by reaction (6:5). The energy of an excited state, in 
(6:5) O = =o a8, tt © 


excess of the ae lowest ionization potential (12.06 eV) would 
be required to break an 0-O bond. The appearance potential of 
17.0 eV measured by Curran corresponds quite closely to that 
of the a excited state a (16.8 eV) measured by Winters 
et al (100). 

Curran observed an appearance potential of 17.0 eV for 


SF : 
Oo, also and suggested it came from the same state. One would 


+ 


expect that this ion could be formed from several of the 05 


states Cn, = 12.06 eV, “n = 16.0 eV, on = 16.8 eV or 


4 ba +%* 
ve = 18.0 eV) as the excess energy of Oo, could be 
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removed by collision. The lowest state, however, would be 
expected to contribute the most since firstly it is the 
state of highest abundance (76 Lindholm) and also it would 


+% 
produce O 


4 with the least excitation energy to be removed 


by collision. The excitation energy would be less than that 
required by reaction (6:6). That is, the high energy excited 


“a + . ’ ‘ 
states could produce O and 0, depending upon the lifetime 


4 


+ 
of the excited complex, but the lower O could only produce 


2 


+ 
oO, . 


: a = aaa pe 
In this work the observed O and 0, intensities were 


extremely small at oxygen pressures Prot O75 (tOrr, waOourt 


+ + , Sans , : : 
Pico es Of The 0, and O, intensities under constant irradiation 


conditions. An attempt was made to measure oO” and O, but their 
intensities were at the lower limit of detection under pulse 
irradiation. Consequently we assumed that (Rx 6:6) did not 
affect the intensity of oa and accordingly only measured the 
oat and chee intensities. 

In a few of the runs an impurity ion m/e 50 (0, (H,0)") 
was observed. This was due to water. In such cases a new 
oxygen sample was prepared. The mass spectrometer gas handling 


plant and ion source were baked out until the water was removed 


such that only the oxygen ions were observed. 
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vii Discussion of Results 

It was very gratifying to get values for AH and AS very 
close to those of Yang and Conway (97). The latter were corrected 
for the effects of ambipolar diffusion as it was felt that this 
did have an effect upon the rate of ion removal from their type 
of ion source. The values of Kp calculated by them were not 
independent of pressure. A close examination of the van't Hoff 
plot in reference (97) reveals that the values increased with 
increasing pressure. Because of this Yang and Conway felt that 
the results obtained at the low pressure limit of their apparatus, 
1 torr gave the most meaningful results when corrected for the 
different diffusion rates of the ions. 

The equilibrium constants calculated from our data were 
constant above a certain pressure of oxygen. The ions were 
probably removed more by mass flow than by ambipolar diffusion. 
In this case there would be no discrimination in the sampling 
of ions and the values obtained for Kp should be correct. The 
ion source exit leaks used by Yang and Conway had a total area 
which is about 1/15 of the area of the exit slit used on the 
present instrument. The outflow of gas would have been smaller 
by this factor. It is, therefore, quite possible to have the 
ion extraction diffusion controlled in the former and mass 
flow controlled in the present apparatus. 


It is usual to equate the AH for a dissociation reaction 
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to the bond energy of the dissociated bond. The AH value 


of the attachment reaction is the reverse for that of 


a 


dissociation. We can then write for the 0, ion: 


(6:xiv) =hHE = D(0, -0 ) 


2 
Thus the dissociation energy of On es° 1OVO. kcalen Alsodet 


the change in volume is small one can equate AH V AE and thus 


: % ps 
(6:xv) AH E. E. 


-10 vy -2 - 7.5 = 9.5 kcals/mole 

(The change in volume PAV gives the difference 0.5 kcal/mole). 

It was observed that the rate constant for the forward 
reaction of oxygen was approximately thirty times smaller than 
that of the corresponding reaction in nitrogen (101). This will 
be discussed, along with the energy transfer theory, in the next 
chapter. Since data indicates that the reaction is third order 
we may then write the reaction as below (6:1b). 


k 
+ f + 
A 
(6:1b) 0. + 20, EGA 4 2:2, 


The two rate constants ke and ke may be composite rate constants. 


= + 
633-7 At Oleh N Si A 
he tachment o itrogen, N, N., Ny 


i Previous Work on this Reaction 


A reaction which parallels the oxygen reaction studied above 


‘ Ty : F 
is the attachment of No molecules to N, ions. This is represented 
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by reaction (6:2) 


+ + 
: _—> 
(6:2) N. + N, Ny 


This reaction is considered to cause the low mobility of No ions 
in nitrogen gas and it has been studied extensively with respect 
to mobility problems (102,103,104). The equilibrium constant and 
the rate constants have been measured at energies higher than 
thermal in electric drift fields and at thermal energies in the 
absence of an electric field. 

Varney studied the mobility of No” and Ny in nitrogen in a 
drift field apparatus (92,105-108) operated at pressures up to 
35 torr and field strength 200 volts/torr. He calculated the 
extent of dissociation of Na from the change of ionic mobility 
with the eet ion, of E/P, (electric field strength/normalized 
pressure where 25 = px273/T). From the extent of dissociation 
he deduced a value for an equilibrium constant K. He also 
proposed a relationship between tie ion temperature T, and 
the E/Po ratio (T a Ties + ak /P, where a is a constant). 
A van't Hoff plot (2nK versus 1/T,) was used to calculate a 
value for the dissociation energy D(N, -N,) = -AH. His value 
for AH was -0.50 eV (VV 11.5 kcal/mole). 


An apparent equilibrium constant was obtained under thermal 


conditions at a pressure of 1 torr N 


> by Knewstubb (109). His 


value was 75 (standard state 1 torr ) at 300°K. If Varney's 
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van't Hoff plot (108) is extrapolated to a value of 1/T, = 1/300°K 
the Kp value obtained is about sitor? This is extremely large 
compared to the values previously obtained for the oxygen reaction 
(e21)*: 

A value for the dissociation energy can be calculated from 


ob 
another set of data. N, is also formed by a second reaction, 


. 


the Hornbeck-Molnar process (6:6). Asundi, Schultz and Chantry 


+ 
(103) measured the appearance potentials of N, and No in nitrogen 


+ 
(6:6) N +N it +e 


by means of a fer een ease electron beam technique. They took 
the appearance potential of No to be 15.6 eV as the reference 
point of the electron energy scale. (The spectroscopic value 
is 15.58 eV. Saporoschenko (110) obtained a value of 15.54.2 eV). 


“f 
The appearance potential of Ny was 15.04.1 eV). This value is 


* 

the energy of the excited neutral N, . It was also measured by 

Curran(11ll1) as 15.04+.5 eV. The difference between the energies 
+ 

of the two methods of formation of Ny gives a minimum for the 


bond energy DIN,” - N,)- This is given by (6:xvi). 


cca pes hes teh hel A a eG OS 


(6: xvi) D(N., 2 


] eilae0 ti2.s KCals/mole 


The two values of -AH, 11.5 and 13.8 kcals/mole agree within 


the experimental error of 0.1 eV (2.3 kcals/mole). 
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The rate constant for the forward reaction (6:2) has been 
measured under drift field and thermal conditions. Since the 
kinetic order has been under contention(112), we will write the 


reaction as (6:2a) in accordance with the reaction (6:la) for 


oxygen. 
a 
(6:2a) N + nN oa IN + n-l N n= 1,2 


Nei toOr a. second order reaction ena n=2,. Lor a third order reaction. 


by means of a drift field technique 


29 


Warneck (101) measured Ke 


and obtained a third order value of 8.5xl0~ ap pele. secs 


The pressure was varied from .09 to 0.2 torr. The constant 
appeared to be independent of E/P, as well as pressure. McKnight 
et al(104) measured drift field velocities of the nitrogen ions 


+ + + + 
N , N, , N, and Ny at pressures trom, ).o. to 1.0, torr: and 


calculated a third order constant for reaction (6:2a) which 


was dependent upon E/P. Their values ranged from son open at 


49 


low E/Po to 2x10 ofpu gto Prasunssecyes at high E/Po: 


Knewstubb (109) reported a third order value of a meet 


Be Vaclecute. sec, measured under thermal conditions at pressures 
up to l torr. He reported a temperature coefficient of "close 

to -l percent per degree at 298°K." This gives an activation 
energy of about -3.6 kcal/mole. Fite et al (45) observed the 


formation of N,~ by means of a pulsed afterglow technique. 


4 ’ 


They estimated a second order k, value of approximately 
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Santen cc/molecule sec. and this corresponds to aiGane pea eiecn ea 
sec. (third order) at the pressure of 0.5 torr used. _ More recently 
Fehsenfeld(113) reported a reaction (6:7) similar to (6:2a), 

observed by means of the flowing afterglow technique. Since the 


pressure of helium (about 1 torr) was approximately 100 times 


+ + 
: > 
(6:7) Ny + N, + He Ny + He 


that of nitrogen, helium acted as the third body. Rate constants 


oe nea male tee ee: at 280°K and elec 


were reported, 1.9x10- 
2 2 2 : i 
cc /molecule sec. at 80°K. The activation energy can be 


calculated to be -0.37 kcal/mole from these two values. 


ii The Purpose of This Study 

We decided to look at the nitrogen reaction under the 
same conditions as used for oxygen for several reasons. Firstly, 
the measurement of the rate constant for the forward reaction 
in nitrogen would provide a check on the value for oxygen. It 
has been noticed that these values were different by a factor of 
approximately thirty which would appear to be exceptionally 
large. The order of the reaction occurring in the pressure 
region 0.2 to 5 torr was under contention. Secondly both 
Varney's values for the equilibrium constant Kp and Warneck's 
value of the rate constant ke for the forward reaction were 


calculated from ion intensities which had been measured 


under the influence of an electric field, in which the ions 
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were not at thermal energies. In the present work, the ion source 
was essentially field free and thus the ion temperature should 
have become the same as the gas temperature after the ion had 
undergone several collisions. Knewstubb's value for the equilibrium 
constant was very different from that obtained by an extrapolation 
of Varney's data (75 compared to Bx10>)*. We hoped to check these 
data. 

We also hoped to obtain the temperature dependences of both 
the rate constant k_ and the equilibrium constant Kp so that they 


£ 


could be compared with those obtained for the oxygen reaction. 


aa : ' oF + 
iii Other Reactions Involving N. and Ny 


wa de : : : 
Ny. ions are formed in nitrogen by two mechanisms. At low 


3 ‘ fe : 
pressures N, is formed from an excited state of N. , (reaction 
6:8), whereas, at high pressures, it is formed by a third order 


+ 
reaction from N , (reaction 6:9). The relative effect of reaction 


+ 
(6:8) compared to reaction (6:2) on the removal of N, ions depends 


+ 4. + + 
: > 
(6:8) N, ( qT ) + N, Ny + No 


+ + 
: > 
(6:9) Noy oN EE IRE a 


+ ; 4.. + ‘ 
upon the amounts of Ny in the excited state LY and in the 
ground state. Asundi et aZ(103) showed that the excited ion is 


deactivated by collisions with N, molecules at pressures above 


10° torr and thus one would expect that reaction (6:8) would 
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not be important in the pressure range 0.4 to 3.6 torr used 
dir wits wore. Reaction (6:0jewill be considered latergin 
section’ (6:5). 


One would not expect that N 


+ 
3 ions would be formed by 


+ . 
collision with N' as the reaction (6:10) is endothermic 


(6:10) N‘*+ N > N+N 


+ a i 4 
(R.E. N = 14.54 eV) ome N, = 15.58 eV). Since the concentration 


of N atoms would be small the reverse of this reaction is probably 


insignificant in depleting the N, concentration. 


4 
The ion N, may be formed from excited N 


4 molecules by the 


2 


Hornbeck-Molnar process reaction (6:6) as mentioned earlier. The 


importance of this reaction compared to reaction (6:2a) depends 


* 


+ * 
upon the ratio of N, toN, . Asundi et al showed that N, was 


deactivated by collision with N, at pressures greater than 107 


2 


torr, and that when the electron energy was greater than 17.5 eV, 


+ 
the production of N, was essentially only by the third order 


+t 


: = ; ; 
reaction from No (reaction 6:2a, n=2). In the present experiment 


the gas pressure and electron energy were high and the contribution 


+ ; 
of the H-M process to Ny formation was expected to be small. 


iv Results 


4 + + 
ard and N 


: + 
The intensity curves for the ions N , N 3 4 


2 


were measured for pressures ranging from 0.4 to 3.60 torr and 
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temperatures from 300°K to 420°K. Most of the results were 


done at pressures 0.6 to 1.50 torr. In some of the runs, only 


+ 
N i and N were measured. The results for 2.65 torr N. at 


2 4 2 


380°K are shown in figure (6:13). The ion intensities replotted 
as fractions of the total ionization are shown in figure (6:14). 
Similar plots for different pressures at the same temperature 
are shown in figures (6:15) to (6:17). 


; ee + aR 
The sum of the fractional intensities of No and Ny LS 


shown on figures (6:14) to (6:17). The sum remained reasonably 
constant as the time increased. This was found to occur in 


many of the runs in which all ions were measured. The sum 


a+ Be, t+ , : ee 
(i, + 1 ) also remained essentially constant. Deviations 
3 


from constancy did not exhibit any trends. It was therefore 


‘ J * ; 
concluded that the ions No and Ny reacted independently from 


+ + 4 
the ions N- and N, ; that N. was lost only by reaction (6:2) 
+ 


to give Ny 


and that the concentration of Ny in the excited 
state ee was low. The intensity curves obtained when only 

+ + , , 
N, and Ny were measured are shown in figures (6:18) to (6:21). 


When there were small amounts of water or oxygen 


+ + + + 
Hee Ono end. 0 were 


i i ti : 1 N 
impurities present the ions 2 2 3 2 


observed. Since all the nitrogen ions have ionization potentials 
+ + 
much larger than the impurities Ho = 12.61 eV) Loe = 12.06 eV) 


the impurities would remove nitrogen ions by charge exchange. Much 


effort was expended to remove the impurities. If the impurity 
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Figure 6:13 Observed ion intensity curves of N No Ny Ny 
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ion intensities were large the results were discarded. In 


* : 
figure (6:14) one can see that oO. ion appeared slowly with 


time. It is mainly responsible for the decrease of the 


; F 3 : : 
normalized intensity of N at reaction times greater than 


4 


50 micro seconds. The oxygen concentration was very low, less 


than 10 4 times that of nitrogen. 


v Evaluation of the Rate Constants 


The forward rate constant Ke of reaction (6:2a) was 


alae : 
calculated from the maximum decrease of in, from the normalized 
curves (both those which were expressed as fractions of 


ay des Sa : ae : 
(i +i ) and those as fractions of the total ionization). 


N2 Ng 
The equation used is similar to that used for oxygen (6:vii) = 


(6:xvii) where n=l for the second order reaction and n=2 for 


mi Pi aus + 
(6:xvii) ke = 2,.50348.00 (in, t J in, ts) 


n 
Le oeae [N.] 


the third order reaction. 

The order of the reaction was first determined from the 
data at 380 and 381°K, since at this temperature the data was 
measured over the widest pressure range. Figure (6:22) shows 
the values of ke calculated for n=l and n=2. The values for 


the second order rate constants increased linearly with 


increasing pressure. The third order values were constant 
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within the experimental scatter and thus reaction (6:2a) is 
third order up to 3.5 torr of nitrogen and may be written as 


follows. 
(6: 2b) N + 2N,, aN +N 


Third order rate constants were calculated for all other 


temperatures. An Arrhenius plot (log k, versus 10°/T), figure 


E 


(6:23), was prepared. The line is a least squares fit to the 


data it gives a value of k, of peepee molectie: sec. at 


3 


298°K and an apparent activation energy Ee = -2.3 kcals/mole. 


These values are quite close to those of Warneck ke = 8.5x1072" 


29 


and of Knewstubb, k, = 8.0x10 ~~, E. = -3.6 kcal/mole. 


re 
It.can be seen from most of the figures that in,” decayed 
to very small values and that the effect of the reverse reaction 
was negligibly small. Thus the rate constants calculated by 
equation (6:xvii) should then be close to the true values. It 
was found that at low temperatures (300°K) the results at 
pressures less than 0.6 torr were unreliable. No ions were 
extracted at pressures less than 0.28 torr and consequently the 


results taken at pressures around 1] torr are the most reliable. 


The technique of correcting for the reverse reaction was 


207 


tried for a few of the results. A plot of 2.303 log lin t (1+C) -C] 


1+C - 


versus t, (where C = in * (eq) /ig TF (eq)), is shown in figure (6:24) 
4 


2 
for the data at 300°K from 0.63 to 1.28 torr. The corrected 
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values for ke are shown on figure (6:23) as squares. It can 
be seen that they gave essentially the same values for ke at 


298°K and do not change Ee This type of plot for other data 


was unsuccessful as the ratio C was extremely small. 


vi Evaluation of the Equilibrium Constant Kp 

The measurements of Kp, from which a value of AH might be 
determined, were unsuccessful. The values of Kp calculated from 
in, and in, at long reaction times from figures (6:14) to (6:21) 
were very dependent upon pressure and also not particularly 
reproducible. Under the low et conditions the total 
ion. signal at reaction times greater than 100 micro seconds 
was small and this contributed to the lack of reproducibility. 
Attempts were made to obtain more sensitivity by setting the 
delay time t" to 100 micro seconds and the ion gate width At. 
to 200 micro seconds. This method gave slightly better results. 
The Kp values calculated from the ion intensities at reaction 
times of 100 ee seconds were found to increase rapidly with 
pressure (figure 6:25). It appears that at low pressures the 
total ion signal decayed in a time (100 micro seconds) much 
shorter than that required for equilibrium to be established. 
Thus the small ion signals detected at low pressures would be 
those of ions not at equilibrium. Consequently one would 


expect the values of Kp to increase with pressure, since the 
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extent of reaction would increase. However at a certain point, 
when equilibrium would be established in less than 100 micro 
seconds the values of Kp would become constant. This behaviour 

was observed for the measurement of the oxygen reaction equilibrium 
constant (see figure 6:9). However it was not observed for the 
nitrogen reaction. This was probably due to the large value of 

Kp. If one takes Varney's extrapolated value of Kp of Exd0" at 
300°K one can calculate the ratio of in, to in, at a pressure 


of 1 torr from (6:xviii) as 0.00001 to .99999. Unfortunately 


(6:xviii) RD) =. = O10 


the 2S of the mass spectrometer was not great enough 
to detect differences of this magnitude from the very small ion 
currents. The limit was perhaps 0.001 to 0.999 and thus the 
largest values of Kp measureable were 107. 

A calculation was carried out to see how the Kp calculated 
from ion intensities at 100 micro seconds would change with 
pressure, and at what pressure equilibrium would be established in 
less than 100 micro seconds. A theoretical Kp' was calculated 
by substituting the appropriate equations for in 35 and iy y into 


7. 4 


(6:xviii). The resultant equation is (6:xix). The value for ke 


was taken from the results figure (6:23) ky was calculated from 


Varney's Kp (5k0° at 300°K) and k- The value of Kp' was found 
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300 °K fk 38 


2 
Ri (torr) 


‘ 7 : 2a 
Variation of Kp of reaction (6:2) N + N 


Liege <Kp 


0°K 


2 a oan 4 


versus pressure. Solid lines experimental. 


lines calculated from equation (6:xix). 
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24°3 
2 -(k [N.]. + k_IN 17) t 
k [NJ] E - e re 2 £ 22 
£5 = 3c100) Kp'= 
"Ny fe, ON) + k,1N,1° ne eared K gIN1°) 4 

to increase rapidly up to pressures of 1.2 torr after which it 
remained constant at Bxc10Ry Part of this curve is shown in 
figure (6:24). The calculation was repeated at 380°K with 
appropriate values of ky and k This is also shown. One 
can see that the theoretical values for Kp at any pressure are 
larger than the experimental values. This is most noticeable 
in the 380°K curves. One might assume that at the higher 
pressures, 3 to 4 torr, the reaction was close to equilibrium. 
That is, that the equilibrium constant at 380°K is larger than 
could be observed with the present instrument but perhaps 
smaller than the value obtained from an extrapolation of 
Varney's results (107). 

An attempt was made to measure Kp, by having a low pressure 
of Nos and a high temperature so that the ratio N,/N. was 
within the range of the instrument detection (not less than 
0.01), and yet have the total pressure high so that ions 
would be extracted efficiently by mass flow. Neon was added 
to increase the pressure. Neon was chosen since its ionization 
potential is greater than that of nitrogen (I,Ne™ = 21.56 eV, 

+ 


a 
Polo = 15.58 eV). It was expected that No ions would be produced 
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214 
by charge transfer and be essentially the only ions produced 
2 * 
(as aa = 24.3 eV > Tene ). Argon was not chosen even though 
+ : + 
Ar has an ionization potential slightly higher (15.8 eV) than Noy 


ae 
as Knewstubb (109) had observed the ion ArN. in mixtures of 


Argon and nitrogen at high pressures. Kaul and Fuchs (114) 


proposed this was formed by a Hornbeck-Molnar Process Reaction (6:11). 


+.) 


* / 
(6:11) Ar + No ms AYN, +e 


Munson et al (115) suggested that the following reaction (6:12) 


would also occur in N, - Ar mixtures. 


+ 


+ 
(6:12) AXN, +N, via Te Ny + Ar 


+ 
It was felt that the presence of ArN, would complicate the kinetics 


2 
and hence neon was used. 

Figures (6:26) and (6:27) show two typical runs of the 

; tae 
N,-Ne mixture. Apparently the Ne ion was unreactive towards 
N. as one can see that its intensity changed only slightly. In 
retrospect, one could explain this by means of the Massey 
adiabatic hypothesis (80) (section 4:4); that the cross sections 
of charge exchange reactions are small for large differences in 
ionization energy ( 6 eV in this case). Figure (6:27) shows also 


that impurities such as O, became the major ions at long 


2 
reaction times, and thus obscured the results. 


The new results for Kp showed no improvement. These are 


shown as filled-in points on figure (6:25). Thus it was impossible 
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Figure 6:26 Normalized ion intensity curves for N,-Ne mixture at 


297°K, total pressure 4.0 torr, nitrogen 0.73 torr. 
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Figure 6:27 Normalized ion intensity curves for N.-Ne mixture 


at 376°K, total pressure 4.1 torr, nitrogen 0.84 torr. 
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a} 


to measure Kp under the pressure and temperature conditions 
available in this instrument. One concludes that the value 
for Kp must be almost as large as that calculated by 

Varney even though his ion temperature scale may be somewhat 


uncertain. 


vii Further Evaluation of Ke 


The data of the N,-Ne mixtures was used in an attempt to 


calculate rate constants for the forward reaction (6:13) in 


+ + 
e > 
(6:13) N, + N, + Ne N, + Ne 


which the third body was not N The equation used was (6:xx). 


Os 


2.303 log (i, ee sitet git) 
2 t, 2 ty 
(6: xx) ke = 


(t,, ~ t,) [N.] [N. + NA] where [Ne]>[N.,] 


Three results are shown as solid points on figure (6:22). The 
results are of the same order of magnitude as those in pure nitrogen. 
: 3 : , 
The two points at 10 /T = 2.6 were taken from runs in which the 
. ; , ‘ f + : 
concentrations of impurity ions 0. and HO were large. Since the 
+ 
rate constant was based upon the decay of i the values are 
2 


probably high since there would be a contribution from charge 


transfer from No” to 0, and H,0. 


29 


One might refer to the value 5.9x10_ eof /molecule” sec. 


obtained at 297°K to propose that Ne is a less effective 
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third body than is N,- The value of 1.9x10 7 cc/molecule sec. 
measured by Fehsenfeld et al 280°K for the rate constant for 
reaction (6:7) in which He is the third body would tend to 
support this proposal. 


A summary of the results is given in Table 6:3. 


TABLE 6:3 


Constants for the Attachment (6:2) 


EXPERIMENTAL RESULTS 


ke 8. 3x10 7? cc” /molecule’sec. Bte29o°R 


= cpa (app BRe aetna ane: at 400°K 
E = -2.3 kcals/mole 
Kp >> 107 (standard state 1 torr) at 298°K 


> 107 (standard state 1 torr) at 380°K 


PREVIOUS RESULTS 


Warneck (101) 


k = Ss AROE ees Fic Maen ene at 298°K 


f 


Knewstubb (109) 


ke = Bip from re eran Sen restated at 298°K 


Ey < <-3.6 kcal/mole. 
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6.4 Comparison of Oxygen and Nitrogen Reactions 
i Theoretical Considerations of the Reaction Mechanism 
The effect of Temperature and Bond Energy on Ke 
The ion-molecule attachment reaction, which requires third 
body stabilization is very much like the recombination of two free 
radicals in the presence of a third body. The recombination has 
been discussed in terms of two mechanisms, the energy transfer 
mechanism and the radical-molecule complex mechanism (116). In 
both, the formation of an excited bimolecular complex occurs . 
In the energy transfer mechanism the complex is _ formed from 
the two radicals and it is stabiriieed by collision with the third 
body M, as is shown in reactions (6:14) and (6:15). The radical- 


* 
molecule complex mechanisms involves the formation of RM which 


reacts with a further radical to form R, and M. 


(6:14) Rat eR R,* 


(6:15) RK oa. Sr at 


The two mechanisms were first proposed by E. Rabinowitch 
in 1937 (117) and have been investigated by many other workers 
(116,118) . 

When the third body is the same species as the reactant, 
the two major mechanisms become the same, and may be treated by 
the theory of the energy transfer mechanism. 


A mechanism for the ionic attachment reaction may be 
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written as shown in the reaction sequence (6:16) and (6:17). 


(6:16) 0," . 0, oo Ses (combination) 
? (decomposition) 
+* 52) ang (stabilization) 
(6:17) oO, + 0, = O, vi 0, 
F (activation) 


The sequence for oxygen is shown. The same reactions apply for 


nitrogen. 


One may find equations for the rate constants by applying 


+%* 
the steady state assumption that the concentration of 0, is 
: +* } +* . 
constant i.e. alo, ]/dt = 0. The concentration of O, is then 


given by equation (6:xxi). 


+e + + | 
(Gtxxipe [oO ] = k ,[0, ] [0.] + kK [0, ] [0.1] 


4 
+ 
Ky k [05] 
* 
The rate Re of the forward reaction is given by the loss of ag 
+ 
to form 0, . 
<: 7% 
(@FXRiT) Re S Ke [0, ] [0] 
Sta if) + Z + 2 
(6: xxiii) Bestar tea LO vale (00). 084 Je ak 25 [0°81]: [02] 


ka + Kk. [0.] 


Similarly the rate R. of the reverse reaction is the rate of 


+* + 
decrease of 0, to form 0. ‘ 
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(6: xxiv) R. = k,[0, 
iaccpiths ) Ra TO al? (0. b& kak toot) (623 
ee re a eae 22 2 dita‘ 74 2 
+ 
Ko - [0] 
+ 


™ ; 
When [0, ] = 0,R, is given by (6:xxvi) and when [0 To] OO; 


£ 


is given by (6:xxvii). 


ri i ey 2 
(6:xxvi) Ry = Ko kK [0., ] [0.] 
ka + kK. [0.] 
© n + 
(6:xxvii) R.. = kK. Ks [0, ] [0.] 


ks ie Kk. [0.] 


Now if ka >> Ke [0], the low pressure situation, the concentration 


* 


of O i is controlled by Ko and k 


4 The forward reaction is third 


a’ 


* 


order with rate given by (6:xxviii) and the reverse reaction is 


second order with rate given by (6:xxix). 


(6:xXxviii) 


Ds) 
i 
nm 
a 
7 
Yn 
O° 


2 , 
] [0.] ik, = kK. k 





(6:xXxix) R 


I 
nw 


+ 
kK. [0, ] [0.] i Kk. 


- : +* 
Conversely, if the concentration of O 


4 is controlled by Ke and 


koe 126 bu 35s kK. [0.1], the forward reaction is second order and 


d 


the reverse reaction first order with rates given by (6:xxx) and 


(6:xxxi). 
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(6:xXxx) R 


I 
a= 
° 
Nh 
ob 
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NO 
axe 
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C6 soc) R 
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nw 
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Under the conditions of pressure and temperature used in 
the present study, the reactions displayed third order kinetics 
and thus Ko >> Ke [0.1]. The rate constants are given by (6:xxviii) 


and (6:xxix) . 





(6:xxviii) Ke = kK. kK. 
oF 
(6:xxXix) k = k 
a a 


Since theories for the effect of temperature upon Ke and Kk. 
have not previously been developed, it might be possible to geta 
qualitative prediction by taking simple equations for koe ar, 

Ss 


d 


ani. 
a 


The rate constant kK is that of the combination process of 


reaction (6:16). 


k 
+ +% 
(6:16) op haart ——- 0, (combination) 


Since this involves the collision of an ion and a molecule to form 
an excited complex one might use the Gioumousis-Stevenson equation 
(l:vi) to predict values of koe 


(1: vi) k= 2mefa 
u 


This equation does not predict any effect of temperature upon ko 
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The rate constant Ka for the decomposition of 0,” in 


reaction,.(6:16).is similar to that of a general unimolecular 


+ RG. 4k 


(6:16) 0.5 + 0. ——- 0, (decomposition) 


reaction, and one might approximate it by application of the 
Rice-Ramsperger-Kassel model. The simplest R-R-K equation is 


given by (6:xxxii). The R-R-K model assumes that the molecule 


te r-l 
° ii) = - 
(6:xxxii) Kg 24 : E 





contains r coupled oscillators all of the same frequency. The 
probability that the total energy E may be transferred into one 
of the oscillators with a critical energy E* is given by the 
term ae) Sa The factor 24 is a proportionality constant 
related Fs the frequency. The critical energy E* may be 
approximated to the bond energy D(0,"-0,) and the total energy, 


to the bond energy D plus some contribution from the thermal 


energies of the reactant species. E is given by (6:xxxiii). 


(6: xxxiii) | eh Ie oy 4 Bal 


It is hard to assign a value to the thermal contribution nRT, but 
since the vibrational levels of oxygen and nitrogen molecules and 
ions have low populations at the temperatures used here (300 to 
400°K) we might consider only the translational and rotational 
levels (i.e. nRT = ans RT + ce RT). The R-R-K equation may 


be rewritten for Ka as -(6:*%xxiv). 
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(6: xxXxXiv) ke =p 7 
: ¢ D + nRI 


Since D > nRT, we may neglect the thermal term in the denominator. 


Equation (6:xxxiv) indicates that ka 
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would increase with increasing 


temperature and decrease with increasing bond energy D. Increasing 


the number of oscillators r which could contribute energy would 


also decrease ka Since (nRT/D)< 1. 


The rate constant kK. is for the stabilization step of (6:17). 


(6:17) O 5 ete C) + O (stabilization) 


; : tie on, ‘ : : 
Deactivation of 0, will only take place if the kinetic energy 


tA é ; 2. ; 
of O and O. is less than a certain critical energy Ee Husain 


4 2 
and Pritchard (119) gave an equation (6:xxxv) for such a reaction. 


They used it for the recombination of iodine atoms. The energy Eo 


(6:xxxv) kK = pa. (1- Pot PS 


was small, less that 1 kcal/mole. One might apply a similar value 
Fomureaction «(G:17) . 2 is the collision frequency and p the 
efficiency of each collision for removal of energy. The effect of 
an increase of temperature would be to decrease ke: Since E is 
small, much less than the bond energy D one might not expect it to 
vary with D. 

All these approximate equations may be combined so that k 


£ 


is defined by equation (6:xxxvi). 
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(6:xxxvi) k = kK PZ. (1 - 


The effect of temperature change would be as follows. Since Ko 


increases, k 


is independent of temperature, K decreases and k £ 


d 


must decrease with increasing temperature. This then agrees 
qualitatively with the observed data, as the Arrhenius plot 


gave an apparent negative activation energy for kee 


One may compare the values for k. of the oxygen reaction 


£ 
(6:lb) to that of the nitrogen reaction (6:2b) by comparing the 


+ 4 
bond energies of the complex ions 0 anc No, (D(0,"-0 


4 4 alk 


2 
+ 
10 kcals/mole D(N,, -N,) ~ 14 kcals/mole. As the bond energy 


D is increased, equation (6:xxxvi) predicts that k. would increase, 


£ 


since ka would decrease by a factor (Do,/Dny)* 


Equation (6:xxxvi.) may be used to qualitatively predict 


how the rate constant ke would change as the third body was 


changed. ke and k, are independent of the third body. Ke 


d 
would be affected by the collision efficiency factor p (equation 
6:xxxv). One might expect that if the third body were an atom 
such as Ne or He, instead of a molecule, that p would be smaller. 
This qualitatively agrees with the findings for the nitrogen 
reaction. At 300°K Ke for reaction (6:2b) from the present 


results was about al Oe’ cote eanles sec. When the 





© 
i ‘ 
m yt? ; 
Tok 
f | 
[~ fon ™ 
ae | 
é » & j 
, «s eT 
i 2744771, BS : : 
- ‘ fy ‘ x 
} 
af 1% ] o git i : 
n e4 +e | t 7 
4 ~ = 
_ my 
' ‘a 
+ 
m4 1aP YAO rp 
i 
‘ 
_ an fei 
i M { 
, ; : 
( £ svi?) 
Ve a . | 
s 
74 : & xc) : 
~ = 
binge, 
“a 
‘ 1. tf } 
poibeiqg vievigse> 1 or B 
» ybod mils ops ‘ 
1 
iy od Dri ef> 10 


ASpoltine gay Wor. ee i edd rade! aawrpe J: 


sacgeta bit 


= 


4 - 

J J) 
+ 

mons (1S 

as 

eae Te 


fog 


a 


otsew yBod Sxire 203 


,sluselom « 26 beevend mek 


od Elvew 4 Jett 


i 


~ "7 





iil 


* stazotow\" 


os on 


‘ 


hee 


east AL pale ee elog\ete 2 


“/ 


+ 







Afi 1gsioqued te eae: a: 


utsxeques to i qebneqsbn 
= 
putwesxont doiw pasardeb 4) 


_ 
_ 
a 
a 
ut 
Saal 
ant 


bevisede a? dow Vewidadt fs up 


; ; z LE 
‘aot el ques odd to poip1sag) riod 


od yam (2vx2%%3 2) noitaups i 


re.iD 


ae 


> 
* 


: data) sabia - 


me 
3% s yd eusetosS bluow nt = ae 


; 7 =! — 
Bhucw , x Inetanod eset on 














4 


r a 5 7 a 
evisapaa Stereqgs a 7) 
ow n 
7 - 


few elt ‘sxsqmoo ‘yem on0 en 
i 

se yO 
aspossin ont te ‘gant od (ait ie ) 


a 
= 


3 i 7 


Fal ¥ - 


a. 


“ oe Z 
3) nolteupe Boek sexta | ta 


oe 


T 


eS 7 
a i 


aera 


hi 













ays? 


nowsbat ore 2 bas oo tk 


; 
ee . 
x «3 
STi 


4 * . Pe. ms ae mts = re 
sagen? ywoaefor®te noleatitoo seal besdetis ed SIwow 


‘apc eae 
sans sone alts eet a 


po 


226 


third body was neon (reaction 6:13) ke at 297°K was 
-29 2 2 | 
about 5.9x10 cc’ /molecule’ sec. Fehsenfeld et al measured 


baa aidbtecite sec: at 280°K for 


a value for k, of 1.9x10 
the reaction with He as the third body. One might expect the 
efficiency to decrease with decreasing complexity and size 
N.?NerHe. | 

Equations (6:xxxii) and (6:xxxv) may qualitatively show 
when the overall attachment reaction such as (6:1) or (6:2) 
displays third or second kinetic order. The case for third 
order was that ky77k [05]. — This would occur at low pressures, 
and at high temperatures for complex ions with a small bond 
energy. REARS the second order kinetic situation k 477k [05] 


would occur at high pressures and low temperatures and when 


the complex bond energy was large. 


ii The Effect of Temperature and Bond Energy 
on k.. at Low Pressures 
When the attachment reaction displays third order kinetics 
the rate constant ke is given by ko the constant for the activation 
of O a by collision. This activation is similar to the initial step 


4 


+* Ze + : : 
(6:17) 0, + 0, Seeger > O, + 0, (activation) 


of a unimolecular decomposition reaction and the same equations 


may be used. 


"i 
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The most appropriate equation is the energy transfer 
equation developed independently by Hinshelwood and Lewis. It 
is discussed by Ritchie (120) and by Fowler and Guggenheim (121). 
It is derived from the simple collision theory. The simple 


collision theory equation (6:xxxvii) assumes that the molecules 


(6:xxxvii) k & Ln uate des 


are hard rigid spheres, and that upon collision (along the line 
of centres), translational energy is transformed into internal 
energy. The probability of gaining an energy greater than Eo 


is given by the Boltzman factor ig ae 


- For a complex molecule 
the approximation of rigid spheres is not satisfactory and the 
contribution of vibrational and rotational energy must be 
included. The bond which requires the energy ES to break, 
gains this eharay not only by eraneeoraarion of translational 
kinetic energy into internal energy, but also by transfer of 
energy from the vibrational and rotational modes. The probability 
of E being available is thus larger than that given by the 

Ss 


E 
Boltzman factor. A term a 3 is added to account for the increase 


s'\ RT. 


in probability. The integer s depends upon the number of square 
terms from which energy is acquired. This number is given by (2s+2). 
The number of square terms for the degrees of freedom are, one 


2 ; 
term (4IW) for each rotational degree and two terms (nurs kinetic 
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energy, and iene potential energy) for each vibrational degree. 

Since the vibrational degrees of freedom will generally provide 

the largest number of square terms one might approximate s as 

the number of vibrations. The energy transfer equation for kK. 

is thus given by (6:xxxviii) in which.Z is again the collision 
s ~E /RT 


a a eel ee 
(6:xxxviii) Kae au (= ; 


number and ES the activation energy. 

The energy transfer equation was used successfully by 
Palmer and Hornig (122) to describe the dissociation of Br. into 
Br atoms in the presence of either bromine or argon under low 
pressure conditions. The log of the rate constant plotted versus 
1/T fitted over a range of foen, 

The value of Kk. measured for the oxygen reaction may be 
used to determine s. Z is taken as the Gioumousis-Stevenson 
collision frequency which is of the order 10 °cc/molecule sec. 

T£ we use T=300°K, Ey = 10 kcals/mole ands = 4, kK. is ToeKiOMeeCa/ 
molecule sec. This corresponds very well with the present 
experimental value k, = 2 e4dn ines cc/molecule sec. At 400°K 


k mento ocymolecite sec. compared to the experimental 


a 


k. = Peis ia eoenislacuic sec. 


- + 
If it is assumed that 0, and N, have the same number of 


vibrations contributing energy and thus the same value for s, 


Kk. may be calculated for nitrogen. The bond dissociation 
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energy E may be taken as 14 kcals/mole as the value of 11.5 
kcal/mole from Varney (107) had a large error, and the value of 
13.8 from Asundi et al (103) should be a minimum value. At 
300°K with S = 4 equation (6:xxxviil) gives a value of 
-15 z -14 

kK. = 9.06x10 cc/molecule sec. and at 400°K KS = 9.8x10 cc/ 
molecule sec. 

A value for s may be obtained from a method described by 


Porter (118). He first found an equation for k_ in the following 


manner. ke may be defined by equation (6:xxxix) which is (6:xi) 


(6 =xxx15) k = k K 


rewritten. Kk. is given by equation (6:xxxviii) for k.: The 


equilibrium constant K is given by (6:xl). The energy required 


(6:x1) x = ef AH/RT) , AS/R 
4% 


4 ) 


id : : : + 
by O (or N ) to dissociate ES is the bond energy D(O, -O 


4 2 
which is -AH. k, is then given by (6:xli) as the activation 


£ 
Sie s AS/R 
(6:x1i) k. = * (2) = 


s!\RT 
energy exponentials cancel. Porter assumed that the change of As 


AsS/R 


with temperature was small so that the effect of e was small 


A 
compared to (1/T)*. Zande e7R 


Equation (6:xli) may then be rewritten. 


E 
(6=x1i) k. = se (“3 
ey Wace 


Ss 
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are combined into a new constant C. 
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Porter stated that an Arrhenius plot (log ke versus 1/T) 
of equation (6:xlii) gave an apparent negative activation energy 


with the form given by equation (6:xliii). 


(G2:x12137) E = -SRT 
app 


The apparent negative activation energies obtained from the two 
Arrhenius plots, figures (6:13) and (6:22), were for oxygen 

Ey = -2.0 kcals/mole and for nitrogen Ey = -2.3 kcals/mole. If 
these values are substituted into (6:xliii), with T = 300°K 

the end of s for oxygen is 3.4 and for nitrogen 3.9. These 
values agree with the value of 4 used for the calculation of 
kK. by means of equation (6:xxxviii). 

The number of vibrational modes of four centre ions are 
given by 3N-5 = 7 for a non-linear ion or 3N-6 = 6 for a linear 
ion. s should be less than either of these numbers as it is the 
effective number of vibrational modes contributing energy to the 
dissociation and not all will contribute with total efficiency. 

The statement that equation (6:xliii) gives values of Ke 
which plot versus 1/T to give a negative activation energy was 
= a. Oxte eeeenateeule aad. and 


= 
-30 
Bm seehh.=+10 .kcals/moles” If-s.5.3, G.=y3. 76x10 Se a fe 


checked. At 300°K, k 


400°K a new value of ky was calculated Poem eenat Ainiecui es sec. 


From an extrapolation of figure (6:13) ke = Le27x100 2° eenoalacute 


sec. The calculation was repeated for s = 4. The agreement 
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with the experimental results was not as good. This shows 
that s can be approximated from the value of the apparent 


negative activation energy. 


iii Comparison of the Equilibrium Constants Kp 
for Oxygen and Nitrogen 
Since the equilibrium constant for the nitrogen reaction 
(6:2) could not be measured experimentally, it would be instructive 
to evaluate a theoretical value by means of the energy transfer 
theory. kK, can be calculated from equation (6:xxxviii) and used 


with the experimental value for k. so that Kp may be calculated. 


f 
Since the nitrogen reaction came closest to reaching 
equilibrium in the 380°K runs 400°K will be chosen for convenience 
of calculation. First the theoretical equilibrium constant for 
oxygen was calculated at 400°K with s = 4 and Eo = 10 kcals/mole 


=12 
and thus ke = 3.80x10 cc/molecule sec. The experimental values were 


k. = Be eelOur cc/molecule sec. and Ke = Gere ale take Ea iaiece sec. 


Thus the theoretical value for Kp is given by (6:xliv) from (6:xiii). 


Se 
(6:xliv) Ko => 


ies T1.034x10 7? 


Kp (theoretical) = AOxth ae compared to Kp (experimental) 
= 4-910. > and the difference is less than a factor of two. 
This procedure may be repeated for nitrogen where 


T = 400°K s = 4 and E = 14 kcals/mole. kK, = Ar okipiaertat ree sec. 
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M 
the theoretical value of Kp = 7.4x10 . The significance of this 


lue i i 28 © ee ion } 
value is seen by comparing *Pu5 Oo Po, he ration Py / Fo, ae 


400°K is: 


KPy 3 

2 9x10 (Theoretical) 

Kp — 
5 





If we take the experimental value of Kp. at 400°K to be > age 
2 


and if the experimental Kp = 4.7107 the ratio is then: 
: 2 


Kp 
‘2 





za 2x10°7 (Experimental) 
Kp, 
74 


The ratio may also be obtained by an extrapolation of Varney's 


(107) results. 


Kp 
N 
> 7x10> (Varney) 





Kp 
or 


Since the predicted value of Kp. using the experimental ky 
2 


and the theoretical ke calculated from the bond energy by means of 


the energy transfer theory equation (7:vii) is smaller than either 
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of the experimental values (even though they may be only approximate), 


+ + 
the greater stability of Ny compared to O, must be due not only 


to a larger bond energy but also to a more favourable entropy. 
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+ + 
6:5 The Attachment N + N. = N, 


i Previous Work 


, : Dae wae 
It was stated in section (6:3:iii) that N is produced 


3 
by two reactions (6:8) and (6:9). The first is important at 


low pressures and the second at high pressures. The formation 


“fb 
of N, was studied at pressures up to 0.2 torr by Saporoschenko (110) 


, + 
(6:8) N Ne ote N +N 


+ ~ 
: > 
(6:9) ING 2N., N, aF N, 


who observed that the reaction was second order and that the 


ob 
appearance potential for N, was 22.1 eV which corresponded to 


3 
+ (4 


the excited state N. ED. Thus N % was formed by reaction (6:8). 


3 
This mechanism was supported by observations of Cermak and 
Herman (123), Keller et aZ(102) and by Munson et aZ (115). 


+ 
Dreeskamp (124) observed N at higher pressures (0.8 torr) 


3 
and ascribed its source to reaction (6:9). This was also reported 
by Luhr (125). Knewstubb and Tickner (126) studied glow discharges 
in nitrogen at a pressure of 0.4 torr. They also concluded that 
Ne resulted from the third order reaction (6:9). McKnight et al 
(104) recently studied the mobilities of nitrogen ions in N, at 
pressures up to 0.9 torr and observed that Nn" and Ny had very 
similar mobilities. They calculated a rate constant for the 


= sien? NEN sec. for 


third body reaction (6:9) of 3x10 
conditions of low drift velocity. More recently Fehsenfeld et al 


(113) reported reaction (6:18) in which the third body is helium. 
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+ 
(6:18) nN” +4 N,+He + N, + He 


ii Purpose of This Study 
The reaction conditions in the present experiment were 
Similar to those used earlier by Dreeskamp, Luhr and Knewstubb 


and Tickner,who showed that N Je was formed from NY by the third 


3 


order reaction (6:9). It was shown earlier (section 6:3:iv) 


+ + + 
that the ions N- and N. reacted independently of No and Ny ° 


+ 
N is formed from N, by reaction (6:19) with an appearance 


(6:19) N, +e > Noe gee N 


+ 
potential of 24.2 eV. The intensity of N could be decreased 


: ae : 
by reaction with N. to form N, , Or by charge exchange reactions 


with impurities. The rate of appearance of impurity ions was 


: 3 ae 
too slow to account for the rapid rate of disappearance of N . 


4 


+ 
The rate of formation of N, agreed quite well and thus it 
+ * 
appeared that N- reacted only to form N, - The rate constant 


+ 
for the formation of N, was measured so that it could be 


Ey 


4 (reaction 6:lb). 


compared to that of N 


iii Results 


, : ; aa 4 
The time dependence of the ion intensities of N and nN, 


+ + 
had been measured with those of N. and Ny . The intensities 


7 + : ° 
of N- and Ny were normalized in two ways. Some were plotted 


as fractions of the total ionization. The plots for pure 


+ - 
an + u + '-6H ¢ 
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nitrogen at pressures from 0.96 to 3.60 torr at 380°K are 
shown in figures (6:14) to (6:17). Two plots for mixtures 
of nitrogen and neon at total pressures of 4.1 torr are 
shown in figures (6:26) and (6:27). The temperatures were 
297°K and 376°K respectively. One can see that in these 


. set Pe Ta: : 
figures the sum (i, + i ) remained reasonably constant. 
3 


ee ae aaeoy 
For some runs only N= and Ny intensities were measured. 


: ; ; + eet: 
They were normalized as fractions of (1, cE iv ). These 


3 
plots are shown in figures (6:28) to (6:30). 


For most of the data the rate constant for the forward 
; ; ; ~ + 
reaction was calculated from the decay of i by means of 


equation (6:xlv) where n=l for a second order reaction and 
toot ; + 
| 2.508 POP Cy | ys “Ny, ) 
(6:xlv) k_= 2 1 


n 
(t.-t,) [NJ] 


n=2 for a third order reaction. For the data obtained from the 
N.-Ne mixtures (figures 6:26 and 6:27) the increase of in . was 
3 
aa 
used as the decrease of iy was rather unreliable. The rate 


constants are shown in table (6:4). Plots of Ke (n=2) and 


RE 3) (second and third order constants) are shown in figure 
(6:31). The plot indicates that the reaction is second order 
since the second order values are constant with pressure, and 


the third order values form a hyperbolic curve. The overall 


reaction should then be written as second order as .(6:9a). 


(6:9a) NN ee: EON 
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The second order rate constant appears to be temperature 
independent. The reactions did not appear to reach any 
equilibrium in the short time that they could be observed 


(O to 60 micro seconds). 


TABLE 6:4 


se 


Rate Constants for the Reaction 


+ + 
+ > - 
N nN, N, + (n-])N, 


N. Pressure 2nd order k 3rd order k> 
Temperature (torr) cc/molecule Bee oe 
Eye n= bree ae 
300°K 0.63 Leal 8.84 
0.83 i Aral 6.38 from 
120 2.06 6.39 
1.28 ix. 86 Aeb1 
300°K 0.61 Te24 “eye 
0.79 Dione ok A.77 Grom 
1.06 I-43 4.19 
378°K 0..50 nai 5.49 
1.02 0.68 2263.57 om 
Lene 1.06 ey 
380°K Jam 9 be 2.10 4.04 
from 
2.05 1.41 2.09 
3.4 1.64 1.89 
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n=2 
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iv Discussion of the Results 
: a 5 
A mechanism for the attachment of N, to N may be written 


exactly like that for attachment of N. to N."- 


2 
k 
+ c (combination) 
: uae +% 
(6:19) N + N, . N, 
da (decomposition) 
k 
+% s 4 (stabilization) 
° —_> : 
(6:20) N, + N. : N, + N, 
a (activation) 


The same conditions for order apply here. The overall forward 


reaction was second order when k3S<k, [NL]. 


The experimental rate constants are given by (6:xxx) and 





(6:xxxi). 
(6: xxx) Ke = Ko 
(6: xXxxi) Kk. = Kk. Ks 
< 
s 


ka and Ke were defined earlier by (6:xxiv) and (6:xxxv). 


r-1 
Z RT 
(6:XxXxXiv) ae = d ae 
nRT 
(6: xxxv) Rema pe (lee co) 


The conditions that k <<k_ [NJ] are that the pressure be high, 


d 


the temperature be low and the bond energy of the complex be 


high. The pressure and temperature under which the reactions 
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(6:2) and (6:9) were observed in the present study were identical. 
The bond energies and the number of vibrations are different. The 
bond energy D(N,-N,) was determined from appearance potential 
results of Asundi et al (103) be greater than 0.6 eV (13.8 kcals/ 
mole). These workers also measured the appearance potential of nt 


; ot ° 
(from reaction 6:19) to be 24.2+.1 eV and that of N, (from reaction 


6:8) to be 21.1+.1 eV. The difference between these (3.14.2 eV) is 


the minimum energy for the bond N’-N .- Thus D(N’-N,) ~ 71 kcals/ 


2 


: tity : 
mole. Since the bond energy of the complex N, is five times 


+ + 
one would expect that k, for N would be 


larger than that of N d 3 


4 


+ 
much smaller than kg for N, and this is probably the reason that 


the reaction (6:9) appeared to be second order in the present 


+ 
experiment. The greater bond energy of N is probably due to 


3 


: + : 
its structure. Ny has an even number of electrons in the 


-f 
valence shell whereas N, has an odd number. Simple Lewis 


+ + 
structures for Ny and N, are shown below. 
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The lack of temperature dependence for Ke (second order) 
agrees with the theory presented earlier as ke = koe and ky 
was shown to be the Gioumousis-Stevenson constant which is 
temperature independent. 

: -29 2 2 

The value of McKnight et aZ(104) of 3x10 cc /molecule 
sec. measured at 300°K and at 0.9 torr No can be converted 
: : : -12 
into a second order constant with the value of 0.87x10 Gos 
molecule sec. This is very close to the average of the 

-12 
present results, 1.5x10 - cc/molecule sec. McKnight et al 
suggested that a reverse reaction took place under drift 
field conditions at high ion energies. This was not observed. 
One would not expect it to be, as the reverse reaction would 
+ 

be smaller than that for Ny which also could not be 


observed. 
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SUGGESTIONS FOR FURTHER RESEARCH 


TL) introduction 


The most abundant molecule in the Earth's atmosphere after 


N, and 0. is Ho 


O which is present in trace quantities. 


Attachment 


reactions between these molecules and their ions are expected 


to occur in the D region of the ionosphere where the pressure is 


high. 


produced by ultraviolet and x-ray radiation, and by charge 


ob 
exchange from N. 


Cis8) 


7:2 Attachment R 


The first suggestion would be to study reactions in oxygen 


containing traces 


; . ; ' ; + : y 
One of the major ions in the D-region is 0. which is 


(127) by reaction (7:1). 


eactions in O,-H,0 Mixtures 


of HO. 


While performing the experiments 


described in chapter 6 the presence of two masses corresponding 


t+ 
to the hydrates 0," (H,0) and 0. (H,9) 5 were observed when a 
-3 
trace of water, about 1x10 Lorzper, torr,.© 


could therefore propose the following series of reactions: 


(6:1) 
(7:2) 


(7:3) 


2 


Toe a, EN eh Sea 
ig ea att im a 2 2 


0.°(H,0)_ +H,0 ——> 0,°(H,0)_,, +0 
aoe oe eee On aed g 


was present. 


One 
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One would expect the reactions to require third body stabilization, 
and since one would expect the water to be in only trace concen- 
trations 0. would be the third body. 

This instrument could be used to measure the rate constants 
for the forward (and reverse) reactions, to arp the time 
required for equilibrium, and to measure the equilibrium constant 
Kp. From temperature studies the activation energy Ey of the 
forward reaction and the enthalpy AH could be determined. It 
would be interesting to see if the reactions had small negative 
activation energies of 1 to 2 kcal/mole as did the attachment 
reactions of oxygen and nitrogen (reactions 6:1 and 6:2), and 
to see if there was any trend in E; as the cluster size 
increased. 

The bond energies (-AH) of other clusters H (H,0)_, 
K"(H,0)_, and NH,” (H,0) have been previously measured (31). 

The bond energies of hydrates of other ions me Na’, Rb’, 

Cs", are presently being measured(128). It would be interesting 
to see how the stability of the oxygen hydrates 0,” (HO), 
compared with these hydrates of atomic ions. The bonding of 

the hydrates with atomic centres has been explained as 
attraction of dipoles to a spherically symmetrical charge, 

and the stability of the ions found to change systematically 


+ 
with the charge radius (31). The molecular ion 0, would not 


have a spherical charge centre. One might be able to modify the 


cbs - 
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theory discussed by Searles (31)to fit this situation. Finally 


hydrates of O ~ have been observed (129). A comparison of 


2 


the positive versus the negative ions could be made. 


7:3 Preliminary Investigation of 0,-H,0 Mixtures 

The investigation of the attachment of HAO and ot was 
undertaken by Dr. Good and myself while this thesis was being 
written. Since it was not part of the thesis project, only a 


brief description will be given. 


i Experimental 

A flow system was set up whereby oxygen was admitted 
through the variable leak into the manifold. The gas passed 
through the ion source, and then through a capillary into a 
mechanical vacuum pump. At the atmospheric side of the pump 
the gas passed into a flow meter. Water was admitted to the 
oxygen between the variable leak and the manifold, through a 


small capillary set into the gas flow line. The total pressure 


of oxygen and water was read off the M.C.T. manometer. The ratio 


of oxygen to water was calculated from the weight of water added 


to the oxygen, from the total pressure and from the oxygen flow 


rate. 


ii Results 


, + 
Many more ions than expected were observed. Besides 0. ’ 


+ + + + 
Oo; P 0. (HO) and 0, (HO) 5. the water ions HO : H, 


+ + 
O,H (HO) ,, 
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+ 
H (H,0), etc. were observed. Figure (7:1) shows the data for 


a mixture of 3.x10-° torr H,O in 2,2 torr O, at 308°K. A 


kinetic system much more complex than proposed earlier was 


suggested by the data, this is shown below. 


(6: 1b) Sek, WaT ea 

(7:4) 0, +H,0 + 0,° H,0 +0, 
C725) Qu HO + HO * H,0” + OH + oO. 
(726) ae HO + HO +0, + oy (H,0) + 0, 
(7:7) HO” + HO +0, > a0 (H,0), + 0, 
(7:8) oO,” (H,0) + H,0 > H™ (H,0) + 0, + OH 


~ + 
: —>- 
(7:9) H (H,0), +/H,0™. oO, H  (H,0), + 0, 


+ + 
. ==> 
(7:10) H (H,9) , re HOE H (H.0) , +20 


Za 2 2 


It appeared that under the present conditions the reaction 


: + 
(7:2) had too slow a rate to compete with the formation of 0, by 


+ + 
J > 
(7:2) oO, ao HO + 0, 0. HO = oO, 


reaction (6:lb). The concentration of 0. was about 100 times that 


of HO. Thus the rate constant for (7:2) must be smaller than 


eae ba aoa lecule. sec. as ke for (6:lb) was measured as 


30 


3x10 ernvecuica sec. 


2 - 
Another observation was the transformation of 0. (H,0) 
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Figure 7:1 Normalized ion intensities for 0,-H,0 mixture at 308°K, 
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into eee It appears that the hydrated protons are much more 
stable than the hydrated on ions. One would expect the 
electrophilic nature of HY to be much larger than that of 00h 
This may explain Perea ce for the transformation. 

The formation of H (H,0) via reactions (7:5) and (7:7) 
which display second order kinetics appeared to be much larger 
than via reactions (7:6) and (7:8) which would display third 
order kinetics. Shahin (130) also observed the transformation 
but ascribed its formation to reactions (7:6) and (7:8). 
However the pressure in his apparatus was very high (30 torr), 
about 10 times that in the present apparatus and thus the 
third order mechanism might be expected to predominate. 


The rate constants for the successive water hydration 


reactions were measured. They are shown in Table 7:1. 


TABLE 7:1 


Rate Constants for the Forward Reactions for Hydration Reactions 


H’ (HO) + =O) +0 peer (EG) + O 
PDE S | 2 2 2 ntl 2 
: a2 2 
Reaction n ke cc /molecule' sec. 
3 f WESSON ae 
7:9 2 tear 
-27 
7:10 3 Lox 10 
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A previous value for the rate constant ke for the 


‘ ee : : 
hydration of HO (reaction 7:7) but with N, as the third body 


28 


has been previously reported as 7x10 ce’ fmolecule” sec, (131). No 


details of the experimental conditions were given but one would 


expect that the presence of N, instead of 0, as the third body 


would not affect the rate constant. 


. 


Since the reaction has so far been studied at only one 
temperature the parameters discussed earlier Eee Eo AH, and Kp 


have yet to be measured. 


7:4 Attachment Reactions in N,-H,0 Mixtures 


A second set of reactions would involve reactions in No 


containing traces of HAO. One would not expect hydrates of 

the form N,” (H,0) to be observed however, as the ionization 
potential of No” 15.58 eV is much larger than that of HO 

(12.61) and if one expects such a hydrate to have a recombination 
energy less than 12.61 to remain stable, 3 eV excitation energy 
would have to be removed by collision with third bodies. Rather 


: + : 
one would expect charge exchange to occur to give H,0 im lois 


system has been studied somewhat previously. Shahin has observed 


the ions from a corona flow of N,-H,0 mixtures and measured 


some reaction cross sections (132). He proposed a reaction sequence 


+ + * 
which is shown below. We also observed ions HO ; HO - NOH 


+ “ aa 
N and N when traces of 


+ a 
and H (HA) 5 as well asN , N, rN, A 
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HAO were present in the N.- 
+ 

6:2 

( ) N, + 2N., 

Gray n.+ + HO 
; 2 2 

+. 

7:12 + 

( ) N, + H0 

Bey Nett HeO 
; 4 2 

(7:14) No aes 
: 4 2 

Cr. 15) NH’ + HO 
; 2 2 

(7:16) H.O' + H.O 
: 2 2 

ea ieeh eh O-tAN 
: 3 2 2 

-f- 
(7:18) H’ (H,0)+H,0+ N 
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OH 


2 


a a 
a 
H (H,0) , + No 


Shahin (132) measured the cross sections for reactions (6:2 


and 7:11,12,15) but did not calculate rate constants. 


Further studies on this system are warranted. Accurate 


values for the rate constants instead of cross sections (obtained 


by Shahin by making many approximations) might be obtained by 


means of the pulsed mass spectrometer in spite of the apparent 


complexity of the reaction system. 


Neither rate constant nor 


, : . : 7 
cross section have been measured for the reactions involving Ny : 


(7:13) 


(7:14) 


N, 7 Ho 


Ny + Ho 
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The irradiation of N,-H,0 mixtures might be more useful 
than irradiation of 0,-HA0 PE ay: to study the kinetics of 
the hydrates of water as no other hydrates e.g. 0," (H,0) are 
formed. The discrepancy between the rate constant measured 
in the present instrument for reaction (7:7) and that reported 
for (7:17) might be resolved. The forward reaction, combined 
with the Kp previously obtained by Searles (31) and Kebarle (133) 
would give the rate constants for decomposition of the hydrates. 
It would be interesting to establish the order of the forward 


reactions in the presence of N Searles (31) stated that one 


2° 
might expect that, as the cluster size increased, the reactions 
would not need stabilization since the increased number of 
vibrational and rotational degrees of freedom could absorb 

the excess energy due to bond formation and therefore would go 
from third order to second order. 


The prediction of the equations for k, from Section (6:4:ii) 


£ 


is not the same. The third order rate, occurs when k 477k, [0]. kg 


and Ke are given by (6:xxxiv) and (6:xxxv). 


r-l 
nRT 
(6:xxxiv) ; ka = Z 4 [ner] r = number of 
vibrational modes 


pz. eieaethoy 


(6: xxxv) ko 


As the cluster size increases, the bond energy decreases, which 


would increase Kg and thus decrease k However the effective 


f° 


number of vibrational modes r increase also and since nRT/(D+nRT)<1l 
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this would decrease Ka and increase Kye The addition of one 


HO molecule to the hydrate increases the total number of 


vibrational modes. If r is always proportional to this number 
(3N-6) a change in N would have a greater effect than changes 


in D and thus 36 would decrease. The measured rate constants 


k, for further reaction of the water hydrates H (HO) m0 


£ 2 


decrease which indicates that changes in D have larger effects 
than changes in r through addition of more HO clusters. Thus 
the fraction of the effective number of vibrational modes 
contributing to the decomposition of the excited state, out 
of the total number of vibrational modes decreases. Thus the 
third body requirement would increase. This is in contradiction 
with Searles statement. Further measurements of the NO 
system might clarify this. The ideal situation would be to 
find a pressure-temperature region in which the addition of 
the first HO molecule to H,0” were second order and then 
test whether the addition of further H,0 molecules to the 
cluster were second or third order. 

The measurement of these reaction rates and their temperature 


dependencies would be useful to help explain the processes occurring 


in the D region of the ionosphere. 


+ + 


7:5 Further Work on the Attachment N, + N, a Ny 


Improved values for Kp and AH and thus the bond energy 


D(N, -N,) should be obtained by studying the attachment reaction (6:2) 
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at temperatures from 500 to 600°K. These could also be used 
to correlate Varney's ion temperature (107) with the thermal 
temperature. The ion source of the present instrument would 
have to be modified so that it could operate at such high 
temperatures as the present limit is 450°K. 
7:6 The Attachment Reaction co* 12 So) ee ed eee 5 

Ss salar male a celia > ares Ry 

The attachment of CO to co” Hes been observed by 
Saporoschenko (134). Since co" has a bond length and vibrational 


+ 
frequency intermediate to those of O and N,” (135) one might 


2 
by the study of the formation of acc be able to correlate 


all three attachment complexes. Recently Conway (136) has 


oe ; : 
O which one might expect to have a 


studied the complex No 3 


+ : +. 
similar structure to C20, - The formation of C50, in, CD 


would probably be third order (reaction 7:1la). 


+ 
(71a) cot + 2c0 > C0," + CO 


+ 
Saporoschenko estimated the bond energy D(CO -CO) to be Vv 0.81 eV. 


It can also be estimated from the appearance potential of C,0, 
produced by the Hornbeck-Molnar reaction (7:20) (12.840.3 eV). 


TV 
The ionization potential of CO is 14.00 eV. Thus the bond 


* + 
(7:20) CO” CO C,°. +e 


energy should be about 1.240.3 eV. Since this value is large, 
3 
a van't Hoff plot of log Kp versus 10 /T would have a large 


slope. It might be necessary to operate at high temperatures. 


254 


v6 18.0 ¥ 


i 


a 


~ 


4 
a” 


< 2 
. (We £828. RO) eet) Mpktoses santon-sose 


pele BE ‘ely emt *068 02 00 ott 2a! 
‘le Shy oe "O08 . ork ot 7 


“et Of (CO- "oo)G ¥ sor giye © 




















ay 
; 8 ; : she 
; toge (DL) exitersg angst ‘aod! ewan 


‘ ar 
> au 
te « 2avo8 ee: vt _ : 


~ iP . 


i tege dread of 


i? fi 
ened ai jimif saseesg oad a 
i; ; = it 
" Py =) aly 
‘ig «co + “oo nobgudet Snem jaumioaasa oct 2 ant 
, : es 4 es ae aa nf 
‘ 
od Bevesedo mead 21 OD OF an to snare ott 
. : i = re . “) 
, "a ’ se og 
cot gael Sood s ast OO som (PEL) nteae 
, + y aL 7. 
c ang (8E2) * ft bas 3 t¢ oeors of ataideansgat yeonsy 
4 al aed - = ¥ 
rsfetsoo oF aida od "0,9 io aes jersos a to youse & 
7 ae 
ers wmseun yt xo fghor snemivstas % 
J (Obs? y ~~ ae 
‘ ‘ ‘ "2 + 
6 wat oF Boar feist SR 4a Lflw oOo" x Lame edd Be. 
/ rf # 
o>: wt ga,.0 2e sot : off «. ,0.9 09 pustusg® 
. ¢ $ » % " i 
fal:< aoisoses! tebto bald od idee ig bis Ww 
» “4 
| wei ZA 
4 : - - 
e+ 70.9. + OE + ee a 
> i a, ' i 
-) i “si ; 


band asd Sovamiseo © ny ruc 
“ecqgs o word apr 


E 7 
fakensgoy soastecucs 
ated 


al 





Rudolph and Lind (137) and Munson et al (115) postulated that 


the following reactions produced the observed ions CO bi and 


C0 . 
C722 1) 
(7:22) 


(7 3.23) 


(7:24) 


The presence of these reactions would complicate the observation 


of the kinetics of reaction (7:la). 


co’ + co” 
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© 64+5Co 
+ + 

CO + CO 


oe 
Gets CO 
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APPENDIX I 
Equation for the 


Correct Value for the Forward Rate Constant for a Reversible Reaction 


Consider the reversible reaction (1) where n=1 or 2, for 


= “f = 
(1) vs + nX eee xX, 4 (n-1) xX, 


k 
xr 


a second order or third order reaction. If we assume that the 
ion concentration is proportional to the measured fractional 
intensity we may write the change in intensity (i): 

(i) rs a -k A [X.] + kK. xe] 


Xo = £ X, 2 r X, 2 


dt 





at equilibrium diy = O. k,. and k_ are related to the equilibrium 
eee 3 


£ 
constant, and dt to the ion concentrations at equilibrium. 
+ 
(ii) kee kei iy eq 
+ 
kK. i eq [X,] 
X 2 
2 
Thus 
ey ae 
(iii) kL = ke [X_]¢ where C = ae (eq) 
+ 
i, (eq) 
4 


+ PME. ¢ : 
Since the two intensities iy and i, are fractions of the 
2 4 


total ion intensity we have (iv) 


(iv) a3 = l-i 
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we may substitute (iii) and (iv) into (i). 


. + n ' + 
(v) ae se -k, [0.] BT Clo oy 





; ‘ ae 
This may be rearranged and integrated from o tot. Leti = xX 


Xo 
for convenience. 
me t 
(vi) dx z =k, [x,1" at 
x(1+C)-C 
fe) fe) 
D4 
ee — ah n _ 
(vii) Qn [x (1+C) - C] = ke [x] (t ty) 
1+C 
ers 
re) 
at t=0, xo = land att=t, x= i, We The final equation is 
2 


then written: 


(viii) =e 2.303 10944 bs (1+C) -C] 
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X2 





(i+C) 
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APPENDIX II 


Ionization Potentials 





Ion Source Ionization Potential Reference 
+ 
He He +e 24.58 15 
+ * 
He. He + He 23013 Ap. 138 
ob 
Ne Ne +e 21. 56 15 
+ * 
Ne, Ne + Ne 20.86 Ap. 138 
ar’ “Ps, Ar +e 15.76 15 
2 
20, Ar+e aie PA i 1 Wes* 
aA 
2 
+ * 
AL, Ar + Ar 14.71 Ap. Wie BS, 
Kr’ 2p | Kr +e 14.00 LS 
3/, 
2p, Kr™+ e 14.67 Lo 
He: 
+ * 
Kr, Kr™ + Kr 13.00 Ap. 139 
Yet “p Yale whe! 15 
3/, 
20. Xe +e 13.44 15 
+ * 
Xe, Xe + Xe i1.16 Ap. BA SS 
nt N +e 14.54 138 
nt Rote 24.2 103 
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APPENDIX II (Cont'd) 


Source 


CH +e 
C,H, +e 


a 
CoH, a 


Ap. = appearance potential. 


Ionization Potential 


aod 


Gee 


el keaal 


fi. © 


ems RE, 


Leo 7 


16.2 


a2 661 


10.51 


13.93 


Eka 


Ap. 


Ap. 


Ap. 


Ap. 


Ap. 


Ap. 


Reference 
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